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PREFACE 


There  exists  a  massive  body  of  literature  on  aircraft  wake  vortex 
research.  As  an  aid  for  further  work,  a  new  bibliography  has  been  pre¬ 
pared.  The  base  of  this  list  has  been  taken  from  the  excellent  listing 
given  by  McCormick,  1971:  Aircraft  Wakes:  A  Survey  of  the  Problem 
(FAA  Symposium  on  Turbulence,  Washington,  D.  C.  ).  This  base  has 
been  complemented  and  updated  by  AeroVironment  and  is  believed  to  be 
very  comprehensive.  The  actual  report  refers  to  only  a  portion  of  these 
papers,  and  contams  other  references  which  are  not  directly  concerned 
with  vortex  wakes.  Thus,  for  ease  of  reading,  a  separate  reference  list 
has  been  prepared,  containing  only  those  papers  quoted  in  this  report. 

It  is  pointed  out  that  this  is  a  first  report,  and  by  the  natui  3 
of  this  particular  subject,  contains  hypotheses  and  suggestions  for 
which  conclusive  experimental  data  is  still  lacking.  A  second  report  will 
appear  later  in  1973  by  which  time  more  analytical  and  experimental 
data  will  be  available.  It  is  believed  that  this  new  :esearch  should  settle 
some  of  the  crucial  questions  identified  in  this  report. 
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1.  INTRODUCTION 


The  hazards  posed  to  following  aircraft  by  the  vortex  wake  system 
of  a  large  airplane  are  well-known.  In  order  quantitatively  to  assess 
the  hazard,  as  an  aid  in  designing  wake  avoidance  systems,  it  is  evidently 
necessary  to  know  where  the  wake  is,  what  its  strength  is,  and  what  the 
danger  to  an  encounter  aircraft  is.  These  three  aspects  are  addressed  in 
this  report. 

The  hazard  determination  depends  upon  the  strength  of  the  wake  and 
the  response /control  characteristics  of  the  encounter  aircraft.  For  a 
given  wake  this  becomes  a  problem  in  aircraft  dynamics,  and  the  specialized 
aspects  of  this  are  not  treated  here,  but  in  Section  3  simple  danger  indices 
are  derived.  The  technology  is  ,~vailable  for  the  refinement  of  these  in¬ 
dices,  but  in  view'  of  the  uncertainties  of  wake  strength  and  position,  such 
extra  precision  may  not  be  required. 

The  main  thrust  of  this  report  is  directed  at  the  determination  of 
wake  strength  and  position.  These  characteristics  are  fundamentally  con¬ 
nected  with  the  initial  parameters  of  the  wake.  Section  2  describes  the 
mechanism  of  vortex  sheet  generation  and  roll-up.  It  is  believed  that  this 
aspect,  essentially  an  aerodynamic  problem,  is  reasonably  well  under¬ 
stood  and  quantified. 

After  roll-up  ,  we  are  concerned  with  the  descent  and  decay  of  the 
wake  in  a  real  atmosphere.  This  atmosphere  contains  omnipresent  turbu¬ 
lence,  buoyancy,  wind  shears  and  wind  gusts.  Effects  of  these  dynamic 
characteristics  of  the  atmosphere  on  the  vortex  wake  are  very  poorly 
understood  and  experimental  data  is  limited.  The  major  portion  of  this 
report,  Sections  4  and  5,  addresses  itself  to  the  principal  elements  of 
of  this  interaction,  wdlh  the  goal  of  arriving  at  quantitative  expressions 
for  these  effects  so  that  operational  predictive  equations  can  be  derived. 
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Finally,  to  specify  the  atmospheric  inputs,  it  is  necessary  to  be 
able  to  infer  or  predict  the  turbulence,  shear  and  other  atmospheric 
dynamic  properties  from  the  basic  meteorological  observables.  Section 
4  shows  how  this  may  be  done  by  using  various  earth  boundary  layer 
models  and  the  large  body  of  statistical  meteorological  data. 


2.  LIFT  AND  VORTEX  GENERATION 


The  process  of  lift  generation  and  vortex  shedding  on  a  continu¬ 
ous  wing  is  well  understood.  Thwaites  (I960)  gives  an  excellent  re¬ 
view  of  modern  win^  theory.  We  note  that  the  determination  of  the 
spanwise  and  chordwise  loading  on  a  wing  is  not  of  concern  here,  but 
that  a  vast  volume  of  literature  treating  this  subject  exists  a-'d  m  d- 
ern  lifting  surface  theory  (as  employed  by  all  the  major  ai  .raf:  com¬ 
panies)  gives  a  prediction  of  the  wing  characteristics  ar.d  )  erformatices 
which  correlates  very  well  with  experimental  tests.  For  cur  pi  rpo»j 
we  are  concerned  with  the  way  in  which  the  trailing  vortex  core  system 
is  connected  to  the  bound  vorticity  on  the  wing.  For  the  moment,  we 
will  discuss  a  wing  in  attached  flow  with  a  continuous  trailing  edge;  that 
is,  where  there  are  no  special  disturbances  due  to  flaps  or  propulsion 
or  fuselage  attachments. 

The  vorticity  on  the  wing  itself  is  principally  of  i(.mw'se  orien¬ 
tation  and  associated  with  the  wing  boundary  layer.  If  the  wing  is  c.t 
zero  lift  this  vorticity  is  shed  oriented  approximately  parallel  to  the 
trailing  edge  and  results  in  a  low  energy  planar  wake  in  which  the  vel¬ 
ocity  perturbations  are  mainly  in  the  direction  of  the  free  stream,  and 
represent  a  reduction  in  the  free  stream  velocity  associated  with  the 
viscous  drag  of  the  wing.  The  momentum  loss  in  this  wake  can  be  di¬ 
rectly  connected  tolbat  in  the  boundary  layer  and  is  represented  in  an 
integrated  form  in  the  profile  drag  coefficient  of  the  wing,  Cjpo  .  The 
downstream  development  in  the  non-lifting  case  consists  of  entrain¬ 
ment  of  outer  flow  air  and  a  general  coalescence  of  the  wake  cross- 
section  from  a  sheet  to  an  ellipse  and  finally  a  circle.  However, 
these  details  are  not  important  since  we  are  interested  in  the  lifting 
case,  when  the  downstream  development  is  dominated  by  the  lift  in¬ 
duced  effects. 

If  the  wing  is  lifting,  an  additional  vorticity  component  is  gener¬ 
ated  at  the  trailing  edge.  This  vorticity  is  shed  parallel  to  the  local 
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velocity  vtTtor,  that  is,  approximately  in  the  direction  of  the  free 
stream.  It  can  be  shown  that  the  vorticity  shed  at  each  spanwise  sta¬ 
tion  is  exactly  equal  to  the  gradient  in  spanwise  lift  at  that  station.  It 
is  this  vorticity  which  rolls  up  to  form  the  trailing  vortex  pair  ard  with 
which  the  induced  drag  D.  is  associated.  This  vorticity  creates  (or  is 
created  by)  the  sidewash  and  downwash  fields  downstream  of  the  wing 
and  does  not  directly  produce  streamwise  flows;  however,  as  is  shown 
later, viscous  effects  control  the  rollup  of  this  vortex  system  and  ulti¬ 
mately  cause  streamwise  velocity  perturbations  which  couple  with  the 
profile  drag  vorticity  to  create  axial  flows  in  the  vortex  core. 

Continuing  the  analysis  of  classical  wing  theory,  a  number  of 
significant  (and  exact)  results  can  be  developed  for  the  relationships 
between  spanwise  loading,  induced  drag,  and  the  final  vortex  position. 

We  write  the  spanwise  circulation  loading  as 

T  =  2  Ub  ?  A  Sin  n6  ,  n  odd 
1  n 

where  cos  9  -  r,  is  the  normalized  span  coordinate,  U  is  the  flight 
speed,  and  b  the  wing  span.  For  symmetrical  loadings  only  odd  n  is 
required.  Then  the  center  section  circulation  Po  is  given  by 

n- 1 

r  =  2  Ub?  (-1)  2  A 
o  ! 

while  the  lift  L  and  induced  drag  are  given  by  the  formulas 

L  =  q  b?n  Ax 
D.  =  q  l^TT?n  A* 

i  i  n 
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y 

where  q  is  the  dynamic  pressure  l/2p  U  .  A  significant  term  in 
the  downstream  core  development  is  the  vortex  span  b  .  This  can 
be  shown  to  be  equal  to  the  centroid  of  the  shed  vorticity  (for  conserv¬ 
ation  of  lift),  and  thus  we  obtain 

.  n_b  _  ..  4  .  L 

^v  "  4  n- 1  ~  rr  pUT 

EC-1)  2  An 

These  simple  results  are  usually  derived  from  lifting  line  theory, 
but  it  is  worth  noting  that  they  are  exact  for  all  planar  wings  (even  of 
low  aspect  ratio)  providing  that  vorticity  is  shad  only  from  the  trailing 
edge. 


It  is  usual  to  consider  only  elhptically  loaded  wings  implying 
=  0  „  n>  1  .  We  remark  that  for  any  continuous  loading  the  additional 
coefficients  may  be  found  by  Fourier  analysis  and  the  more  precise  re¬ 
sults  determined.  However,  deviations  from  the  elliptical  loading  for 
aircraft  wings  in  the  clean  condition  are  usually  small.  As  an  illustra¬ 
tion  we  quote  results  from  Thwaites  for  a  number  of  different  untwisted 
planforms.  We  exhibit  these  as  factors  for  comparison  with  the  ellipti¬ 
cal  results.  Thus  for  an  untwisted  elliptical  wing  we  get 


D. 


L3 

q  rr  b8 


b 


v 


rrb 
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while  for  the  other  wings  we  write 


D. 

i 


K 


D 


L3 

qrr  b3 


b 

v 


I< 

v 


TT  b 
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Typical  extreme  cases  are  listed  below  for  a  wing  of  aspect  ratio  6. 
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Straight  Constant  Chord  Wing  1.05 

Straight  Pointed  Wing  (taper  ratio  1)  1.  14 

45°  Swept  Constant  Chord  Wing  ~  1.08 

We  note  that  these  cases  give  quite  different  spanwise  vorticity 
loadings.  Compared  with  the  elliptical,  the  pointea  wing  has  increased 
inboard  loadi  ig.  The  swept  wing  has  inboard  ic-ading  significantly  lower 
than  that  at  abou .  LOfo  span  (the  so-called  saddleback  loading),  while  the 
rectangular  wing  has  relatively  high  tip  loading.  However,  even  in  these 
extreme  cases  (within  norma?  wing  proportions)  there  is  not  very  much 
difference  in  the  final  vorte  .  span  or  the  induced  drag. 

It  should  also  be  recalled  that  most  aircraft  wings  are  aerodynam¬ 
ic;  ’> /  t  .isted  so  that  the  distribution  at  cruise  will  be  nearly  elliptical 
thus  they  will  not  exhibit  such  extreme  loading  distributions  as  the 
untwisted  planforms  noted  above. 

On  these  grounds  we  note  that  the  elliptical  assumption  is  certainly 
adequate  for  most  practical  wings  for  our  purposes,  especially  in  view 
of  the  much  cruder  assumptions  inherent  in  analysis  of  the  vortex  wake 
development.  However,  if  required,  the  precise  determination  of  the 
global  vortex  properties  can  readily  be  made  by  the  formulae  above. 

When  the  wing  is  in  the  dirty  condition,  with  flaps  extended,  the 
remarks  about,  elliptical  loading  no  longer  apply.  However,  the  three 
equations  for  L,  and  bv  still  obtain  and  the  calculations  can  still 
be  conducted.  With  part  span  flaps  the  circulation  distribution  contains 
gradients  which  are  logarithmically  singular  at  the  discontinuity,  which 
means  that  many  terms  in  the  sine  expansion  must  be  used  for  an  ac¬ 
curate  determination  of  D.  and  bv  .  This  can  be  avoided  by  introducing 
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an  aaditional  circulation  function  capable  of  exactly  matching  the  dis¬ 
continuities  in  airfoil  effective  angle  at  the  flap  extremities.  Using 
this  stratagem  a  very  accurate  solution  can  be  obtained  with  only  a 
few  terms.  A  modern  development  of  this  technique  applicable  to 
chord  and  twist  discontinuities  is  given  by  Lissaman  (19? 3).  This 
method  has  been  programmed  and  is  available  at  NASA  (Langley). 

It  represents  an  exact  solution  for  the  problem  discussed  here. 

The  effect  of  inboard  flaps  is  greatly  to  increase  k^  and  to  re¬ 
duce  kv  .  These  changes  are  such  that  the  elliptical  estimate  is  cer¬ 
tainly  no  longer  valid,  and  each  flap  type  circulation  loading  should  be 
analyzed  as  described  above. 


As  an  illustration  of  this  effect,  take  a  wing  of  aspect  ratio  5  with 
30%  chord  half- span  inboard  flaps  deflected  at  8°.  Taking  the  lift  distri¬ 
bution  from  a  similar  arrangement  given  by  Bilanin  and  Widnall  (1973) 
we  obtain:  . 


Straight  Constant  Chord  Flapped  Wing  ^  ^ 
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The  circulation  loading  distribution  of  a  wing  defines  all  the  pro- 
Arties  of  the  inviscid  vortex  development.  However,  as  pointed  out 
above,  for  most  vortex  development  purposes  only  a  few  characteristics 
of  this  distribution  are  required;  the  zero  moment  for  the  lift,  and  the  cen 
ter  section  bound  vorticity  for  the  centroid  of  lift.  In  addition,  we  re¬ 
quire  the  induced  drag  which  theoretically  involves  all  moments  and  is 
actually  quite  strongly  dependent  on  higher  terms  in  the  sine  expansion 

One  further  moment  will  prove  useful;  the  moment  of  momentum, 
that  is.  the  angular  momentum  of  the  fluid  rotating  around  the  vortex 
sheet.  It  has  been  shown  by  Betz  (1933)  that  this  can  be  conveniently 
evaluated  by  taking  the  second  moment  of  the  shed  vorticity  about  any 
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spanwise  position.  A  convenient  position  is  the  centroid  of  this  vor- 
ticity.  This  term  can  also  be  expressed  compactly  in  the  Four¬ 

ier  coefficients  as 
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Physically,  Mm  has  the  dimensions  of  a  moment,  and  repre¬ 
sents  the  moment  applied  to  the  air  about  the  centroid  of  vorticity  by 
one  side  of  the  wing.  It  can  quite  easily  be  evaluated  mechanically  if 
the  spanwise  loadmg  is  known  by  taking  the  moment  of  this  load  about 
the  cenfo'd  of  vorticity. 


We  now  notj  thnt  we  have  developed  expressions  for  the  six  major 
parameters  of  the  trailing  vortex  sheet.  These  are  the  lift,  the  induced 
drag,  the  moment  of  momentum,  the  vortex  span,  the  total  bound  vorti¬ 
city,  and  the  profile  drag.  The  first  five  of  these  can  be  expressed  di¬ 
rectly  from  an  analysis  of  the  spanwise  loading  curve.  Only  four  of 
these  are  actually  independent  in  thet  can  be  determined  from  by 
and  L  .  The  profile  drag  can  readily  be  determined  from  standard 
drag  determination  methods.  The  total  drag  D^,  =  Dq  +  . 

Thus  we  can  consider  the  vorticity  sheet  shed  at  the  wing  trailing 
edge  to  be  determined  (in  a  global  sense)  by  L,  Drp,  M  ,  by,  P0  ■ 
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Rollup  and  Development  of  the  Vortex  Sheet 


After  leaving  the  v/ing,  the  shed  vorticlty  undergoes  a 
series  of  complicated  convolutions  until  it  is  finally  diffused  into  the 
atmosphere.  We  can  arbitrarily  define  four  stages  to  this  process, 
in  which  fairly  distinct  interactions  occur.  These  are  shown  in  Fig¬ 
ure  2.1  and  can  be  thought  of  as  the  Planar  Sheet  Stage,  the  Rollup 
Stage,  the  Viscous  Vortex  pair  stage,  and  the  Distributed  Impulse 
stage.  Our  main  interest  i8  in  the  latter  two.  They  are,  of  course, 
controlled  by  the  initial  conditions,  which  come  from  the  Planar  Stage 
discussed  in  the  previous  section. 

After  leaving  the  wing,  various  mutual  inductions  cause  the 
vortex  sheet  to  curl  up  at  the  tips  and  to  form  a  scroll-like  shape  and 
initiate  the  rollup  stage.  Within  about  20  spans  this  has  developed  into 
a  viscous  vortex  pair  of  strength  and  span  ro,  by  and  having  a  finite 
core  within  which  the  vorticity  is  contained.  Many  authors  (see  Bib¬ 
liography)  have  attempted  to  describe  this  rollup  process  but  the  pro¬ 
per  analysis  is  still  highly  controversial.  We  will  not  specifically  dis¬ 
cuss  these  analyses  here  except  to  state  that  the  results  of  numerical 
free  vortex  schemes  should  be  regarded  with  caution.  While  these 
techniques  in  principal  exactly  treat  the  self-induced  development  of 
the  trailing  vortex  sheet,  the  necessary  discretization  for  computer 
operation  may  be  expected  to  introduce  instabilities  and  motions  re¬ 
lated  more  to  the  discretization  and  computer  scheme  than  to  the  ac¬ 
tual  deformation  of  the  vortex  sheet.  This  is  particularly  true  if  the 
vortex  sheet  is  modeled  by  a  set  of  discrete  line  vortices.  For  our 
purposes  we  are  not  concerned  with  the  details  of  the  rollup  process 
because  this  occurs  so  close  to  the  wing. 

The  next  stage  is  the  viscous  vortex  development,  which  is 
discussed  in  more  detail  in  a  following  chapter,  and  finally  the  decay  to 


Stages  of  Wake  Vortex  Development 


the  distributed  impulse  (also  considered  in  more  detail  in  the  Trans¬ 
port  Section,  4 ) . 

To  assist  in  the  discussion  of  this  development,  we  de¬ 
scribe  the  various  global  invariants  in  the  process.  Evidently  the 
lift  and  drag  must  be  conserved  throughout  the  development.  Now 
we  note  the  other  terms;  shed  vorticity  and  moment  of  momentum 
are  actually  null  for  the  entire  wing  and  can  be  defined  only  for  one 
half  wing.  However,  provided  there  is  no  viscous  interaction  be¬ 
tween  the  two  sides  (left  and  right)  of  the  vortex  system,  then  even 
in  a  viscous  fluid  we  will  find  that  and  rQ  are  conserved.  We 
note  by  the  kinematic  condition  of  the  solenoidal  quality  of  vorticity, 
that  even  though  the  vorticity  is  diffused  into  originally  irrotational 
flow,  ro  must  be  conserved  until  vorticity  from  the  one  core  inter¬ 
acts  with  the  other  on  the  center  line.  The  moment  of  momentum  is 
conserved  according  to  the  same  arguments.  A  physical  interpre¬ 
tation  of  this  is  that  no  external  torque  can  be  applied  to  the  vortex 
core  and  its  surrounding  flow  system  until  appreciable  vorticity  has 
reached  the  center  line.  We  comment  here  that  as  long  as  the  flow 
on  the  center  line  is  symmetrical  even  introducing  viscous  transfer 
in  the  equations  will  produce  no  shear  along  the  center  line  since 
velocity  gradients  there  are  zero. 


Thus  during  the  rollup  and  the  viscous  pair  stage  a  good  as¬ 
sumption  may  be  to  conserve  and  Fo  .  In  the  final  decaying  im¬ 
pulse  stage  we  will  expect  and  to  decay  according  to  transfer 
processes,  both  of  molecular  viscosity  and  turbulence. 

We  have  not  discussed  the  variation  of  H  =  p/p  +  i-  Va  the 
total  head  of  the  flow.  This  provides  a  local  flow  property  (as  opposed 
to  global  properties)  which  is  useful  to  us.  We  note  that  at  the  wing  trail¬ 
ing  edge  there  is  already  a  sheet  of  low  energy  air  shed  from  the  boundary 
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layer,  having  dimensions  scaled  by  the  span  and  the  trailing  edge  bound¬ 
ary  layer  displacement  thickness, and  containing  a  total  head  deficit  re¬ 
lating  to  the  boundary  layer  momentum  thickness.  This  wake  will  en¬ 
train  external  flow  in  the  way  normally  associated  with  a  non-lifting 
velocity  deficit  wake. 

Of  more  interest  is  the  entrainment  and  head  loss  due  to 
the  rotational  flow  associated  with  the  sidewash  and  downwash  flows 
connected  to  the  induced  drag.  This  effect  features  in  the  swirling  vis¬ 
cous  interaction  which  develops  the  axial  flow  in  the  vortex  core.  For 
convenience  we  define  H.  as  the  head  loss  connected  with  the  lifting 

J_> 

flow.  It  seems  that  this  term  becomes  significant  only  in  the  vortex 
pair  and  decay  stages. 

It  is  noted  that  the  global  force  and  local  energy  relation¬ 
ships  are  not  sufficient  to  solve  for  the  details  of  the  vortex  develop¬ 
ment.  However,  they  provide  a  useful  insight  into  the  processes  oc- 
curing  and  can  be  very  helpful  in  assessing  different  theoretical  analy¬ 
ses,  most  of  which  contain  direct  or  implicit  statements  relating  to 
assumptions  concerning  these  quantities. 


3.  THE  ENCOUNTER  HAZARD 


Evidently,  if  conditions  are  such  that  an  aircraft  encounters  a 
vortex  wake  while  the  wake  is  shill  in  the  organized  state,  consisting 
of  two  near-axial  viscous  vorticuB,  the  hazard  will  be  a  function  of 
some  parameters  describing  the  strength  and  geometry  of  the  wake, 
coupled  with  further  factors  relating  to  the  control  power  and  design 
load  factors  of  the  encounter  aircraft.  Other  parameters  in  this  situ¬ 
ation  are  the  flight  path  relative  to  the  wake  of  the  generating  aircraft, 
and,  of  course,  the  proximity  to  the  ground  of  the  encounter. 


Initially,  we  shall  consider  the  case  of  the  encounter  aircraft 
entering  the  eye  of  one  of  the  trailing  vortices  on  a  flight  path  approxi¬ 
mately  parallel  to  the  vortex  axis.  In  these  c  ircumstances  a  severe 
rolling  moment  will  be  induced  on  the  encounter  aircraft.  It  would  at 
first  appear  that  the  magnitude  of  the  rolling  moment  would  be  closely 
related  to  the  core  development,  that  is,  the  vortex  age.  However, 
when  one  assumes  that  the  core  will  be  less  than  2/10  of  the  gener¬ 
ating  span,  and  that  outside  this  the  flow  very  closely  approaches  that 
of  a  horseshoe  vortex,  we  see  that  much  of  the  encounter  aircraft  wing 
for  an  encounter  span  of  about  1/2  the  vortex  span,will  be  in  the  invis- 
cid  portion  of  the  flow.  Coupling  this  with  the  fact  that  the  contribu¬ 
tions  to  roiling  moment  vary  directly  with  the  spanwise  position,  we 
may  expect  that  the  actual  size  of  the  core  will  not  be  too  important 
providing  it  is  less  than  a  certain  fraction  of  the  encounter  span.  On 
this  basis,  Crow  ( 1 97 1 )  has  shown  that  the  induced  rolling  moment  will 


be  proportional  to  f, /b  a  ,  where  F  is  the  strength  of  the  generator 
aircraft  vortex  and  the  span  of  the  encounter  aircraft.  If  it  is  then 
assumed  that  there  is  a  certain  maximum  aileron  induced  roll  rate  of 


the  encounter  aircraft,  P  ,  one  can  assert  that  hazardous  condi- 

rnax 

tiong  will  be  obtained  when  the  ratio  of  vortex  induced  roll  to  aileron 


rolling  power  exceeds  a  certain  value.  Writing  the  trailing  vortex 
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strength  in  terms  of  lift  coefficient,  ,  span  b  ,  and  aspect  ratio,  A 
we  define  this  ratio  below. 

If  the  critical  rolling  rate  of  the  encounter  aircraft  is  P  (in 
radians  per  second),  the  non-dimensional  danger  factor  D  ,  expressing 
the  hazard,  can  be  written  as 

D  = 

where  the  suffice  g  applies  to  the  generator  aircraft  and  e  to  the  en¬ 
counter  aircraft.  It  has  been  estimated  by  Crow  from  Boeing  tests  by 
Condit  and  Tracy  (1971)  that  a  value  for  D  of  about  unity  is  a  good 
measure  of  critical  hazard.  If  D  exceeds  this  value,  dangerous  bank 
attitudes  may  be  expected.  At  low  altitudes  a  smaller  D  value  will  be 
appropriate  since  recovery  time  is  reduced. 

The  danger  factor  has  a  very  compact  formulation  and  it  will  be 
noted  that  it  contains  no  term  relating  to  the  distance  behind  the  gener¬ 
ator  aircraft.  This  is  a  significant  simplification  and  relates  to  the 
slow  growth  of  the  core  and  the  invariance  of  the  total  circulation.  A 
more  detailed  but  still  idealized  analysis  of  the  induced  rolling  moment 
conducted  by  Condit  and  Tracy  (1971)  shows  the  separation  distance  de¬ 
pendence.  Taking  the  results  of  their  numerical  calculations,  we  see 
that  the  danger  factor  is  in  fact  almost  invariant  with  separation  dis¬ 
tance.  For  example,  using  dimensions  appropriate  to  a  B707  we  find 
that  the  more  precise  analytical  calculations  will  give  ?.  normalized  re¬ 
duction  in  induced  rolling  moment  of  about  2.5*%  per  nautical  mile  be¬ 
hind  generator  aircraft.  These  numbers  are  substantiated  by  fligM 
tests  on  Boeing  aircraft  reported  in  the  sar.  paper. 

It  is  noted  that  at  first  this  appears  to  contradict  the  results  of 
Bisgood,  Maltby  and  Dee  (1971)  who  use  Squire’s  Theory  to  compute 
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wake  dissipation  and  exhibit  a  theoretical  curve  showing  a  rather  marked 
reduction  in  rolling  rate  with  vortex  age.  This  curve  is  for  an  assumed 
eddy  viscosity  of  -  0041*  {1^  being  the  vortex  strength).  However,  the 

paper  chows  eddy  viscosity  inferred  from  flight  measurement  of  core 
velocity  peaks,  indicating  that  the  effective  eddy  viscosity  according  to 
this  theory  is  more  like  .  0002ro  with  some  readings  as  low  as. 000181^. 
Using  the  value  of  .  0002  f^and  conducting  a  similar  analysis  would  show 
the  reduction  in  rolling  rate  to  be  about  5%  per  nautical  mile,  which  is 
more  in  line  with  the  more  modern  and  comprehensive  Boeing  tests. 

Now,  we  should  consider  the  other  real  factors  associated  with 
vortex  entry  in  a  near  axial  direction.  These  relate  to  control  re¬ 
versals  as  the  encounter  aircraft  traverses  the  vortex  laterally  or 
enters  the  other  vortex, as  well  as  the  coupled  yaw-roll  dynamics  as¬ 
sociated  with  the  effect  of  the  vertical  stabilizer.  Such  excursions  due 
to  random  oblique  encounters  will  introduce  hazard  elements  of  more 
significance  than  those  computable  by  simple  axial  vortex  decay  and 
thus  justify  ignoring  the  latter  effect. 

On  these  grounds  it  is  believed  that  the  danger  factor  Li  repre¬ 
sents  a  useful  criterion  for  the  vortex  hazard  due  to  near  axial  en¬ 
counters,  and  the  assumption  that  D  is  invariant  with  separation  dis¬ 
tance  ia  a  significant  and  valid  simplification.  The  critical  value  of  D 
will  depend  on  the  aircraft  ground  proximity. 

It  must  be  remarked  here  that  w-hilo  the  danger  factor  is  invari¬ 
ant  with  vortex  age  or  separation  distance,  its  existence  depends  upon 
there  being  an  organised  viscous  vortex;  thus  once  the  regular  motion 
has  been  disrupted  by  sinuous  or  core  instability  the  danger  factor  will 
change.  This  illustrates  the  importance  of  the  understanding  and 
prediction  of  the  catastrophic  instabilities. 

The  other  critical  mode,  the  encounter  normal  to  the  trailing  vor¬ 
tices,  is  somewhat  less  simple  to  handle.  In  pi\icsple  it  involves  the 


same  /actors.  The  disturbance  syBtem  will  be  described  in  terms  of  the 
peak  vertical  velocities  and  their  spacing,  while  the  response  system 
now  relates  to  the  longitudinal  transfer  function  of  the  encounter  air¬ 
craft.  Confining  the  analysis  to  the  encounter  with  only  one  core  of  the 
vortex  pair.it  has  been  shown  by  Houbolt  (1971)  that  the  maximum  norm¬ 
al  acceleration  N  of  the  encounter  aircraft  can  be  written  in  a  form 
involving  the  wing  loadings  and  speeds  of  both  aircraft,  the  aspect  ratio 
of  the  generating  aircraft  and  the-  transfer  function  of  the  encounter  air¬ 
craft.  This  transfer  function  (k^)  is  of  course  related  to  relative  mass 
of  the  encounter  aircraft  as  well  as  the  ratio  of  its  chord  to  the  vortex 
core.  A  curve  for  this  is  given  by  Houbolt,  which  shows  fairly  weak 
variation  with  relative  mass  and  core  radius.  »Ve  write  Houbolt's 
equation  in  the  same  terms  as  the  Crow  axial  danger  factor,  and  as¬ 
sume  a  relatively  high  value  of  k^  of  .  85  (to  round  off  numbers),  then 
obtain  a  transverse  danger  factor  D^.  as 

O-  =  20(UC.  /ttA)  ,"VUCt  ) 
t  g  JL>  e 

where  N  is  the  limit  normal  acceleration  factor  of  the  encounter  air¬ 
craft,  On  this  basts,  a  value  of  D^.  exceeding  unity  may  be  regarded 
as  hazardous. 

We  now  note  that  the  ratio  of  these  terms  will  indicate  the  relative 
severity  of  axial  or  normal  encounter.  Thus  we  obtain 

D/Dr  =  b  (NUC.  /20  b* 

T  g  C  max  o 

bP 

Noting  that  £  normalised  roll  rate.  — s  „  of  about  .  07  i#  usually  at¬ 
tained.  we  car.  write  ths*  as 
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D/DX  -  0-35(NCL)e  bg/b 


Now,  taking  representative  transport  aircraft  values  of  N  -  4.  5  , 

~  .  6  we  see  that  the  ratio  of  the  danger  factors  is  approximately 
equal  to  the  span  ratio,  showing  that  in  all  significant  case«  (smaller 
encounter  air  ' raft)  D>  so  that  the  axial  encounter  can  be  expected 
to  be  more  hazardous. 
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4.  VORTEX  TRANSPORT 


4  .  1  GENERAL 

The  trailing  vortex  pair  creates  an  induced  field  which 
causes  the  vortex  system  to  convect  downwards.  This  vertical  speed 
is  relatively  small  compared  to  the  flight  speed.  Making  the  classical 
assumptions  of  elliptical  loading  and  inviscid  rollup  to  a  finite  vortex 
pair,  we  obtain  for  the  downwash  velocity  U* 

U*/U  =  4C,  /r^A 


We  note  for  a  typical  wing  of  aspect  ratio  6  that' the  downward  motion 
may  be  about  1%  of  the  flight  cruise  speed,  increasing  to  3  or  4%  at  the 
landing  state. 

This  can  be  compared  with  the  same  ratio  before  the  sheet 
has  rolled  up,  which  is  about  {C^/rr)[l/A  +  1 /n}  at  the  trailing  edge, 
taking  into  account  the  bound  vortex  term,  and  about  2C^/rrA  a  few 
spans  downstream.  Thus  for  the  first  few  spans  the  downwash  may  be 
about  four  times  that  of  the  final  state;  however,  our  interest  is  with 
the  motion  after  the  vortex  pair  has  developed.  Field  observations 
have  indicated  that  apparently  in  all  cases  the  vortex  dercsnt  rate  is 
approximately  constant  for  30  -  40  seconds,  then  begins  to  reduce. 
Condit  and  Tracy  (1971)  report  that  in  no  cases  have  descents  of  more 
than  about  1000  ft  been  observed,  while  the  recent  experiments  of 
Tombach  (1972)  show  a  distinct  reduction  in  transport  velocity.  How¬ 
ever,  we  stress  that  definitive  measurements  of  vortex  descent  rate 
are  lacking.  We  discuss  this  further  in  the  vortex  decay  section.  It 
appears  clear  that  after  about  60  seconds  the  descent  rate  is  not  ade¬ 
quately  described  by  the  classical  equations. 
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If  the  pair  is  immersed  in  a  uniform  crosswind,  standard 
calculations  indicate  that  the  pair  should  simply  be  laterally  convected 
with  the  crosswind,  with  no  change  in  cell  shape  or  vortex  dynamics. 

This  appears  to  be  the  case. 

For  a  vertically  non-uniform  crosswind  the  situation  is  less 
simple.  The  pair  is  now  immersed  in  a  field  containing  a  directed  or¬ 
ganized  vorticity  and  one  would  expect  various  unsymmetrical  effects 
such  as  vortex  tilting  to  occur. 

The  simplest  case,  that  of  a  uniformly  sheared  crosswind 
(linear  wind  profile),  can  be  solved  exactly  in  the  inviscid  case. 

This  pass  -.»es  a  steady  solution  in  a  frame  of  reference  fixed  in  the 
vortex  cores,  and  for  the  inviscid  model  gives  no  banking  or  change  in 
trans  .ort  speed.  This  analysis  is  developed  in  a  later  section.  It  is 
simply  noted  here  that  the  reason  this  problem  can  be  solved  is  that 
although  the  external  vorticity  is  redistributed  by  the  vortex  field,  this 
does  not  introduce  my  unsteadiness  in  the  flow  since  the  redistribution 
does  not  change  the  uniform  vorticity  field  in  which  the  cell  is  immersed. 

If  one  considers  a  field  having  a  shear  gradient,  then  even 
the  inviscid  problem  becomes  unsteady  because  of  the  redistribution  of 
the  non-uniform  ambient  vorticity.  An  unsteady  model  was  utilized  by 
Burnham  (1972)  ir.  attempting  to  model  the  ground  effect  case  with  a 
non-uniform  crosswind.  A  finite  set  of  discrete  vortices  was  introduced 
into  the  ambient  field  and  permitted  to  convect  with  the  flow.  This  model 
is  an  aUempt  to  explain  the  observed  rising  of  one  of  the  vortices  near 
the  ground.  It  is  believed  that  this  model  would  give  different  quantita¬ 
tive  vortex  paths  for  different  or  more  detailed  models  of  the  same  cross 
wind  field.  We  note  too,  that  Bisgood,  Maltby,  and  Dee  (1971)  modeled 
the  vortex  trajectory  in  ground  effect  with  a  non-uniform  crosswind 
simply  by  computing  the  convection  at  each  height  and  taking  no  account 


of  voracity  redistribution.  This  simple  kinematic  model  gave  excell¬ 
ent  correlation  with  the  observed  vortex  paths,  at  least  for  the  few 
reported  tests. 

We  note  here  that  a  large  proportion  of  the  tests  con¬ 
ducted  by  TornbacH  exhibited  unsymmetrical  vortex  response  under 
crosswinds.  This  was  normally  observed  in  the  vortex  pair  banking, 
and  occasionally  in  the  collapse  of  one  vortex  due  to  core  bursting. 

Such  behavior  might  be  expected  under  the  ambient  conditions  described, 
where  interaction  between  the  ambient  and  core  vorticity  could  cause 
strong  anti-symmetric  effects  both  with  respect  to  vor'ax  transport  and 
core  development.  Changes  in  core  development  could  accelerate  the 
instability  of  one  element  of  the  pair. 

One  of  the  more  important  effects  of  uiir.ymmetrical  core 
bursting  might  be  expected  to  be  the  development  of  the  solitary  vortex. 
In  a  number  of  tests  reported  by  Tombach  it  was  observed  that  when 
only  one  vortex  remained  (or  at  least  was  visible)  this  vortex  appeared 
to  ha/e  a  much  longer  life  than  usual,  and  did  not  show  any  signs  of 
sinuous  Crow  Instability. 

It  appears  that  if  these  unsymmetrical  effects  are  due  to 
wind  shear,  then  they  should  certainly  be  primarily  attributable  to  the 
uniform  shear  case  and  that  it  should  not  be  necessary  to  postulate  a 
non-uniform  shear  to  describe  them.  The  basis  for  this  is  twofold, 
the  first  reason  being  that  these  effects  have  been  observed  at  altitu¬ 
des  at  which  the  uniform  shear  is  weak  and  the  shear  gradient  even 
weaker;  the  second  that  any  viscous  model  will  indeed  give  anti¬ 
symmetric  behavior  when  immersed  in  a  uniform  shear  field. 

Thus  while  an  inviscid  uniform  shear  model  will  not  give 
any  vortex  tilting,  it  is  believed  that  the  uniform  shear  is  the  dominant 
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effect  and  that  a  combined  uniform  shear/viscous  decay  model  should 
be  sufficient  to  account  for  vortex  behavior  in  a  general  crosswind. 
This  is  described  in  Section  4 . 3 


A  further  factor  complicating  the  descent  dynamics  is  the 
effect  of  flow  stratification.  The  previous  paragraphs  have  described 
the  situation  in  a  homogeneous  fluid  with  the  main  effects  due  to  en¬ 
trainment  or  ambient  organized  vorticity.  If  the  atmosphere  is  strati¬ 
fied  as  is  frequently  the  case,  one  would  certainly  expect  this  to  have 
some  effect  on  the  descent  rate.  At  the  moment,  both  the  sense  of  the 
effect  and  its  magnitude  are  controversial.  Contrary  to  intuition,  it 
appears  that  descen:  does  not  necessarily  slow  down  in  a  stable  strati¬ 
fication;  that  is,  the  cell  is  not  necessarily  retarded  on  penetrating  the 
denser  medium.  Physically,  this  anomalous  behavior  is  accounted 
for  by  the  contraction  of  the  vortex  span  --  the  problem  is  whether 
the  span  does  in  fact  expand  or  contract  and  at  what  rate.  In  a  later 
section  we  develop  a  new  theory  which  shows  that  an  inviscid,  strati¬ 
fied  model  does  i~  ’.y  a  span  contraction  and  consequent  acceleration 
of  the  downward  descent.  More  details  and  an  evaluation  of  these  as¬ 
sumptions  are  made  in  the  section. 

li  should  always  be  recalled  that  in  practically  all  cases 
the  organized  vortex  pair  motion  is  terminated  by  sinuous  instability 
or  by  core  bursting.  This  is  the  element  which  makes  experimental 
data  on  descent  rate  so  incomplete. 

In  the  following  sections  we  assume  that  an  organized  vor¬ 
tex  pair  motion  exists  and  discuss  in  more  detail  three  important  cases: 
the  descent  in  a  homogeneous  Lrrotational  fluid,  the  descent  in  a  homo¬ 
geneous  fluid  of  uniform  rotation,  and  the  descent  in  an  irrotational 
fluid  of  uniform  stratification. 
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Thus,  if  we  consider  the  descent  rate  to  be  characterized 
by  U*  =  U*(Z,  l/T)  where  Z  is  the  ambient  vorticity  and  T  the  Brunt- 
Vaisala  time  (a  measure  of  the  stratification)  we  attempt  to  describe 
the  problems 


U*  =  U*(0,0) 

U*  =  U*(Z,0) 

U*  =  U*(0,1/T) 

Evidently  a  full  analysis  of  the  general  problem  cannot  be  achieved  un¬ 
til  the  special  cases  above  are  understood  and  quantified.  However, 
we  note  that  a  validation  of  the  above  cases  may  be  difficult  because  it 
must  depend  on  one  of  the  atmospheric  parameters  approaching  zero 
in  the  field  test. 
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4.2  NEUTRAL  IRROTATIONAL  FIELD 


After  the  wing  vortex  sheet  has  rolled  up, the  trailing  s/s¬ 
tem  consists  of  a  pair  of  vortices  of  finite  rotational  core  area,  but 
with  the  core  radius  a  relatively  small  fraction  of  the  vortex  span.  If 
this  vortex  pair  is  immersed  in  a  still  homogeneous  inviscid  flow,  a 
well  known  classical  solution  exists  under  which  the  pair  is  convected 
downwards  at  a  velocity  p  /2 nb.  The  classical  analysis  shows  that 
there  is  a  closed  recirculating  mass  of  air,  of  roughly  oval  propor¬ 
tions,  associated  with  the  concentrated  vortex  pair,  and  that  this  cell 
is  convected  downwards  at  a  uniform  speed.  Flow  exterior  to  this 
cell  never  enters  it.  On  this  basis  a  long  vortex  pair,  which  may  be 
regarded  as  substantially  two-dimensional,  will  move  downwards  in 
an  unbounded  fluid  with  constant  velocity  for  all  times. 

In  real  flows  this  situation  does  not  persist  indefinitely, 
and  most  experiments  appear  to  show  the  rate  of  descent  to  reduce 
and  finally  approach  zero.  This  is  certainly  due  to  diffusion  of  core 
vorticity  by  some  combination  of  laminar  and  turbulent  viscosity  and 
will  occur  even  in  homogeneous  (unstratified)  flows.  Much  effort  has 
gone  into  explaining  and  quantifying  this  effect,  but  the  subject  still  re¬ 
mains  controversial.  Recently  a  new  rational  interpretation  of  this, 
coupled  with  careful  observations,  has  been  put  forward  by  Maxworthy 
(1972)  which  greatly  assists  explanation  of  the  effect. 

Maxworthy  conduct’d  experiments  with  vortex  rings  in 
water,  using  various  visualization  techniques  to  identify  where  the 
flow  went.  When  the  vorticity  was  relatively  well  distributed  in  the 
ring  he  observed  that  the  outer  flow  was  entrained  into  the  back  of  the 
cell,  causing  an  increase  in  the  cell  volume.  At  the  same  time  a  por¬ 
tion  of  the  cell  vorticity  was  shed  into  the  wake,  removing  both  vorti¬ 
city  and  momentum  from  the  cell.  The  combined  effect  of  this  is  to 
increase  the  cell  size  and  to  reduce  its  propagation  velocity. 
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The  mechanism  of  entrainment  is  important  for  our  further 
development.  Figure  4-1  shows  a  sketch  of  this  flow  field,  in  coordin¬ 
ates  fixed  at  the  core  centers,  so  that  the  outer  flow  is  represented  by 
a  uniform  but  unsteady  flow  from  below.  We  note  that  the  cell  has  a 
well  defined  stagnation  point,  A  ,  and  over  the  front  portion,  a  well  de¬ 
fined  cell  boundary,  A-B  .  Across  this  boundary  the  pressure  and 
velocity  fields  of  inner  and  outer  flows  are  continuous,  the  only  discon¬ 
tinuity  being  between  the  inner  vortical  fluid  and  the  outer  irrotational 
flow.  Due  to  both  laminar  and  turbulent  effects,  this  vorticity  is  trans¬ 
ferred  to  the  outer  flow,  and  as  a  consequence  the  total  head  of  the  flow 
is  reduced.  Thus,  after  passing  the  maximum  velocity  point  near  B  , 
the  outer  flow,  contained  approximately  by  the  stream  tube  CD  ,  is  un¬ 
able  to  recover  sufficient  velocity  to  rejoin  the  outer  flow  at  the  rear, 
but  remains  as  part  of  the  stationary  cell.  Thus  the  cell  size  is  in¬ 
creased.  At  the  same  time,  a  neighboring  stream  tube  EF  acquires 
a  smaller  amount  of  vorticity  and  suffers  less  head  reduction,  so  that 
is  does  depart  from  the  cell  at  the  rear,  but  at  a  lower  than  free  stream 
total  head.  This  portion  develops  into  a  wake  behind  the  cell. 

Thus  the  same  process  causes  entrainment  of  the  outer 
flow  into  the  cell  and  a  detrainment  (removal)  of  some  of  the  cell  vor¬ 
ticity  and  momentum.  A  further  process  occurs  on  the  centerline  of 
the  cell,  AX.  Here  vorticity  is  annihilated  by  diffusion  from  the  left 
and  right  cells.  Thus  three  vorticity  transfer  mechanisms  occur  and 
the  overall  effect  controls  the  cell  dynamics. 

Maxworthy  showed  that  initially  the  vortex  shedding  to  form 
the  wake  was  extremely  weak,  since  the  cell  vorticity  at  the  boundary 
was  quite  weak.  Thus,  although  the  cell  grew  in  size,  it  did  not  lose 
momentum,  and  the  impulse  I  was  conserved.  In  these  circumstances 
the  main  vorticity  loss  occurred  along  the  centerline  and  was  small,  and 
there  was  minimal  wake  momentum  loss.  During  the  later  stages  in 
growth  history,  when  more  vorticity  is  present  near  the  boundary,  the 
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wake  develops;  thus  the  impulse  in  the  cell  reduces  while  the  cell  size 
increases.  Both  these  effects  contribute  to  the  reduction  in  speed  and 
final  complete  annihilation  of  the  cell  momentum. 

It  must  be  noted  that  Maxworthy's  experiments  were  con¬ 
ducted  with  vortex  rings  at  extremely  low  Reynolds  Numbers,  when  the 
flow  was  certainly  laminar.  Maxworthy  has  indicated  in  a  private 
■  munication  that  further  flow  visualization  tests  with  finite  wings 
aioo  exhibited  a  detrained  wake.  These  experiments  were  also  at  very 
low  Reynolds  Numbers.  It  is,  however,  possible  that  during  the  later 
stages  of  development  of  an  aircraft  trailing  vortex  system  that  similar 
processes  of  mass  entrainment  and  momentum  detrainment  occur.  For 
laminar  transfer  the  time  scales  would  be  too  long  to  be  of  interest,  but 
if  the  transfer  is  assumed  turbulent,  it  may  be  possible  to  account  for 
some  of  the  observed  effects.  Thus  it  appears  very  probable  that  the 
later  development  of  a  vortex  pair  follows  qualitatively  the  stages  de¬ 
scribed  by  Maxworthy;  with  an  additional  initial  stage  which  we  will 
discuss.  We  therefore  postulate  three  stages  as  shown  in  Figure  4-2. 

Stage  I  -  The  Inviscid  Cell.  Here  the  vorticity  ;•  confined 
to  well  within  the  cell  boundary.  On  the  boundary  itself  there  will  be  no 
laminar  transfer  (since  there  is  no  distortion)  and  turbulent  transfer  will 
have  no  net  effect,  since  both  inner  and  oviter  flow  have  the  same  total 
head.  In  these  circumstances  the  inviscid  cell  model  will  be  a  good  re¬ 
presentation  of  the  dynamics  and  we  find  that  the  time  rates  of  change  of 

•  9 

cell  size,  b  ,  and  impulse  I  are  zero;  so  that  the  propagation  velocity 

_  .*  .  v 

U  ts  constant. 

Stage  II  -  The  Entraining  Cell.  As  the  core  vorticity  dif¬ 
fuses  and  approaches  the  cell  boundary,  the  first  process  (of  mass  en- 

•  • 

trainmcnt)  occurs,  represented  by  bv  >  0  ,  1  =  0  .  We  find  that  the  pro¬ 
pagation  velocity  reduces  from  the  inviscid  value. 
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Stage  III  -  The  Decaying  Cell.  During  the  later  8tage6, 

substantial  momentum  and  vorticity  shedding  occur,  causing  a  wake  to 

*  * 

develop  and  giving  by>  0,  I<  0.  Of  course,  various  catastrophic  insta¬ 
bilities  may  have  developed  before  the  complete  decay  has  occurred. 


1 .  Ordering  of  the  Process 

Maxworthy  has  made  a  laminar  analysis  of  the  last  two 
stages  for  a  ring  vortex,  and  found  excellent  correlation  with  his  ex¬ 
periments.  He  also  includes  an  analysis  for  the  two  dimensional  vor¬ 
tex  pair  for  laminar  diffusion.  Here  we  continue  this  for  turbulent  dif¬ 
fusion. 


To  establish  a  turbulent  transfer  coefficient,  we  note  that 
there  will  presumably  be  two  turbulent  scales  at  the  cell  boundary. 

During  the  initial  stage  we  assume  that  ambient  turbulence  is  the  major 
term,  since  the  inner  flow  near  the  cell  boundary  has  very  weak  vorticity. 

If  this  is  the  case,  the  only  turbulent  parameter  is  „  ,  the  turbulent 

dtsstoation.  and  we  c“in  define  an  rms  turbulent  velocity  V  for  all 
*  s 

scales  smaller  than  the  coll  scale  bv  by 

v  ~  e’*/3b  1/3 

€  v 

In  a  later  stage,  we  expect  mechanically -gene  rated  turbu¬ 
lence  due  to  the  interior  cel!  motion  to  dominate.  No'v  u*e  w’iU  tine  a 
turbulent  velocity  coupled  to  the  celt  velocity,  eo  define  another 
turbulent  scale  velocity 


V 

m 


U* 
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a.  The  Entraining  Cell.  Here  we  assume  vorticity  is  en¬ 
trained  linearly  with  time  as  the  layer  progresses  from  A  to  B  . 
Thus  the  height  is  given  by 


6 


e 


V  t' 

e 


where  t1  is  the  time  of  travel.  An  estimate  of  f  is 

t'  =  TTbv/(2U+) 


for  the  inviscid  cell  proportions.  Now  the  unit  density  mass  flow 
through  BD  is  .  and  taking  the  cell  area  to  be  about  r*  h ^  and 

assuming  all  this  mass  is  entrained,  we  obtain 

b  ~  V  h  ~ 

v  c  '  v  e  1 

Taking  the  total  impulse  to  be  that  of  the  inviscid  cell,  I  -  Tb 
and  the  propagation  speed  U*  =  r/2rrb  .  Thus  for  I  invariant  we  obtain 

U*  =  I/(2net3) 

W'e  compare  this  with  Maxworlhy’s  result  for  laminar  diffusion  of  kine- 
matic  viscosity.  .  .  where  he  gives 


U*  -  (l 


>  2,1/5  -4/5 

<  l  ( 


The  time  dependence  and  the  viscous  terms  appear  in  a  slmstar 
fashion  in  both  the  turbulent  and  laminar  cases,  however  is  the  turbulent 
cate  the  departure  from  the  :nviecid  celt  results,  as  represented  by  the 
exponents,  are  much  more  distinct. 
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b.  The  Decaying  Cell.  In  thin  case  we  admit  both  entrai  tent 
of  mass  and  detrainment  of  momentum;  thus  the  impulse  is  no  Ion,  con¬ 
served.  We  assume  a  *urbulent  velocity  scale  V^,  which  will  be  defined 
later.  Then,  as  before,  we  obtain  the  cell  growth  equation 


where  a  is  some  constant. 

Now,  to  find  the  ehange  in  impulse,  we  estimate  the  reduc¬ 
tion  in  circulation.  We  ignore  vorticity  annihilation  on  the  interior  cell 
centerline  but  attempt  to  compute  the  vorticity  transferred  to  the  wake 

by  turbulence  through  the  cell  bounda  ry.  Dimensional  arguments  sug- 

2 

gest  that  the  cell  vorticity  goes  like  r/b  and  the  turbulent  flux  of  vor- 

v  2 

ticity  through  the  boundary  then  becomes  of  the  order  (r/b.  1  VT  . 
Thus  wc  obtain 


r  =  a  r  b  v 


V  'T 


where  jj  is  another  constant.  Combining  these  equations  with  §/•>  -  y 
gives  the  result 

hr  --  ybv/bv 

independent  of  .  implying  that  r  -  bv  v  . 

if  we  now  assume  V_.  -  Us  we  obtain  she  results 


o  -  t 


'  - 


■>r> 


,-Sr,? 


I  4s?:; 


.  ~i 

m 


^v-V 


■SiSS?:*'--?  ; 
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F  •=“;  !>’ 
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We  observe  that  for  the  impulse  to  decay,  v  >  1  ;  however,  the  reduc¬ 
tion  of  propagation  velocity  will  be  leoa  than  that  in  the  case  of  assum¬ 
ing  turbulence  coupled  to  the  ambient  dissipation.  We  note  also  an  in¬ 
herent  contradiction  in  that  if  the  mechanical  turbulence  is  coupled  to  the 


U  ,  then  eventually  it  will  reduce  so  that  ambient  turbulence  again  dom¬ 


inates. 


<1.  Comparison  of  the  Assumptions 


We  note  that  in  this  simple  model  we  can  assume  that 
the  impulse  is  conserved  or  not,  and  that  the  turbulent  scale  is  either 
the  ambient  turbulent  velocity  or  the  propagation  velocity.  Thus  we 
actually  obtain  four  cases  for  similarity  solutions.  The  ordering  of 
the  response  is  shown  in  the  following  table. 


TABU.  4  1.  TUT  UU'KRINu  OF  THl.  RhSFONSK 
’•'OR  FOUn  ASFS  OF  SIMILARITY 


Turbulent 

Mechanism 


Ambient 


Mechanical 


No 

Wake 


W  a’  o 


,~3{y<U/2 


t-(Y-l}/tY*2 


We  note  that  even  on  this  simple  basis,  four  different 

similarity  solutions  can  be  developed.  It  is  possible  that  the  stages  with 

time  progress  from  the  ambient/no  wake  to  mechanical  turbulence /no 

wake  case  through  the  mechanical  turbulence/wake  case  to  the  ambient/ 

turbulence/wake  case.  However,  this  assumption,  although  sounding 

physically  plausible,  gives  a  non-monotonic  variation  in  the  propagation 

speed  decay  rate,  which  appears  surprising.  For  example,  if  y=  1  +  ft  , 

we  would  find  the  propagation  speed  varied  like  t'3,  t"2^3  t"(2+6) /(3+6 ) 
3(1+26) /2  .... 

t  as  the  different  stages  occurred. 

We  note  here  that,  in  fact,  transfer  must  be  due  to 
some  combination  of  the  ambient  and  mechanically  generated  turbulence. 
A  similarity  solution  does  exist  for  the  case  =  g1^3  by1  ^3  +  y  U*  but 
only  for  y  =  -4/3  which  gives  an  increasing  impulse!  It  is  clear  from 

the  above  equation  that  if  the  cell  size  increases  with  time  and  the  vel¬ 
ocity  reduces  (as  appears  to  be  the  case  for  a  descending  vortex  pair), 
then  the  ambient  and  mechanically  generated  turbulence  cannot  main¬ 
tain  the  same  ratio.  This  illustrates  that  similarity  solutions  can  only 
be  expected  when~one  form  of  transfer  dominates. 

3.  Dimensional  Analysis 

If  we  attempt  simply  to  model  the  problem  dimension- 

ally,  we  find  that  there  are  three  initial  parameters,  b  ,  U  =  r  /2rib 

o  o  o  o 

and  r  .  A  functional  equation  connecting  U  and  t  may  be  written  as 

u*/uo  *  fm<e>'3b'0'S/u.  •  t'bo'ro'> 


where  we  have  used  a  mechanical  time  scale  to  normalize  t  and  retain 
the  parameter  ^  ^ 3b  ^3Uq’  *  .  This  would  provide  a  collapsing  of 

data  where  ambient  turbulence  was  insignificant.  However,  where  am¬ 
bient  turbulence  dominated,  one  would  use  an  ambient  time  scale  and 

Write  U*/Uq  =  fa(el,3bi,s/U0  ,  t/^3c“1/3) 
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Thus  formulation  will  collapse  data  in  regions  where  the  ambient  turbu¬ 
lence  is  important.  A  representation  of  this  is  shown  in  Figure  4-3 
with  the  ambient/mechanical  ratio  1-5  as  a  parameter. 


1 


If  data  collapse  of  this  sort  is  observed  from  proper 
analysis  of  flight  tests,  then  the  exponents  of  t  described  in  the  previ¬ 
ous  sections  can  be  determined. 

We  note  in  passing  that  there  is  yet  another  time  scale, 
that  of  viscous  transfer.  According  to  Maxworthy's  experiments,  this 
has  the  form  (bQ^/U  It  appears  that  this  scale  is  too  large  to  be 

of  concern  to  us;  that  is,  the  time  involved  for  the  problem  to  become  a 
purely  viscous  one  is  beyond  our  range  of  interest. 
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4.3  NEUTRAL  ROTATIONAL  FIELD 


Statement  of  the  Problem 

Various  explanations  of  the  motion  of  vortices  in  a  croaswind 
with  shear  have  been  given  (Burnham,  1972;  Harvey  and  Perry, 

1971).  These  analyses  were  concerned  with  the  vortex  transport 
near  the  ground  plane,  where,  typically,  both  vertical  shear  and  ver¬ 
tical  shear  gradients  are  present  when  there  is  a  crosswind.  The  ob¬ 
jective  of  these  works  was  to  analytically  model  the  observed  unsym- 
metrical  motion  of  the  vortex  pair  as  they  descended  into  the  wind 
shear  region.  Instead  of  noth  vortices  leveling  off  at  the  expected 
height  above  the  ground  of  h  =  b/2  ,  the  upwind  vortex  sometimes  drops 
below  this  value,  and  the  downwind  vortex  levels  off  and  tnen  begins  to 
rise  to  a  distance  greater  than  b/2  .  Thus,  a  tilting  of  the  vortex  pair 
results,  superimposed  on  the  lateral  motion  of  the  vortices  caused  by 
the  uniform  component  of  the  crosswind,  and  the  proximity  of  the  ground 
plane. 

Such  analysis  near  the  ground  plane  is  unsteady  both  because  of 
the  redistribution  of  vorticity  in  the  shear  gradient  of  the  external 
flow  and  the  change  in  separation  of  the  aircraft  vortex  pair  due  to  the 
presence  of  the  ground  plane. 

However,  vortex  tilting  has  been  observed  at  sufficient  heights 
above  the  ground  (Tombach  1972),  that  ground  effect  need  not  be  con¬ 
sidered.  Under  these  conditions,  the  problem  can  be  treated  in  the 
classical  sense  of  a  descending  vortex  pair  in  steady  flow.  Also,  at 
altitude,  any  severe  vertical  crosswind  shear  graiients  would  not  be 
expected  to  be  present,  and  the  shear  can  be  represented  as  a  uniform 
vorticity  field,  again  a  steady  flow  situation.  The  problem  may  now  be 
examined  in  a  steady  coordinate  system,  fixed  to  the  vortex  pair.  The 


solution  to  this  problem  will  yield  the  streamlines  of  the  external  flow 
as  well  as  the  recirculating  flow  within  the  vortex  cells.  Since  the  vor¬ 
tex  tilting  phenomenon  has  been  observed  at  altitude,  such  an  analysis 
is  expected  to  give  a  good  inviscid  representation  of  the  effects  of  the 
flow  asymmetry  produced  by  a  sheared  crosswind  interacting  with  the 
vortex  pair  and  should  form  the  basis  for  a  more  detailed  viscous  solu¬ 
tion. 

Analysis 

A  vortex  pair  with  circulation  +PQ  are  located  at  y  =  0  and 
x  =  jlb/2  respectively.  A  vertical  component  of  velocity  v  =  r  /2n-b 
is  directed  upwardfrom  below  the  vortex  pair.  The  wind  shear  is  re¬ 
presented  as  a  horizontal  velocity,  u  =  f(y)  with  u  -  f(0)  =  0  to  remove 
the  effects  of  the  contribution  of  the  uniform  component  of  the  cross- 
wind  (for  the  steady  coordinate  system  chosen). 

Assuming  a  linear  crosswind  shear  profile,  u  =  f(y)  =  Ky  ,  from 
which  f(0)  =  0  .  This  satisfies  the  condition  of  the  removal  of  the  uni¬ 
form  component  of  the  crosswind.  The  details  of  the  combined  flows 
described  above  are  presented  in  Figure  4-4, 

The  stream  functions  for  each  of  the  flows  may  easily  be  written 
separately  as: 

vortex  flow: 


Yi  = 


4tt 


Sn 


(x-b/2)8  +  y3 
(x+b/2)3  +  y2 
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uniform  vertical  flow: 


*3 


and 


shear  flow: 


Y3 


from  which  the  stream  function  ^  for  the  complete  flow  field  may  be 
obtained. 


By  superposition: 


Y  =  Yi  +  +  Ya  , 


and 


r' 

o 

Z  TT  b  x 


4rr  [(x+b/2)3  +  y3J 


(4-1)) 


Introducing  the  following  dimensionless  parameters. 


X 


x 

"TJz 


Y 


a 


Krrb3 

~ZT~  ’ 

o 
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Eqn.  4-1  becomes: 


Y*  =  ay2 


x  - 


[(xtI)3  +  y8  j 


(4-2) 


Here,  a  is  a  non-dimensional  shear  parameter  which  relates  the  scale 
of  the  wind  shear  to  the  vortex  ; 
are  located  at  Y  -  0  ,  X  =  +  1  . 


of  the  wind  shear  to  the  vortex  parameters,  ro  and  b  .  The  vortices 


Before  proceeding  further  it  would  be  useful  to  calculate  the  order 
of  magnitude  of  the  shear  parameter,  q,  for  subsequent  evaluation  of 
the  streamlines,  ,  from  Eqn.  4-2. 

For  a  727  aircraft, 

T  =  34CG  ft2  /sec 
o 

b  =  108  ft 
o 

b  =  (rt/4)b  =  IT  (108)  /4  =  84.  78  ft 

o 

a  =  =  3.  31  K  . 

Typical  values  of  shear  near  the  ground  (Zwieback,  1964)  give  a  velocity 
of  18  ft/sec  at  44  ft  above  ground,  so  that  K  =  u/y  =  18/44  =  .4  sec  *  . 
Stronger  shears  a  :e  possible  closer  to  the  ground  or  in  gusty  situations, 
so  that  it  is  not  unreasonable  to  assume  a  range  of  0  sK  s  1  ,  This 
gives  an  upper  limit  on  a  of  o  -  3.  31  .  Eqn.  4-2  will  be  evaluated 
over  r.  range  of  a  from  zero  to  3. 
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Computation  of  y* 


Equation  4-2  gives  the  coordinates  X,  Y  for  the  family  of 
streamlines  Y*  at  constant  a  •  Although  the  equation  is  in  closed 
form,  it  is  not  easy  to  solve  for  X  and  Y  for  any  given  y#  ,  especi¬ 
ally  when  the  value  of  y*  ,  the  stagnation  point  streamline,  which  also 
describes  the  bounding  streamline  of  the  recirculating  vortex  cells,  is 
not  known  apriori. 

The  expression  for  y*  was  programmed  on  a  Hewlett-Packard 
Model  9820A  mini- computer  using  an  x-y  plotter  for  data  output.  After 
some  experimentation  with  various  inputs  of  y*  and  o  ,  the  pattern  of 
the  structure  of  the  streamlines  became  clear.  Figures  4'- 5,  -6,  -7, 
and  -8  show  the  geometry  of  the  stagnation,  or  dividing  streamline  of 
the  vortex  recirculation  cells  as  well  as  several  streamlines  of  the  ex¬ 
terior  flow  for  values  of  -  0,  1,2,  3  respectively.  In  addition,  Figure 
4-9  is  a  "magnified"  view  of  the  streamlines  in  the  immediate  vicin¬ 
ity  of  the  stagnation  point  for  the  case  a  =  3  . 


Discussion  of  Results 

The  obvious  (and  very  striking)  conclusion  from  Figures  4-5  through 
4-8  is  that  as  the  strength  of  the  wind  shear  increases,  the  size  of  the 
upwind  vortex  recirculation  cell  increases  and  the  size  of  the  downwind 
cell  decreases.  Presumably,  for  a  large  enough  (a  not  very  realistic 
case),  the  downwind  cell  would  approach  zero. 

The  problem  treated  here  was  assumed  steady  in  all  respects,  in¬ 
cluding  constant  a  ■  That  is,  the  vortex  pair  was  generated  in  a  cross- 
wind  shear  field  of  constant  o  which  remains  constant  as  the  vortex  pair 
translates  downwind.  It  is  not  clear  from  the  present  analysis  what  the 
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Figure  4-8.  Streamlines  in  the  Vicinity  of  a  Vortex  Pair  in  Feavy  Shear.  Value  of  the 
Shear  Parameter,  a  =  3.0.  Crosswind  from  the  Left. 


y 
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Figure  4-9.  Streamlines  in  the  Immediate  Vicinity  of  the  Top 
Stagnation  Point  for  a  Vortex  Pair  in  Heavy  Shear,  Value 
of  the  Shear  Parameter,  a  -  3.0. 
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effects  of  changing  n  would  be  on  the  geometry  of  the  flow,  once  es¬ 
tablished.  However,  it  might  be  reasonable  to  speculate  that  Figures 
4-5  through  4-8  give  an  instantaneous,  steady  view  of  the  flow  for 
each  value  of  a  ,  and  if  o  changes  slowly,  the  flow  field  would  change 
roughly  in  accordance  with  these  figures.  Thus,  as  o  increases,  the 
upwind  cell  would  gradually  increase  in  size  and  the  downwind  cell 
would  gradually  shrink. 

Further  analysis  must  be  performed  to  verify  this,  but  if  the 
above  assumption  is  reasonable,  a  method  of  predicting  the  flow  geo¬ 
metry  of  a  vortex  pair  descending  into  a  crossflow  with  shear  gradients 
would  be  possible. 

It  is  not  clear  from  the  present  analysis  of  the  exact  effects  that 
a  sheared  crosswind  has  on  vortex  tilting.  However,  it  is  unquestion¬ 
ably  the  case  that  the  wind  shear  produces  asymmetry  in  the  recircu¬ 
lating  cells  surrounding  the  two  vortices.  This  asymmetry  may  explain 
the  process  whereby  vortex  tilting  occurs.  It  should  be  remembered 
that  the  asymmetry  always  acts  to  increase  the  size  of  the  upwind  vor¬ 
tex  recirculation  cell  and  to  reduce  the  size  of  the  downwind  cell,  and 
that  this  effect  is  more  pronounced  as  the  strength  of  the  wind  shear  is 
increased. 

Section  4-2  describes  the  outward  diffusion  of  the  vortical  core 
as  the  age  of  the  vortex  pair  increases.  The  core  diffusion  is  presum¬ 
ably  related  to  vorticity  gradients.  Due  to  wind  shear,  the  downwind 
vortex  recirculation  cell  will  be  smaller  than  the  upwind  cell.  The  dif¬ 
fusing  vorticity  of  the  downwind  vortex  will  thus  reach  the  boundary  of 
the  inviscid  cell  before  the  similar  situation  occurs  for  the  upwind  vor¬ 
tex.  Thus,  vorticity  associated  with  the  downwind  vortex  will  diffuse 
into  the  free  stream  flow  and  be  swept  away,  reducing  the  vorticity  and 
hence  the  circulation  of  that  vortex.  An  asymmetry  will  then  be  produced 
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in  the  downwash  velocities  of  one  vortex  with  respect  to  the  other. 

Although  the  effect  of  these  asymmetrical  induction  velocities 
should  take  into  account  the  exact  distribution  of  vorticity  within  the  al¬ 
tered  vortical  cores,  the  sense  in  which  the  vortices  translate  with  re¬ 
spect  to  each  other  should  be  straightforward. 

If  the  circulation  of  the  downwind  vortex  is  decreased  over  that 
of  the  upwind  vortex  due  to  detrainment  of  vorticity,  then  its  downwash 
velocity  will  be  less  than  that  of  the  upwind  vortex.  As  a  result,  the 
upwind  vortex  will  translate  downward  less  rapidly  than  the  downwind 
vortex.  Tilting  will  occur  with  the  sense  that  that  the  upwind  vortex 
will  rise  with  respect  to  the  downwind  vortex. 

Thi3  result  contradicts  Burnham's  observation  that  "the  upvind 
vortex  usually  drops  to  a  lower  altitude  than  the  downwind  vortex.  " 
However,  Burnham's  work,  and  the  experimental  observations  which 
he  cites  were  all  performed  in  ground  effect  at  relatively  low  cross- 
wind  velocities.  He  states,  "Both  experimental  observations  and 
calculations  indicate  that  other  phenomena  occur  at  much  higher  cross- 
wind  speeds  (e.g.  ,  30  fps  instead  of  the  10  1'ps  considered  here).  " 

Future  work  needs  to  be  performed  to  adequately  describe  the 
me<  •  anism  which  initiates  vortex  tilting.  It  seems  reasonable  from 
the  p.  esent  analysis  that  the  asymmetry  in  the  flow  fields  surrounding 
the  vortex  cores  forms  an  excellent  point  of  departure  for  such  further 
investigation.  The  assumption  of  a  uniform  shear  was  sufficient  to  in¬ 
dicate  this  asymmetrical  behavior.  Including  the  effects  of  shear  gradi¬ 
ents  in  future  work  would  undoubtedly  refine  the  results  somewhat,  but 
their  effects  should  probably  be  small  and  would  produce  results  in  the 
same  sense  indicated  by  the  present  analysis. 


For  a  sufficiently  large  vortex  age,  the  core  diffusion  process 
would  continue  to  a  point  at  which  vorticity  detrainment  from  the  down¬ 
wind  vortex  (with  shear  flow)  would  result  in  sufficient  reduction  of  its 
circulation  that  for  all  practical  purposes,  it  would  cease  to  exist. 

This  would  not  be  the  case  for  the  upwind  vortex,  considering  its  larger 
cell  size.  This  process  would  proceed  more  rapidly  for  larger  3  .  In¬ 
deed,  for  a  locally  large  value  of  a  such  as  might  exist  in  a  gust  with 
high  shear  rate,  the  above  process  would  occur  catastrophically  --  a 
possible  explanation  for  core  bursting.  For  both  the  case  of  the  slower 
detrainment  rate  and  the  catast.  ophic  situat.on,  the  remaining  vortex 
would  be  the  upwind  vortex. 

Concluding  Remarks 

The  calculation  of  the  family  of  streamlines  y*  for  a  given  value 

of  shear  3  is  tedious.  The  only  really  necessary  streamline  to  obtain 

is  that  which  describes  the  boundary  of  the  recirculation  cells,  yg  , 

However,  the  value  of  yg  is  identical  to  the  value  of  the  stream  function 

at  the  stagnation  point  y*  .  y*  may  be  solved  in  terms  of  3  explicitly 

s  s 

and  presented  graphically  or  in  tabular  form.  Substitution  of  y*  =  Y*(a) 

s 

in  Eqn.  (4-ZJ  will  then  give  the  shape  of  the  recirculation  cells  for 
any  desired  value  of  3  . 

The  conditions  for  the  calculation  of  y*  are: 

s 

**  “  MXg.Ys.o)  -  0  ,  (4-3) 

=  fa(X  ,Y  .0)  =  0  ,  (4-4) 

9  S 
S 

=  fu(Xg,Yg)  =  0  ,  (4-5) 

s 


where  the  subscript  s  implies  that  the  quantities  are  to  be  evaluated 
at  the  otagnation  point.  Equation  (4-3)  is  identical  to  Eqn.  (IV-  2) 
evaluated  at  s  . 


The  simultaneous  solution  of  Eqns.  (4  -  3,  4,  5) 


will 


yield  the  values  of  the  three  unknowns,  X  ,  Y  ,  y*  in  terms  of  p 

8  S  3  v 


of  which  only  y#  is  of  interest.  The  simultaneous  solution  of  these 


equations  is  also  not  exactly  easy,  but  can  be  performed  in  a  straight¬ 
forward  manner  using  computer  techniques. 


The  cell  boundary  streamline  will  be  given  by  the  expression: 

1  -  o  ,  (4-6) 


-y*(a}  +  a Y2  -  X  -  ft!  fe-Ul-i  1 


(x+1)3  +  y8  J 


from  which  the  shape  of  the  recirculating  cells  may  be  calculated  or 
plotted  by  computer.  The  shape,  sir.e,  or  area  of  the  recirculating 
cells  may  thus  be  obtained  as  continuous  functions  of  the  wind  shear, 
r;  ,  within  the  accuracy  of  the  computational  methods  used. 
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4.4  WAKE  TRANSPORT  IN  A  STRATIFIED  FIELD 


If  a  -'ortex  wake  descends  through  a  stably  stratified  atmos¬ 
phere,  as  is  usually  the  case,  then  it  becomes  buoyant  as  it  descends 
because  its  density  increases  at  a  slower  (adiabatic)  rate  than  the  am¬ 
bient  density.  A  wake  could  also  conceivably  acquire  buoyancy  at  the 
time  of  its  generation,  from  hot  engine  exhaust  gases  mixed  into  it.  In 
either  case,  the  buoyancy  thus  generated  acts  to  decrease  the  overall 
wake  impulse  and  (possibly)  circulation,  although  any  mixing  between 
the  wake  and  its  surroundings  would  tend  to  dilute  the  buoyancy  and  to 
decrease  these  effects. 


4.4.1  EXPERIMENTAL  OBSERVATIONS 

Observations  of  wakes  generated  by  full-size  aircraft 
tend  to  suggest  that  stratification  may  have  an  effect  on  wake  transport. 
Specifically,  a  study  of  the  dissemination  of  particles  released  from  air¬ 
craft  (Smith  and  Wolf,  1963;  Smith  and  MacCready,  1963)  included  some 
observations  of  wake  descent  from  a  variety  of  aircraft.  Tests  with  air¬ 
craft  up  to  DC-3  size  gave  typical  descent  distances  of  25  to  100  feet  over 
land,  with  the  wake  descending  initially  at  the  theoretical  speed,  but  then 
broadening  and  slowing  up.  On  the  other  hand,  similar  measurements  over 
the  ocean  sometimes,  but  not  always,  indicated  wake  descents  of  600 
feet.  A  possible  explanation  was  that  turbulence  slowed  the  wake's 
motion  and  helped  spread  it  and  that  a  stable  atmosphere  further  im¬ 
peded  descent.  Both  stability  and  turbulence  are  usually  less  over  the 
ocean  than  over  land,  thus  resulting  in  longer  vortex  life  and  descent. 


Another  study  (FAA  Task  Force,  1971)  found  that  the 
vortices  generated  by  jumbo  jet  aircraft  began  descending  at  400  to  500 
feet  per  minute,  but  that  fully  developed  wakes  were  generally  not  found 
more  than  1000  feet  below  the  altitude  of  the  generating  aircraft. 
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Leveling  off,  combined  with  start  of  breakup  usually  began  after  a  de¬ 
scent  of  800  to  900  feet  and  the  vortex  spacing  appeared  to  remain  con¬ 
stant  until  breakup.  Occasional  descents  of  1100-1200  fee'c,  with  no 
noticeable  slowing  of  descent  speed,  were  measured  during  the  same 
study  by  Andrews,  et  al  (1972).  In  contrast,  in  one  case  the  wake  from 
a  cruising  C-5A  (Mach  no.  =  0.8)  was  measured  at  500  feet  below  the  air¬ 
craft  at  3  miles  behind  it,  and  remained  at  that  level  for  more  than  40 
miles. 

During  a  British  study  of  wakes  behind  jet  transport 
aircraft  (Rose  and  Dee,  1965)  wake  descents  of  up  to  800  feet  were  meas- 
sured  during  the  first  150  seconds  of  wake  life  and  no  general  leveling-off 
trend  was  observed.  Quantitative  measurements  of  vortex  spacing  were 
obtained  in  another  British  study  (Bisgood,  et  al,  1971)  where  the  vortex 
spacing  behind  a  delta-wing  aircraft  was  observed  to  grow  to  about  3 
times  the  initial  spacing  in  about  45  seconds. 

The  descent  of  the  wake  generated  by  a  lightplane  in  a  variety  of  meteo¬ 
rological  conditions  was  measured  by  Tombach  (1972,  1973).  He  found 
that  the  descent  speed  of  the  wake  was  constant  for  about  20  seconds  in 
a  stable  atmosphere,  after  which  the  descent  was  slowed  and  often 
stopped.  There  seemed  to  be  no  noticeable  differences  in  wake  behavior 
at  different  atmospheric  stratifications,  however,  but  the  dominant  fact¬ 
or  governing  wake  descent  was  probably  a  banking  of  the  plane  of  the  vor¬ 
tices  (which  was  discussed  earlier). 

Tombach  also  attempted  to  measure  vortex  spacing 
as  the  wake  descended.  The  same  banking  tendency  complicated  the 
interpretation  of  the  data,  but  he  noted  that  slight  increases  in  vortex 
spacing  could  be  observed  in  the  few  cases  when  the  wake  remained 
relatively  level.  Considerable  increases  in  spacing  were  noted  for 
wakes  which  had  banked  well  out  of  the  horizontal  plane. 


4-34 


The  observations  by  Tombach  showed  that  the  smoke- 
marked  vortex  wake  sometimes  rose  again  after  its  descent  had  stopped. 
Whether  the  cause  for  this  is  atmospheric  upcurrents  or  buoyancy  is  not 
known. 

A  few  laboratory  experiments  have  been  performed  to 
study  the  motion  of  buoyant  vortex  pairs.  Most  of  these  have  been  con¬ 
cerned  with  line  thermals,  where  the  buoyancy  of  the  thermal  is  such  as 
to  increase  the  total  circulation  around  each  vortex.  When  in  a  homo¬ 
geneous  medium  and  well  developed,  i.  e.  ,  when  the  vortex  motion  is 

well  organized,  these  thermals  are  observed  (Richards,  1963;  Tsang, 

1/3  2/3 

1971)  to  grow  in  scale  as  A  p  t  ,  where  t  is  time  (with  t  =  0  at 

the  extrapolated  point  when  the  vortex  pair  would  have  been  of  zero  size) 

and  A  p  is  the  density  difference  between  the  thermal  and  the  ambient 

fluid  at  the  time  of  release.  Their  speed  of  advance  has  been  found  to 
1/3  -1/3 

slow  as  A  p  t  .  This  entrainment  of  the  ambient  fluid  is  observ¬ 
ed  to  erode  the  buoyant  driving  force  and  to  eventually  bring  the  growing 
thermal  to  a  halt.  Extension  of  such  results  to  the  vortex  wake  problem 
is  not  straightforward,  however,  because  (1)  the  buoyancy  opposes  the 
motion  of  the  vortex  wake,  and  (2)  the  impulse  of  the  thermal  is  coupled 
to  its  initial  buoyancy,  whereas  the  impulse  of  the  vortex  wake  is  totally 
unrelated  to  its  initial  buoyancy  (which  is  often  zero). 

One  set  of  laboratory  experiments  to  study  the  motion 
of  a  vortex  pair  with  initial  impulse  in  a  stratified  medium  was  perform¬ 
ed  by  Tulin  and  Shwartz  (1971)  (with  a  related  experiment  by  Birkhead, 
et  al,  1969).  They  generated  an  impulsive  vortex  pair  in  a  stratified 
water  tank  and  observed  its  motion,  which  looked  like  that  shown  in 
Figure  4-10.  (The  figure  has  been  redrawn  to  correspond  to  the  di¬ 
rection  of  motion  of  an  aircraft  wake).  The  total  vertical  travel  of  the 
wake  was  about  4  vortex  spacings  (which  is  of  comparable  magnitude  to 
the  aircraft  results  quoted  earlier)  before  the  stratification  dominated 
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Figure  4-10.  Trajectory  of  a  Vortex  Pair  in  a  Linearly  Density 
Stratified  Medium  (from  Tulin  and  Shwartz,  1971).  The 
Upper  Sketches  Show  the  Observed  Cross  Section  of  the 
Wake. 
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the  motion  and  caused  the  wake  to  reverse  direction  and  collapse.  This 
process  took  place  in  a  characteristic  time  which  is  proportional  to  the 
Vaisala  period  =  (g/pdp/dy  -  g  /a  )  where  p  is  the  atmospheric 
density  and  a  is  the  speed  of  sound,  with  the  proportionality  factor  de¬ 
pending  on  the  initial  strength  of  the  vortices,  which  in  turn  is  repre¬ 
sented  by  another  time  =  b^/V  where  V  is  the  descent  speed  and 
b  is  the  vortex  spacing.  According  to  Tulin  and  Shwartz,  if  T-  <  T  , 
the  vortices  dominate  the  motion,  while  if  T^  >  T^  the  stratification 
dominates  it.  Since  T^  increases  if  the  wake  descends  and  grows, 
their  experiments  show  that  stratification  must  invariably  dominate  the 
motion  unless  some  instability  has  first  destroyed  the  vortices. 

The  Tulin  and  Shwartz  data  appear  to  be  of  good  quality, 
but  some  reservations  exist  about  their  applicability  to  the  aircraft  wake 
problem.  The  method  of  vortex  generation,  with  an  impulsive  motion 
of  a  piston  pushing  fluid  through  a  slit,  does  not  generate  a  pure  vortex 
pair,  but  adds  a  wake  behind  the  moving  vortex  pair.  (This  problem  is 
a  well  known  one  in  vortex  ring  generation  also).  Since  the  total  vortex 
wake  motion  from  the  wall  orifice  is  only  a  few  times  the  wake  size,  the 
effects  of  the  wall  and  of  the  wake  trailing  behind  the  vortices  might  be  sig¬ 
nificant.  It  should  be  noted  also,  that  the  stratification  was  quite  strong, 
with  the  Vaisala  period  T  ranging  from  1-4  seconds  (as  compared  to 
typical  atmospheric  values  around  100  sec.).  However,  the  ratio  Ta/T 
is  initially  of  order  1/10  for  both  their  experiments  and  for  full  scale 
wakes  in  the  atmosphere,  which  indicates  comparable  scaling  of  vortex 
strength  versus  buoyancy. 

In  summary,  the  expei  .menial  picture  is  somewhat 
fuzzy.  The  evidence  to  show  that  wake  descent  stops  in  a  stratified  at¬ 
mosphere  is  not  conclusive,  but  the  suggestion  is  there.  Similarly, 
vortex  spacing  seems  to  vary  little  with  descent,  or  possibly  to  increase, 
but  again  the  experimental  evidence  is  not  extremely  strong.  Using  the 
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knowledge  of  wake  behavior,  limited  and  confused  as  it  may  be  as  a 
foundation,  let  us  now  look  at  efforts  to  approach  the  problem  analy 
ticaliy. 


4.4.2  PREVIOUS  THEORETICAL  MODELS 


The  theories  which  have  been  developed  to  describe  the  motion 
of  a  vortex  wake  in  a  stably  stratified  atmosphere  are  numerous, 
include  various  types  of  modeling  approaches,  and  give  a 
variety  o»  contradictory  conclusions.  The  state-of-the-art  is  sum¬ 
marized  below  in  Table  4-2.  As  shown  in  the  table,  the  models  sug¬ 
gest  two  main  types  of  behavior  --  a  slowing  of  wake  descent  with  an 
increase  of  vortex  spacing,  or  an  acceleration  of  wake  descent  accom¬ 
panied  by  a  decrease  in  vortex  spacing. 

TABLE  4-2.  COMPARISON  OF  THEORETICAL  MODELS 
FOR  DESCENT  OF  A  VORTEX  WAKE 
IN  A  STA3LY  STRATIFIED  ATMOSPHERE. 


Buoyancy- 

Vortex 

Descent 

generated 

Author 

Spacing 

Speed 

Vorticity 

Costen  (1972)** 


Weak  stability 

Decreases 

Increases 

Nro  effect 

Strong  stability 

increases  (? )* 

Stops  (°)* 

Kuhn  Is  Nielsen  (1972) 

Decreases 

Increases 

Pa  rtly 
entrained 

Saffman  (1972) 

Increases 

Stops 

Ent  rained 

Scorer  &  Davenport  (!97i 

0)  Decreases 

increases 

Detrained 

Tombach  (1971) 

Strong  stability 

Inc  rrases* 

Stops® 

Entrained 

Weak  stability 

Dec  reases 

Increases 

Entrained 

Tultn  k  Shwarts  { 1  97 1 ) 

Increases 

Stops 

Entrained  k 
annihilated 

*  See  discussion  in  text. 

**Costcn's  model  is  a  bu 
a  buoyant  wake  oval. 

oy»«t  core  model. 
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Ail  others  are 

models  for 

In  order  to  assess  the  credibility  of  the  various  models  each  of 
them  will  be  discussed  briefly  below.  Since  there  has  been  a  slight 
evolutionary  trend,  they  will  be  presented  in  a  generally  chronological 
order. 


The  first  two  modeling  approaches  (Scorer  and  Davenport,  1970; 
Tombach,  1971)  are  both  based  on  the  assumption  that  the  effect  of 
buoyancy  is  to  change  the  classical  hydrodynamic  impulse  of  the  vor¬ 
tex  pair  and  its  accompanying  fluid.  Both  assume  that  the  buoyancy  ef¬ 
fects  take  place  rather  slowly,  so  that  the  instantaneous  shape  of  the 
oval  cylinder  of  accompanying  fluid  is  the  same  as  that  of  the  classical 
wake  in  a  uniform  medium,  and  hence  the  wake  size  scales  with  the  vor¬ 
tex  spacing.  The  basic  equations  of  motion  are  thus 


and 


dl 

dt 


PoF 


dF 

dt 


-Av  /T8 


where  I=»  2  poTh  is  the  impulse  per  unit  length  of  wake,  pQ  is  the 
initial  density,  F  =  gAAp/oQ  gives  the  buoyancy,  A  is  the  cross- 
sectional  area  of  the  wake  oval,  the  Vaisala  period  T  gives  the  stabil¬ 
ity,  and  v  is  the  descent  speed  of  the  wake. 

Scorer  and  Davenport  assume  the  circulation  to  be  constant.  They 
find  a  single  solution  in  which  the  vortices  converge  and  accelerate  downward 
in  a  stable  atmosphere.  They  compute  numerically  the  internal  streamline 
patterns  for  this  case  and  suggest  that  the  detrainment  is  indeed  a  stable 
process  and  that  the  circulation  is,  in  fact,  constant  because  the  vorticity 
generated  by  the  buoyancy  at  the  boundary  of  the  wake  oval  is  continually 
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detrained  as  a  curtain  above  the  descending  oval.  After  a  time,  how¬ 
ever,  they  state  that  some  of  this  fluid  and  vorticity  is  mixed  into  the 
wake,  which  results  in  its  eventual  destabilization  and  destruction. 

The  idea  of  detrainment  has  some  experimental  support.  Photo¬ 
graphs  of  contrails  made  from  the  side  (Smith  and  Beesmer,  1959)  show 
a  curtain  of  condensed  vapor  extending  above  the  wake  up  to  the  flight 
level  of  the  generating  aircraft.  These  same  contrail  observations,  as 
well  as  low  altitude  measurements  by  Smith  and  MaeCready  (1963), 
suggest  however  that  the  vortex  separation  increases  at  least  slightly 
as  the  wake  descends  instead  of  decreasing  as  the  theory  predicts. 

This  troublesome  point  will  reappear  again  in  Uter  analyses,  so  it 
should  be  pointed  out  here  that  buoyancy  may  indeed  cause  the  vortices 
to  converge  and  accelerate  even  though  such  behavior  has  not  been  ob¬ 
served  experimentally.  Since  none  of  the  analyses  consider  either  tur¬ 
bulence  or  viscosity,  both  of  which  are  omnipresent  in  the  atmosphere, 
it  is  quite  possible  that  other  mechanisms  overwhelm  the  buoyancy  ef¬ 
fect  and  cause  She  observed  behavior. 

Saffman  (1972)  points  out  that  the  Scorer  and  Davenport  analysis 
does  not  properly  treat  the  dynamical  effect  of  the  vorticity  generated 
at  the  interface  between  the  accompanying  fluid  and  the  denser  surround¬ 
ing  atmosphere.  A*  he  shows,  this  vorticity  must  be  included  in  the 
calculation  of  the  impulse  of  the  wake.  The  error  in  neglecting  it  may 
be  small,  however,  if  the  generated  vorticity  is  totally  annihilated  on 
the  centerline  by  interaction  with  its  counterpart  of  opposite  sign  from 
the  other  side  of  the  oval.  This  is  equally  true  for  entrained  or  de¬ 
trained  vorticity. 

There  arc  several  confusing  mathematical  error*  in  the  Scorer 
and  Davenport  paper  which  affect  its  readability,  but  do  not  alter  any 
of  its  conclusions.  Specifically,  the  right  sides  of  their  equations  ($j 
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which 


and  (11)  should  both  be  multiplied  by  (-1).  The  quantity  n  , 
they  call  the  vorticity,  is  really  the  vortex  sheet  strength,  i.  e.  ,  the 
circulation  per  unit  length  of  sheet  (which  is  the  vorticity  integrated 
across  the  sheet  thickness).  Their  equations  (17)  for  the  "non-dimen¬ 
sional"  vorticity"  CD  do  not  give  a  dimensionless  value.  The  corrected 
expressions  should  be 

5cp  _  k  w  5  cp  _  k  _1_ 

~  4ttR  uq  ’  bC  4ttR  q 


If  these  equations  are  used,  then  the  expression  in  Ean.  (16)  is  correct. 

In  contrast  to  the  constant  circulation  postulated  by  Scorer  and 
Davenport,  Tombach  (1971)  considered  the  possibility  that  buoyancy¬ 
generated  vorticity  would  decrease  the  overall  circulation  around  each 
vortex,  by  allowing  their  circulation  to  vary  at  a  rate  proportional  to 
the  ambient  density  gradient.  The  approach  used  followed  that  of  Turner 
(i960)  for  rising  buoyant  vortex  pairs.  In  addition  to  the  geometric  simi¬ 
larity  assumption  made  by  Scorer  and  Davenport,  Tombach  also  assumed 
that  the  overall  entrainment  and  mixing  process  could  be  characterized 
by  a  single  parameter  which  is  represented  as  a  length  in  the  Bjerkness 
circulation  equation,  with  the  rate  of  change  of  the  circulation  of  a  vor¬ 
tex  being  proportional  to  this  length.  He  also  assumed  that  this  length 
scales  with  all  the  other  lengths  in  the  wake,  which  then  results  in  the 
rate  of  change  of  circulation  being  related  to  the  w^h.e  scale  and  the 
atmospheric  stability,  and  suppresses  any  explicit  display  of  the  effect 
of  the  wake  buoyancy.  Although  a  similar  assumption  was  used  by 
Turner  (I960)  and  is  supported  by  experiments  by  Woodward  (1959) 
with  thermals,  it  is  not  a  valid  one  to  make  initially  when  the  vortex 
wake  has  no  buoyancy  and  thus  is  not  generating  any  vorticity.  Thus 
the  Tombach  model  does  not  properly  describe  either  the  initial  be¬ 
havior  of  the  wake  nor  can  it  properly  treat  it  at  very  long  times  when 
the  oval  height  could  be  sufficiently  great  so  that  gravitational  effects 
on  its  shape  have  to  be  considered. 
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In  the  intermediate  time  intervals,  the  model  suggests  that  two 
distinct  types  of  wake  behavior  may  be  possible,  with  the  choice  of 
which  behavior  will  occur  being  determined  by  the  initial  vortex 
strength,  the  atmospheric  stability,  and  the  degree  and  manner  of 
entrainment  (which  at  present  cannot  be  quantified).  If  the  initial  vor¬ 
tex  strength  is  weak  enough  and/or  the  stability  great  enough,  the 
model  says  buoyancy  rapidly  erodes  the  circulation,  and  the  residual 
impulse  can  only  be  accommodated  by  an  increase  in  the  vortex  separ¬ 
ation  and  a  slowing  down  of  the  wake  descent.  The  motion  is  much  like 
that  of  a  vortex  pair  descending  into  ground  effect. 

On  the  other  hand,  if  the  vortex  is  strong  or  the  stability  weak, 
then  there  is  residual  circulation  remaining  when  the  impulse  has  been 
eliminated  by  the  buoyancy,  and  the  kinematical  consequence  is  a  rapid 
convergence  of  the  vortices  with  an  accompanying  downward  acceler¬ 
ation,  similar  to  that  predicted  by  Scorer  and  Davenport.  The  index 
determining  which  behavior  will  occur,  denoted  by  Q  ,  is  given  by 

Q  ~  TT  /b  3 
o  o 

where  T  is  the  Vaisala  period,  rQ  is  the  initial  circulation  of  one 
vortex,  and  bQ  is  the  initial  vortex  spacing.  The  proportionality 
factor  depends  on  the  entrainment  details.  Small  Q  (< rr/2)  corres¬ 
ponds  to  the  first  type  of  motion  mentioned,  while  large  Q  (>rr/ 2)  cor¬ 
responds  to  the  second  type. 

As  one  would  expect,  the  time  scale  of  the  motion  is  found  to  be 

proportional  to  T  ,  hence  the  greater  the  stability  the  more  rapidly 

the  motions  described  above  take  place.  The  vertical  length  scale 

2/3 

is  proportional  to  T  ,  thus  as  the  atmospheric  stability  is  in¬ 
creased  the  vertical  extent  of  the  motion  will  decrease. 


4-43 


In  the  same  way  as  for  the  Scorer  and  Davenport  model,  Saffman 
states  that  the  Tombach  model  is  also  in  error  because  it  equates  the 
change  in  the  wake  impulse  to  the  buoyancy  force,  which  is  correct 
only  if  the  Boussinesq  approximation  holds  (i.  e.  ,  if  the  density  differ¬ 
ence  between  the  wake  and  the  environment  are  small)  and  if  all  the 
vorticity  generated  by  the  buoyancy  is  included  in  the  computation  of 
the  impulse.  He  claims  that  the  latter  requirement  was  not  satisfied 
in  the  analysis.  However,  the  Tombach  analysis  does  take  into  account 
the  change  in  impulse  due  to  the  buoyancy  generated  vorticity  according 
to  the  Bjerkness  equation.  It  does  assume,  however,  that  it  is  possible 
to  approximate  the  impulse  of  the  wake,  including  the  new  vorticity,  by 
the  impulse  of  a  geometrically  identical  wake  with  all  of  the  vorticity, 
old  and  new,  concentrated  in  the  two  vortex  cores.  This  approximation 
will  probably  result  in  errors  after  the  wake  has  acquired  considerable 
buoyancy,  although  by  that  time  the  model  may  no  longer  be  valid  for  a 
variety  of  other  reasons. 

Although  qualitatively  similar,  the  Tombach  iarge-Q  solution 
car.not  be  matched  to  the  Score-  and  Davenport  solution.  Scorer  and 
Davenport  assume  constant  circulation,  but  permit  detrainment  of  mass 
from  the  wake  oval.  On  the  other  hand,  Tombach  requires  that  any  m?<s 
transfer  must  affect  the  circulation,  hence  the  assumption  of  constant 
circulation  (which  requires  that  his  entrainment  parameter  s  =  0)  would 
automatically  require  no  entrainment  or  detrainment. 

Costen  (1972)  points  out  that  the  portion  of  the  Tombach  model 
for  which  the  circulation  decreases  more  rapidly  than  the  impulse  is 
impossible.  This  conclusion  is  derived  by  taking  the  ratio  of  Tom- 
bach's  equations  (3)  and  (9),  which  give  (in  his  notation) 

d  log  T  _  1  dT  .  A  dM  -x  _  2h  .  p-p" 
d  log  M  "  F  dt  *  M  dt  q  p-p' 
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where  hb  is  the  height  of  the  wake  oval  (b  is  the  vortex  semi-span 

2 

for  his  analysis),  the  wake  cross-sectional  area  is  qb  ,  p  is  the 
ambient  density,  p'  is  the  average  wake  density,  and  p"  is  the  av¬ 
erage  wake  density  along  the  axis  of  symmetry  between  the  vortices. 
Numerical  values  of  the  constants  are  h  =  3.46  and  q  «  11.  62  ,  giving 


4420, »  .  6 


d  log  M 


P-P 


and  since  one  would  expect  pa  p1!  a  p‘  ,  this  gives 


d  lQg  L  *  .6 

d  log  M 


Thus  the  only  self-consistent  solutions  of  Tombach's  equations  are 
those  for  which  Q  <  rr/2  ,  i.  e. ,  the  converging  and  accelerating  solu¬ 
tions. 


Another  theory  was  developed  by  Tulin  and  Shwartz  (1971)  in 
which  they  modeled  the  vortex  system  by  separate  velocity  scaling  of 
the  flows  internal  and  external  to  the  wake  oval,  with  a  shear  layer  at 
the  boundary  between  the  two  flow  fields.  Based  on  experiments  they 
performed  in  a  water  tank,  they  conclude  that  the  wake  entrains  the 
vorticity  generated  at  this  shear  layer  (in  contrast  to  the  detrainment 
postu  <ted  by  Scorer  and  Davenport)  but  that  the  ingested  vorticity  is 
mainly  canceled  out  at  the  wake  centerline  through  mixing  with  vorticity 
of  the  opposite  sign  from  the  other  side  of  the  wake.  Consequently, 
they  model  a  turbulent  wake  in  a  stratified  medium  by  assuming  con¬ 
servation  of  volume,  mass,  and  energy  and  neglecting  vorticity  and 
momentum.  They  then  find  completely  similar  solutions  which  depend 
on  four  parameters  --  the  initial  buoyancy  of  the  wake,  the  stability  of 
the  fluid,  the  dissipation  of  kinetic  energy,  and  the  ratio  of  kinetic  to 
potential  energy.  The  latter  two  parameters  are  assigned  values  based 
on  their  experiments,  which  were  dircussed  earlier. 
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Tulin  and  Shwartz  give  formulas  for  vortex  wake  motion  in  a 
stratified  medium,  which  are  too  complex  to  present  here,  and  present 
excellent  correlations  between  their  formulas  and  experiments.  Their 
model  predicts  that  a  wake  will  slow  down  and  stop  its  descent  in  a 
stratified  fluid  and  that  the  vortex  separation  will  increase  as  the  wake 
descends.  The  descant  in  a  homogeneous  fluid  has  the  behavior 
z  ~  t"  ^  and  b~  ^  .  In  a  stably  stratified  fluid,  the  descent  slows 
down  more  quickly  and  the  spreading  is  more  rapid,  but  similarity  no 
longer  holds  when  the  wake  has  come  to  rest.  The  nature  of  the  mo¬ 
tions  they  observed  experimentally  and  their  model  endeavors  to  de¬ 
scribe  was  shown  earlier  in  Figure  4-10. 

In  the  Tulin  and  Shwartz  formulation,  as  well  as  in  those  dis¬ 
cussed  above,  the  time  scale  of  the  motion  is  proportional  to  the 
Vaisala  period  T  .  For  those  models  which  predict  a  stopping  of  wake 
descent,  this  then  means  that  the  time  required  for  the  wake  to  descend 
to  its  lowest  level  is  proportional  to  T  . 

Saffman  (1972)  has  presented  a  model  which  indicates  that  a  vor¬ 
tex  system  in  a  stable  environment  descends  with  the  vortex  spacing 
remaining  constant  if  there  is  no  entrainment,  and  increasing  if  en¬ 
trainment  is  assumed  proportional  to  the  density  difference  between 
the  wake  and  environment.  In  either  case,  the  wake  descent  is  stopped 
by  the  stability  and  it  subsequently  oscillates  vertically.  To  find  these 
solutions,  Saffman  solves  the  Laplace  equation  for  the  potentials  inside 
and  outside  the  wake  volume,  with  th..  constraint  that  there  is  no  change 
in  wake  volume  (except  due  to  entrainment,  which  he  treats  separately). 
His  criticisms  of  the  Scorer  and  Davenport,  and  Tombach,  models  have 
already  been  reviewed,  and  his  formulation  was  designed  to  overcome 
the  deficiencies  he  described.  . 


Saffman  points  out  that  it  is  "by  no  means  certain"  that  turbulent 
mixing,  caused  by  the  Kelvin-Helmholtz  instability  of  the  wake-atmosphere 


interface,  must  occur  when  the  motion  is  generated  aerodynamically 
and  buoyancy  effects  are  (at  least  initially)  a  small  perturbation.  As 
was  pointed  out  earlier,  this  is  a  different  situation  from  the  line 
thermal  case,  where  the  motion  is  generated  by  buoyancy.  Since  the 
experimental  evidence  obtained  to  date  does  allow  not  allow  one  to  assert 
with  any  confidence  that  significant  turbulent  mixing  does  take  place,  he 

considered  both  non- entraining  and  entraining  cases  for  completeness. 

Work  by  Kuhn  and  Nielsen  (1972)  attempts  to  model  the  entrain¬ 
ment  differently  than  Tombach  or  Scorer  and  Davenport.  One  short¬ 
coming  of  the  Tombach  model  was  the  assumption  that  the  rate  of 
change  of  circulation  was  proportional  to  the  scale  of  the  wake.  This 
is  erroneous  at  small  times  if  the  wake  has  no  initial  buoyancy,  since 
then  one  should  have  no  generation  of  vorticity,  Kuhn  and  Nielsen  use 
a  variable  entrainment  length,  proportional  to  the  density  difference 
between  the  entrained  fluid  and  the  ambient  fluid,  with  the  entrained 
fluid  density  being  a  weighted  average  of  the  ambient  density,  the  den¬ 
sity  in  the  vicinity  of  the  vortex,  and  the  density  in  the  region  into 
which  the  outer  fluid  is  entrained.  The  effect  is  to  make  Tombach1  s 
entrainment  parameter  s  into  a  variable  which  has  value  zero  initially. 

Kuhn  and  Nielson  thus  model  the  entrainment  so  that  part  of  the 
buoyancy  generated  can  be  entrained  and  part  of  it  detrained,  with  the 
proportions  of  each  being  governed  by  an  unknown  wake  mixing  para¬ 
meter.  Their  analysis  shows  the  wake  accelerates  as  it  descends  and 
the  vortex  spacing  decreases.  Increased  values  of  the  wake  mixing 
parameter  reduce  the  rates  of  descent  and  convergence,  while  the  ad¬ 
dition  of  heated  air  to  the  wake,  but  outside  the  vortices,  is  found  to 
cause  a  leveling-off  and  spreading  of  the  wake.  Saffman's  criticisms 
about  the  proper  formulation  of  the  impulse  to  include  the  entrained 
vorticity  apply  here  also. 


Yet  another  model  by  Costen  (1972)  finds  that  buoyant  vortex 
cores  (different  from  the  buoyant  wake  considered  in  the  other  models) 
will  accelerate  and  close  together  as  they  descend  in  a  neutral  atmos¬ 
phere.  He  concludes  that  stable  atmospheric  stratification  will  accel¬ 
erate  the  effect.  Costen  also  suggests  that  a  strong  inversion  layer 
might  retard  the  vortex  motion  and  cause  their  spacing  to  increase, 
but  does  not  present  any  support  for  this  inference. 

In  summary,  the  analytical  situation  with  respect  to  the  effects 
of  stratification  is  not  clearer  than  the  experimental  one.  In  fact,  if 
one  arrives  at  a  consensus  based  solely  on  the  number  of  models 
which  predict  the  same  behavior,  one  would  conclude  that  buoyancy 
causes  the  wake  to  accelerate  downward  and  the  vortex  spacing  to  de¬ 
crease.  This  appears  to  be  in  contradiction  to  the  experimental  ob¬ 
servations,  which  suggest  the  opposite  behavior.  Both  the  theories  and 
the  experiments  may  be  correct  however,  with  the  differences  being  due 
to  factors  such  as  entrainment  (due  to,  say,  ambient  turbulence)  which 
were  not  considered  in  the  models.  In  fact,  all  of  the  models  which 
predict  the  slowing  down  and  spreading  out  behavior  do  include  entrain¬ 
ment  as  a  factor. 

There  are  fundamental  questions  which  arise  about  the  formula¬ 
tions  of  these  models,  however,  some  of  which  were  raised  during  the 
discussion.  All  of  the  models  contain  simplifying  assumptions  about 
the  distribution  of  vorticity  and  of  buoyancy,  in  order  to  make  a  tract¬ 
able  problem.  In  many  cases  the  effect  of  such  assumptions  is  not 
really  known,  nor  is  there  empirical  evidence  to  support  or  refute 
their  validity. 

As  part  of  the  current  study,  AV  has  developed  a  more  rigor¬ 
ous  description  of  the  effects  of  buoyancy  on  a  wake  descending  through 
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an  inviscid,  non-turbulent,  incompressible,  stratified  fluid.  It  also 
concludes  that  the  effect  of  buoyancy  is  to  contract  the  vortex  spacing 
and  accelerate  the  wake  motion.  The  analytical  derivation  of  this 
conclusion  is  presented  in  the  next  subsection,  in  which  the  model  is 
described  in  detail.  Extension  of  the  derivation  to  a  compressible 
stratified  medium  is  straightforward,  and  requires  only  that  the  potential 
temperature  or  potential  density  gradient  be  used  in  the  computation  of 
the  Vaisala  period, 
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4.4.3  NEW  THEORY  FOR  TRAILING  VORTICES  IN  A  STABLY 
STRATIFIED  ATMOSPHERE 

Nature  of  the  Theory 

The  new  theory  is  based  on  two  approximations, 
best  illustrated  by  dimensional  reasoning.  We  focus  attention  on  the 
x-y  plane  and  imagine  that  an  aircraft  passes  through  the  origin  of  the 
coordinates  at  time  t  =  0  .  The  aircraft  generates  a  pair  of  vort  es  of 
circulation  _+  T  and  initial  separation  2s^  ,  and  the  vortices  descend 
along  the  y-axis  according  to  some  displacement  history  y(t)  <  0  .  The 
displacement  y (t)  can  be  transformed  into  a  spatial  trajectory  y(z/U)  , 
where  U  is  the  speed  of  the  aircraft  and  z  is  distance  behind,  but  non¬ 
steady  motion  in  the  x-y  plane  is  a  more  convenient  point  of  view  for 
the  present  analysis. 


If  the  atmosphere  were  of  uniform  density,  then 
the  vortex  spacing  would  remain  constant,  and  the  pair  would  descend 
according  to  the  formula 

t 

y  =  -  4ns  0  • 

But  we  are  interested  in  the  case  where  the  ambient  density  pQ  is  a 
function  of  y  ,  say 

dpo 

P0M  =  P0  +  y  if  + . 


where  dpQ/dy  is  negative  in  a  stable  atmosphere,  and  the  first  two 
terms  of  the  Taylor  series  describe  the  density  variation  with  sufficient 
accuracy  near  the  vortices.  Stratification  implies  that  y  may  depend  on 
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pQ  ,  on  dpQ/dy  ,  and  on  the  acceleration  of  gravity  g  ,  as  well  as  on 

the  vortex  parameters  introduced  already.  Thus,  y  =  y  (t,I\  sq.  pQ,  dpQ/dy,  g) 

which  can  be  simplified  by  dimensional  analysis  as  follows: 


y*  (Nt,  NT,  sq/ g  T3  )  . 


N  is  the  Brunt- Vaisala  frequency  for  atmospheric  gravity  waves. 


and  T  is  the  time  required  for  the  vortex  pair  to  induct  itself  down¬ 
ward  a  distance  s  in  the  absence  of  stratification: 
o 


y#  is  the  dimensionless  trajectory  of  the  vortex  pair,  depending  on  the 
dimensionless  time  Nt  ,  and  on  two  dimensionless  parameters  NT  and 
sQ/gT3  .  The  possibility  must  also  be  admitted  that  the  spacing  2s  be¬ 
tween  the  vortices  changes  as  the  pair  descends  into  a  stably  stratified 
atmosphere,  so 

s  =  sos*(Nt,  NT,  so/gTa) 

by  similar  dimensional  reasoning.  Note  that  the  kinematic  viscosity  v 

does  not  appear  in  the  arguments  of  y*  and  s*  .  The  Reynolds  number 

8 

r/v  is  enormous  in  practical  situations,  say  10  ,  and  the  theory  accounts 
for  the  effect  of  stratification  without  reference  to  viscous  or  turbulent 
diffusion. 
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The  parameter  NT  is  a  measure  of  the  strength 
of  the  stratification  relative  to  the  internal  dynamics  of  the  vortex  pair. 
If  NT  is  small,  then  the  vortices  travel  downward  many  times  their  own 
separation  within  a  single  period  of  the  gravity  waves  that  they  generate. 
The  parameter  so/gT®  measures  the  ratio  of  fluid  acceleration  induced 
by  the  vortices  to  the  acceleration  of  gravity.  Local  acceleration  of  a 
fluid  particle  operating  on  the  density  gradient  within  it  generates  vor- 
ticity.  If  so/gTs  is  small,  the  vorticity  arising  in  the  flow  surround¬ 
ing  the  original  vortex  pair  is  mainly  the  result  of  gravity  acting  on 
density  gradients.  To  estimate  NT  and  sQ/g'I^  ,  recall  that 

N  =  /v  &  =  0. 035  sec”  1 
1  c 


for  an  isothermal  atmosphere,  where  y  is  the  ratio  of  specific  heats, 
c  is  the  adiabatic  speed  of  sound,  and  the  numerical  estimate  is  based 
on  y  =  1.4  and  the  sea-level  value  c  =  1100  ft/sec.  An  isothermal  at¬ 
mosphere  is  untypica.lly  stable,  so  0.  035  sec-1  can  be  taken  as  the  prac 
tical  upper  limit  of  N  .  A  Boeing  747  generates  vortices  of  strength 
r  =  9000  ft^/sec  and  spacing  2sq  =  110  ft  during  approach  to  landing 
(Crow  and  Olsen,  1969),  so  the  parameters  entering  the  stratification 
problem  are  of  order 


NT  =  0.  15 


at  most,  and 


So  =  0.10 

gT3 


The  theory  is  designed  to  exploit  the  fact  that  both  parameters  are  much 
smaller  than  unity.  In  fact,  the  acceleration  parameter  SQ/gT3  is  much 


4-52 


less  important  to  the  structure  of  the  theory  than  NT  ,  so  the  trajec¬ 
tory  and  vortex  separation  will  emerge  as  a  one-parameter  family 


y 


r 

4ttNsq 


y*  (Nt,  NT) 


3  =  33*  (Nt,  NT) 
o 


with 


NT  «  1 


The  idea  of  weak  stratification  is  acceptable 
enough,  but  it  leads  to  mathematical  conclusions  that  seem  to  defy  in¬ 
tuition.  The  point  of  view  taken  here  is  that  the  conventional  intuition 
is  appropria'.e,  if  at  all,  to  the  limit  NT  >  >  1  opposite  to  the  limit  Gf 
aeronautical  interest.  Most  previous  theories  imply  that  the  recircu¬ 
lation  cell  surrounding  the  vortices  becomes  buoyant  as  it  descends 
into  heavier  fluid  and  bounces  upward  after  a  time  of  order  N“  1 ,  as 
though  the  cell  were  a  gas-filled  balloon  (cf.  Saffman  1972).  Those 
theories  originate  from  the  work  of  Turner  (I960),  who  was  concerned 
with  vortices  generated  from  rest  by  buoyancy.  But  here  the  wake  is 
assumed  to  descend  too  fast  for  buoyancy  to  generate  significant  local 
vorticity.  The  vortieity  responsible  for  altering  the  trajectory  of  the 
wake  arises  slowly,  primarily  in  a  region  above  the  wake  called  the 
buoyant  upwash.  The  flow  induced  by  the  upwash  draws  the  vortices 
together  and  diminishes  the  volume  of  their  recirculating  fluid.  Ulti¬ 
mately  the  vortices  diffuse  into  one  another,  while  the  upwash  relaxes 
into  Boussinesq  gravity  waves,  which  carry  the  impulse  of  the  aircraft 
wake  to  infinity. 
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Hydrodynamic  Drift 


The  passage  of  a  vortex  pair  leaves  no  resi¬ 
dual  velocity  disturbance  in  an  ideal  fluid.  Fluid  particles,  however, 
suffer  a  net  downward  displacement  during  the  passage,  called  drift 
in  the  hydrodynamic  literature.  Figure  4  - 1 1 .  illustrates  the  drift  of 
fluid  particles  above  the  closed  recirculation  cell  of  the  trailing  vor¬ 
tices.  The  heavy  lines  in  the  drawing  can  be  regarded  as  filaments 
of  ink,  inscribed  horizontally  in  the  x-y  plane  prior  to  tne  start  of 
vortex  motion.  The  filament  BB  below  the  recirculation  cell  as  yet 
is  hardly  disturbed,  but  the^filament  AA  has  been  stretched  around 
the  cell  and  sucked  downward  behind  it.  As  time  passes  the  filament 
approaches  an  asymptotic  displacement  h(x)  ,  which  is  permanent  in 
an  unstratified  fluid.  The  drift  h(x)  was  introduced  by  Darwin  (1953), 
who  showed  that  the  integral 


0 


”h(x)  dx 


-QD 


is  an  apparent  mass,  in  this  case  the  apparent  mass  of  the  fluid  out¬ 
side  the  recirculation  cell.  The  drift  h(x)  diverges  logarithmically 
at  x  -  0  owing  to  the  stagnation  point  of  the  recirculation  cell,  but  the 
integral  is  finite. 


No  analytical  solution  for  the  drift  of  a  vortex 
pair  is  known,  but  h(x)  can  be  computed  numerically  from  the  constru 
tion  of  Figure  4-1?.  For  the  moment  the  coordinates  x-y  are  fixed  to 
the  vortices,  so  the  flow  outside  the  recirculation  cell  is  steady,  with 
an  upward  velocity  V  -  ]f/4-~s  at  infinity.  Fluid  particles  P  and  Q 
move  on  separate  streamlines,  P  passing  close  to  the  recirculation 
cell  ar-.d  Q  remaining  indefinitely  far  away.  Suppose  P  and  Q  come 
from  the  same  stratum  y  far  beierw  the  vortices  and  are  presently 
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separated  a  distance  y  vertically.  Then 


^  =  V  -  v(T,y)  , 

where  v  is  the  vertical  component  of  velocity  at  P  ,  y  is  the  stream 
function  identifying  the  streamline  on  which  P  is  moving,  and  y  is  the 
current  location  of  P  .  Moreover, 

d  A  y  _  d  A  y  dt  _  V  -  v 
dy  dt  dy  v 


so 


(4-7) 


which  gives  the  ultimate  displacement  of  a  particle  on  the  streamline  y 
The  stream  function  for  the  flow  around  a  pair  of  vortices  located  at 
x=_ts,  y  =  0  is  given  by  the  formula 

=  -V  jfc(y  +  (x+af  ]  +*)  ' 

which  can  be  solved  for  x(y,y)  by  iteration.  The  resulting  value  of  x 
can  be  substituted  into  the  formula 

v  =  .  61  =  v[2(x-3)s  _  _  .  2(x+_sJ_s_  +  jl 

t>x  |_y3+(x-s)2  y2  +(x+s)3  J 

to  determine  v(y,  y)  .  The  procedure  is  repeated  for  each  step  of  the  in 
tegration  4-7,  resulting  in  numerical  values  of  the  drift  h(x)  =  h(-¥/V) 

The  drift  is  best  presented  as  a  universal  function  r|(§)  ,  where 
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11 


S(1)T,(S)  ,  § 


(4-3) 


X 

s(t) 


Here  we  take  account  of  the  fact  that  the  separation  2s(t  )  between  the 
vortices  may  change  slowly  with  time,  out  the  shape  of  their  recircu¬ 
lation  cell  and  the  shape  of  the  drift  they  leave  behind  will  be  invari¬ 
ant. 


Figure  4-13  is  a  plot  of  the  dimensionless  drift  p  (?)  •  The 
function  wa  ,  obtained  numerically  by  the  procedure  outlined  above. 
Semi-logarithmic  coordinates  are  used  to  emphasize  the  asymptotic 
behavior  of  p  tor  small  £  ,  ard  in  fact  the  asymptotic  expression 

n=  1.073  2.298^1?  | 


is  rather  accurr .'e  throughout  the  range  0  <  |p|<  1  .  A  quantity  of  major 
importance  in  the  theory  of  wake  descent  is  the  integral 


D 


y  n(?) 


8.  184 


(4-9) 


which  also  was  evaluated  numerically.  It  is  a  curiosity  of  the  theory 
that  the  influence  of  stratification  depends  upon  the  part  of  the  drift  be¬ 
tween  the  vortex  cores,  as  further  analysis  will  show,  rather  than  upon 
the  integral. 


/ 

-CD 


T)(?)d5 


proportional  to  the  apparent  mass  of  the  recirculation  cell. 
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The  Dimensionle 


Buoyant  Upwash 


Drift  would  have  no  dynamical  consequences 
in  a  neutrally  stable  atmosphere.  The  flow  velocity  would  relax  to 
zero  after  the  passage  of  the  wake,  and  fluid  particles  would  remain 
permanently  depressed.  When  the  atmosphere  is  stratified,  however, 
the  drift  results  in  a  density  distribution  that  cannot  persist.  The  den¬ 
sity  distribution  immediately  above  the  wake  is  illustrated  in  Figure 
4 '14,  where  the  darker  bands  represent  the  heavier  strata  of  fluid. 
The  wake  leaves  behind  a  density  field 


<3p0(y) 

p(x,y)  •-  p0(y)  +  h(x)  — —  *  (4-10) 

By  our  sign  convention  h  is  positive,  and  dp^/dy  is  negative  in  a  stably 
stratified  atmosphere,  so  the  density  is  reduced  behind  the  vortex  pair. 


The  next  dynamical  stage  is  that  the  depressed 
light  fluid  oys  upward.  The  initial  buoyant  motion  is  best  treated  in 
terms  of  ti;;  vorticity  equation,  which  can  be  written  as 


VP 

P 


Du 

x  te  -  ~  ) 

~  Dt 


(4-11) 


for  two-dimensional  flow  (vortex  stretching  is  absent).  Here  Q  is  the 
vorticity  aligned  in  the  z-direcfcion,  and  VQ/Dt  is  the  rate  of  change 
of  vorticity  following  a  fluid  particle.  The  right-hand  side  is  the 
Bjerkness  force,  coupling  the  density  gradientVp  to  the  acceleration 
of  gravity  g  diminished  by  the  particle  acceleration  Du/Dt  .  The  ratio 
of  Du  /Dt  to  g  is  of  order  s  / gT3  ,  a  small  parameter  according  to 

v  <v  O 

Section  1,  so  the  particle  acceleration  can  be  neglected  compared  with 
gravity  in  Eqn.  (4-11). 


Suppose  the  current  time  is  t  ,  and  the  trailing 
vortices  passed  level  y(r)  at  some  previous  time  t  •  If  the  time  differ¬ 
ence  (t-r)  is  short  compared  with  the  inverse  of  the  Brunt- Vaisala  fre¬ 
quency  N  ,  then  the  vorticity  Q  is  in  the  initial  stage  of  growth  described 
by  the  linearized  version  of  Eqn.  (4-11)  ,  namely 


81  =.  JL  i£ 

8t  p  dx 
ro 


where  the  cross-product  has  been  performed,  and  pQ  has  been  substi¬ 
tuted  for  p  with  sufficient  accuracy  in  the  denominator.  By  virtue  of 
Eqn.  (4-10). 

BK,  _  _  £_  ffo  dh  n2  dh 

81.  ~  “  p  dy  dx  dx  * 

ro 


and  from  Eqn.  (4-8), 


where  the  prime  denotes  differentiation  with  respect  to  the  dimensionless 
argument  of  p  .  The  vorticity  equation  can  be  integrated  from  the  time 
of  wake  passage  t  to  the  current  time  t  ,  with  the  result  that 

£[x,y(r)j  =  N2(t-T)y[-~-y!  (4-12) 


aa  long  as 


N(t-r)  «  1  . 
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The  ordinate  y  is  given  the  argument  t  as  a  reminder  that  we  are 
dealing  with  the  stratum  y  penetrated  by  the  vortex  pair  at  time  j  , 
when  the  vortex  cores  were  separated  a  distance  2s(t)  . 

The  drift  p(^)  is  positive  and  decreases  sym- 
metrica1ly  around  £  =  0  ,  so  the  vorticity  Q  is  antisymmetric  and  neg¬ 
ative  for  positive  x  .  The  vorticity  distribution  (Eqn.  (4-12)  thus  re¬ 
presents  an  upwelling  of  gathering  strength  behind  the  original  vortex 
pair.  Figure  4-15  depicts  the  vertical  component  of  velocity  with.n 
the  upwelling  flow.  From  the  vantage  point  of  the  vortex  pair,  the 
flow  appears  as  a  concentrated  upwash  whose  strength  increases 
linearly  with  distance  upward,  a  buoyant  upwash.  Figure  4-16 
shows  instantaneous  streamlines  induced  by  the  buoyant  upwash. 

The  induced  flow  tends  to  draw  the  trailing  vortices  upward  and  to¬ 
gether. 


Induction  by  the  Buoyant  Upwash 

The  trailing  vortices  move  under  their  own  in¬ 
duction  and  the  induced  field  of  the  buoyant  upwash.  Suppose  the  right- 
hand  vortex  is  located  at  x  =  s(t)  and  y  =  y(t)  .  Then  it  moves  accord¬ 
ing  to  the  equations 


dt  =  '  4^TT)  +  v  '  (4-  13a) 

ds  _  (4- 13b) 

dt  "  u  ' 

where  r/4rrs(t)  is  the  downward  speed  induced  by  the  left-hand  vortex, 
u  is  the  x-component  of  velocity  induced  by  the  upwash,  and  v  is  the 
y-component.  Both  u  and  v  are  to  be  evaluated  at  x  =  s(t),  y  =  y(t)  . 
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d  of  the  Buoyant  Upwash. 


The  vorticity  distribution  Q  determines  u  and 
v  according  to  the  construction  shown  in  Figure  4-17  An  area  ele¬ 
ment  dA  is  located  at  x  ,  y(T)  a  distance 


r  =  ^[s(t)-x]3  +  [y(-r)  -  y(t)]3 


(4-14) 


away  from  the  trailing  vortex.  A  vorticity  element  QdA  has  evolved 
since  time  j  and  currently  induces  a  velocity 


du  =  £dA 


[y(T )  -  y(t)] 
2rrrs 


(4-15a) 


dv  =  £dA 


[s(t)  -  x] 

Zirr3 


(4- 1 5b) 


The  induced  velocity  u,  v  follows  by  integrating  Eqns.  (4-15)  over 
the  x-y  plane.  It  is  convenient  to  introduce  the  transformation 


dA  =  dx  dy(t) 


dr 


If  N(t-t)  is  everywhere  small,  the  vorticity  field  (4-12)  can  be  combined 
with  Eqns.  (4-15)  to  produce  the  following  results: 


=J  dr  /  dx  N3(t-T)ri'  [■—) 


dt 


V.T)  -  y(t)~] 

/nr3 


djr  [3({)  ~  X1 

dT  2tti* 


where  r  is  given  by  Eqn.  (4-14)  A  partial  integration  over  x  re¬ 
moves  the  derivative  on  q  ,  and  finally 
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U  =  /dr  J°°  dx N2(t-T ) s(r  )r\  [-^ ]  ||X  |  1 


(4- 16a) 


V.  fdr^dxN^t-TWT),^]  ||Z|{to).-y(T)|24-[x-3|.)l2)  (4  l6b) 

°  ~°°  (  2irr  ) 


Now  consider  the  structure  of  Eqns.  (4-13) 
and  (4-16).  The  integrations  over  x  in  Eqns.  (4-16)  could  be 
carried  out  numerically  if  the  current  value  s(t)  and  its  previous 
values  s(t)  were  known.  The  remaining  integrals  over  T  involve  the 
prior  history  of  the  vortex  trajectory  s(T)  ,  y(T)  weighted  by  the  mem¬ 
ory  function  Ns(t-T)  .  Equations  (4-13)  and  (4-16)  are  thus  a  set 
of  integro-differential  equations  for  the  trajectory  s(tl  ,  y(t)  .  They 
could  be  solved  step-by-step  numerically  up  to  a  time  when  Nt  ceases 
to  be  small.  Beyond  that  time  the  memory  function  Na(t-r)  becomes 
suspect  over  earlier  times  T  .  We  should  expect  the  memory  function 
to  become  oscillatory  in  the  manner  of  an  atmospheric  gravity  wave. 

lo  pose  the  correct  gravity-wave  problem, 
we  must  examine  the  integrand  of  Eqn.  (4-  16b)  under  the  simultan¬ 
eous  conditions 


N(t-T)«  1 


(4-  17a) 


Ay  =  >'(t)  -  y(t)  »  s(r) 


(4-1'h) 
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-^Sbi:  ■'  a  .  ''  *'  .. ^  ,  .'■- 


Both  inequalities  can  be  satisfied  if  the  parameter  NT  is  sufficiently 
small,  because 

av  ~  a  . 

Inequality  Eqn.  (4- 17a)  means  that  we  are  dealing  with  a  time  lapse 
(t-r)  such  that  the  simple  vorticity  formula  Eqn.  (4-12)  is  still  valid, 
and  ^qn.  (4- 17b)  means  that 

(ay)2  -  [x-»(t)lz  _  _ 1 

2xr  2r(Ay)2 

lor  any  x  such  that  ti[x/s(t)1  is  appreciable.  Then  the  integrand  of 
the  integral  over  j  is  approximately 

7^2*  [l&i  SW~  ’><«»'*«]  <4' 

2r(Ay)  L  -oc  J 

where  the  term  in  square  brackets  is  the  total  rate  of  volumetric  dis¬ 
placement  per  unit  time.  The  exact  shape  of  the  drift  r(?)  is  imma¬ 
terial  at  large  A  y  ,  which  suggests  that  we  look  for  gravity  waves  gen¬ 
erated  by  a  concentrated  impulsive  displacement  of  fluid. 

Boussinesq  Far  Field 

A  concentrated  impulsive  displacement  can  be 
represented  as  a  product  of  delta  functions; 


The*  meaning  of  Eqn.  (4- 19}  is  that  an  element  of  fluid  of  vanishing 
area  A  is  impelled  downward  an  arbitrarily  large  distance  D  in  such 
a  way  that  the  product  AD  equals  a  finite  volume  Q  .  The  concentrated 
displacement  imparts  an  infinite  potential  energy  to  the  fluid,  which 
serves  as  an  infinite  reservoir  for  the  growth  of  kinetic  energy  later  on. 
The  problem  is  similar  to  the  impulsive  displacement  of  a  free  water 
surface,  as  treated  at  length  by  Lamb  (1932). 

Equation  (4-19)  is  valid  only  near  t  =  0  . 
Throughout  the  course  of  its  evolution,  the  vertical  component  of  vel¬ 
ocity  v  satisfies  the  Boussinesq  wave  equation 

-tvZHN2^  =  -V  V2  [Q6(x)6(y)6(t)]  (4-20) 

8t4 

where  the  inhomogeneous  forcing  function  on  the  right-hand  side  in¬ 
sures  that  Eqn.  (4-19)  is  satisfied  near  t  =  0  .  The  wave  Eqn. 

(4-20)  is  based  on  the  assumption  that  v  is  small  enough  to  render 
nonlinear  convection  unimportant  (cf.  Yih  1969).  We  are  interested 
in  the  response  of  the  stratified  atmosphere  to  drift  imposed  many 
vortex  spacings  away  (inequality  Eqn.  (4-17b),  so  the  lineariza¬ 
tion  implicit  in  Eqn.  (4-20)  should  be  all  right.  Once  Eqn.  (4-20) 
is  solved  in  the  form 


v(x,  y,  t)  =  Qf(x,  y,  t) 

for  quiescent  conditions  prior  to  time  zero,  then  the  solution  for  a  dis¬ 
placement  source  moving  down  the  y-axis  can  be  written  as  a  Duhamel's 
integral: 


(4-21) 
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If  the  theory  is  correctly  formulated,  then  the  integrand  of  Eqn.  (4-21) 
at  x  =  0  should  match  Eqn.  (4-18)  in  the  limit  N(t-r)  «  1*  A  condi¬ 
tion  for  matching  ought  to  be  that 


«  =  |£|  .2(T)/\<eXH  .  (4-22) 

•00 


whereby  dQ/dr  can  be  identified  as  the  rate  of  creation  of  drift.  The 
temporal  behavior  of  f(0,  £  y,t-r)  should  yield  a  uniformly  valid  mem¬ 
ory  function,  matching  the  function  N2  (t— r)  of  Eqn.  (4-18)  for  small 
values  of  N(t-T)  .  and  displaying  oscillations  for  large  N(t-x)  .  The 
approach  resembles  a  singular-perturbation  theory,  with  the  vorticity 
distribution  (4-12)  as  the  inner  solution,  and  the  Boussinesq  wave 
field  as  the  outer  solution. 


The  method  for  solving  dispersive  wave  equa¬ 
tions  like  (4-20)  was  developed  in  full  generality  by  Lighthill  (1965). 
The  solution  is  assumed  to  have  the  form  of  a  double  Fourier  integral, 


oo  oo  i(kx+fy )  p 

v  =  J  dk  J  dl  e 


iwt 


.00  -00 


Fj(k,I)e  +F2 


-iutl 

(k,  1)  e  J  (4-23) 


in  which  u)(k, i)  is  chosen  to  satisfy  the  homogeneous  wave  equations 
for  t>  0  ,  and  Fx  and  Fa  ars  chosen  to  model  the  particular  exci¬ 
tation  that  takes  place  at  t  =  0  .  In  the  case  of  Eqn.  (4-20) 


where 


The  phase  velocity  of  the  waves  and  their  group  velocity  play  an 


important  role  in  evaluating  Eqn.  (4-23)  in  the  limit  of  large  Nts 


w  ,k  ,1  ,  Nk  n  _  .  „  , 

c  =17  e  +  ^  e  )  =  —-5-  (k  e  +  i  e 

~p  K  K  K~y'  ^3  ~x 


and 


%  =  W  ®  x  +  8i  ®  y 


(ie  -  ke  ) 
3  ~x  ~y 


(4-24) 


K 


where  e  and  e  art,  unit  vectors  in  the  x  and  v  directions.  The 
-x  ~y 

phase  velocity  and  group  velocity  are  perpendicular  to  each  other, 


c  *  c  =  0  , 
~p  ~g 


w'hich  means  the  disturbance  radiates  outward  as  a  fan,  with  radii  of 
constant  angle  0  =  tan-i(x/y)  being  lines  of  constant  phase.  Mowbray 
and  Rarity  (1967)  evaluated  Eqn.  (4-23)  in  the  limit  of  large  Nt  by 
the  method  of  stationary  phase  and  presented  some  beautiful  photo¬ 
graphs  of  the  Bou8sinesq  wave  field  around  a  point  disturbance. 


The  method  of  stations ry  phase,  however,  is 
not  applicable  to  the  present  analysis.  In  the  first  place,  we  need  the 
behavior  of  Eqn.  (4-23)  in  the  limit  of  small  Nt  in  order  to  match 
the  inner  solution  (4-18)  and  secondly,  we  must  fine  v  on  the  y- 
axis,  where  the  method  of  stationary  phase  happens  to  break  down. 

The  origin  of  the  breakdown  is  evident  in  Eqn.  (4-24);  waves  pro¬ 
pagating  directly  down  the  y-axis  must  have  I  -  0  ,  in  which  case 

c  =  0  .  We  must  turn  to  the  full  solution  (4-23),  which  takes  the 
~g 

form 

v  .iL  [-H(t)  -8-  /  dk  J  di  £  sin  ^coskx  cosfyj  (4-25) 

8t  L  ?rN  o  0  J 
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for  the  problem  posed  in  Eon.  (4-20).  H(+)  is  the  Heaviside  unit 
step  function,  and  Eqn.  (4-25)  is  valid  throughout  all  time 
-»  <  t  <eo  .  For  t>0  ,  the  time  derivatives  can  be  carried  through 
the  integrations 

v  =  /  dk  J  dj  |r  sia^p  coskx  cosfy  , 

V  o  o 

t  >  0  . 

In  particular 

oo  . 

£(d,  y,t)  =  ""j  /  dk  J  d!  gsia  cosfy  .  (4-26) 

K  O  0 

which  gives  the  reoponae  along  the  y-axis  to  a  unit  volumetric  displace¬ 
ment  at  the  origin. 


Expression  {4*  26)  can  be  evaluated  by  trans¬ 
forming  the  variables  of  integration  from  k,  i  to  K,  i .  Thus 


f(0,  y,  t)  =  /  dK  /  df  sin  V?  -7 

v  o  o  ' 


cos  ty  . 


According  to  entry  3.  876.  7  in  the  tables  of  Gradshteyn  and  Ryzhik  (1965), 


K. 


r~$ — 


J  sin  (P  )cosiy  df  =  -  f  j?  JQ  (K  Vy^+P^)  . 


8p 


Thus 


„T  »  J' 
f(°.y.t)  =  -  -jy  / 

O 


V(Ky)*  +  (Nt)*] 


4 


Nt  K  dK 


Ky)  +  (Nt) 
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Substitute 


=  V( Ky? 


+  (Nt)4 


with  the  result  that 


f(0,  y,t)  =  - 


“4  /  wdK 

2iry*  Nt  0 


N  t  J  (Nt) 

,  2 
2iry 


(4-27) 


Equation  (4-27)  meets  the  objective  of  the  Boussinesq-wave  analysis. 
JQ(Nt)  is  a  Bessel  function  of  the  first  kind,  with  asymptotic  behavior 
such  that 


f(0,  y,t) 


cos(Nt-?r/4) 

y 


as  Nt  “*oo  . 

The  impulse  response  f(o,y,  t)  oscillates  in  time  but  not  in  space,  in 
accord  with  the  disappearance  of  group  velocity  in  the  y-direction. 

The  amplitude  of  the  oscillation  grows  parabolicaliy  with  time,  but  an 
infinite  reservoir  of  potential  energy  is  available  to  feed  the  disturb¬ 
ance  as  in  the  case  of  the  analogous  water-wave  problerri  (Lamb  1932). 


Combining  equations  (4-27)  and  (4-21), 
we  find  the  Duhamel's-integral  representation  for  the  vertical  com¬ 
ponent  of  velocity  beneath  a  moving  source  of  displacement: 


v(0,y,  t) 


t  N2(t-T)Jo[N(t-T)]  4Q  ^ 
o  2»(Ay)2  dT 


(4-28) 


where 


Ay  =  y  -  y(T)  . 


Since 

J  (o)  =  1  , 
o 

the  integrand  of  (4-28)  does  indeed  approach  (4-18)  in  the  limit 
of  small  N(t-r),  provided  dQ/d  is  chosen  to  satisfy  (4-22).  Equa¬ 
tion  (4-16b),  moreover,  preserves  its  character  in  the  limit  (4- 17a) 

and  approaches  (4-28)  in  the  limit  (4- 17b),  if  the  approximate  mem- 

2, 

ory  function  N  (t-r)  is  replaced  with 

N2(t-T)  J0LN(t-T)]  .  (4“2< 

Wake  Trajectory 

The  trailing  vortices  move  according  to  Eqn. 
(4-13),  in  which  u  and  v  represent  the  indication  of  the  buoyant  up- 
wash.  Equations  (4-1 6)  specify  u  and  v  for  small  values  of  N(t-r)  , 
and  the  equations  can  be  generalized  for  arbitrary  times  by  replacing 
N3  (t-r)  with  the  more  general  memory  function  (4-29).  The  result 
is  a  uniformly  valid  set  of  integro-differential  equations  for  the  wake 
trajectory: 


4  -  7  S 


(4-30a) 


(4-30b) 


with  r  still  given  by  Eqn.  (4-14).  The  equations  could  be  solved  on 
a  computer  for  given  values  of  p,  N,  and  the  initial  conditions 


y(o)  =  0  ,  b(o)  =  so  , 


and  the  exercise  might  produce  some  interesting  results.  Implicit  in 
Eqns.  (4-30),  however,  is  the  assumption  that 

6  =  NT«  1  , 

where  T  is  the  inner  time  scale  4jjB^  /p  discussed  on  page  4-51 
otherwise  the  buoyant  upwash  would  commence  before  the  drift  had 
approached  its  asymptotic  value.  So  far  the  full  power  of  the  assump¬ 
tion  5  «  1  has  not  been  exploited. 

To  appreciate  the  simplification  that  takes  place 
for  very  small  values  of  5  ,  we  should  scale  the  variables  appearing  in 
Eqns.  (4-30)  as  follows: 
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where 


r2  =[y(r)-y(t)]2  +  62[|a{T)-a(t)]2  (4-32#) 

and  the  integration  over  x  has  been  transformed  into  an  integration 
over  5  .  An  examination  of  the  order  of  magnitude  of  the  various 
terms  in  Eqns.  (IV-31)  reveals  that 

%  —  m  61  • 
ah°»>  • 

as  §  0  .  The  quantity  ^  fin  6  is  very  much  smaller  than  g  ,  so  the 

buoyant  upwash  moves  the  vortices  horizontally  much  more  effectively 
than  it  retards  their  descent  vertically.  For  the  purpose  of  calculating 
the  trajectory  for  very  small  5  ,  it  is  sufficient  to  retain  the  right- 
hand  side  of  (4-3 l*b)  as  the  only  buoyant  effect. 

If  5  were  set  equal  to  zero  in  the  integrand  of 

(4-31  *b),  then  the  integral  would  diverge  near  (t-r)  =  0.  The  effect 

4 

of  finite  6  is  to  prevent  the  denominator  r  from  going  to  zero  at 
(t-T )  =  0  .  The  functions  s(t)  and  y(r)  influence  the  value  of  the  integral 
only  near  t  =  t  .  With  sufficient  accuracy  for  the  purpose  of  integrating 
in  the  limit  of  small  5  , 


y(t)  -  y (r)  w  (t-r)  ^  , 


r4* 


*$&"-**]  ■ 


and  so  forth.  The  outcome  of  such  approximations  is  that 
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d8  62  3{,v  fcM  ,°°  .r  S . 


(t-T )%({)( 1-8)  _ 

where  the  order  of  integration  has  been  reversed.  The  integration  over 
x  can  be  performed  explicitly.  Set 


°1M> 

dy/dt 


Then 


f*  (t-r)2dT  _  1  _ -1  ,t  , 

j  jy-y-r  =  tan 


Jo  " 26a 


6a'  E(t2+62a2) 


which  approaches 


»  lit 


46|a|  46s  |l-e 


in  the  limit  of  small  $  .  Thus 


(dy/dt)3  M| 
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Since  r)(g)  is  a  symmetric  function,  the  factor  sgn(l~£)  causes  con¬ 
tributions  to  the  integral  from  the  region  £>  1  to  cancel  those  from 
?  <  1  .  Thus 


f  ’d(6)  sgn  (l-£)d£  -  f  *  D 

J  “  -l 


-00 


where  D  is  the  constant  evaluated  in  Eqn.  (4-9).  Evidently  the  buoy¬ 
ant  upwash  outside  the  span  of  the  two  vortices  contributes  nothing  to 
their  rate  of  separation.  The  phenomenon  is  analogous  to  a  spherical 
gravitating  shell,  which  exerts  no  force  on  a  particle  in-  .de. 

The  trajectory  equations  have  reduced  to  a  very 
simple  form  in  the  limit  of  small  5  : 


dy  _  1 

"  a  • 


(4-33*a) 


ds  6D  a 

3F  =  T  "SyTat  * 


(4-33*b) 


to  be  solved  in  conjunction  with  the  initial  condit;  jns 


yto)  =  0,  s(o)  =  1  . 


(4-34*) 


The  solution  is  also  simple.  Take 

d^y  1  ds  dl> 

"  7  ar  ■  Ta^rar 


or 


Integrate  from  0  to  t  ,  with  dy(o)/dt  =  -1  according  to  (4-33*a)  and 
(4-34*): 

$Dt 


or 


d£  _  ...  6Dt  a 

at  --(1  +  T') 


The  vortex  spacing  follows  at  once  from  (4»33*a), 

...  «i  j.  SDt 
s(t)  =  (1  +  -y) 


(4-35*) 


and  after  a  further  integration  with  respect  to  time, 

3/2 

...  4  . ,  ,  6Dt.  . 

y<t)!=“36D  (1  +  T“)  ”  1 


(4-36*) 


The  trailing  vortices  gradually  draw  together  under  the  induction  of 
the  buoyant  upwash,  and  they  travel  downward  under  their  own  induct¬ 
ion  at  an  ever  increasing  speed. 

Concluding  Remarks 


Solutions  (4-35*)  and  (4-36*)  of  the  trajectory  equa¬ 
tions  can  be  written  in  terms  of  dimensional  quantities  as  follows: 


i(t)M  Tt}^  . 

Q  & 


(4-37*) 


4*o  f  D  2  1 

?m  *  - fi  *yN~rt)  -if  (4- 

3i xjstt  1  J 
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where 


and  D-8. 134  . 


Figure  IV-18  is  a  plo.  of  the  vortex  separation  2s(t)  and  displacement 
y(t)  for  the  example  mentioned  on  page  IV- 52,  a  Boeing  74?  during  ap¬ 
proach  to  landing  ia  an  isothermal  atmosphere.  The  effect  of  stratifi¬ 
cation  is  not  very  strong  under  such  conditions.  The  vortex  separation 
is  halved,  for  example,  only  after  140  sec,  by  which  time  the  mutual- 
induction  instability  (Crov.  1970)  would  probably  have  destroyed  the 
vortices  anyway.  The  leveling  of  trailing  vortices  observed  by  Condit 
and  Tracy  (1971)  should  probably  be  attributed  to  instability  rather  than 
buoyancy. 


The  really  interesting  consequence  of  the  solu¬ 
tion  (4-371  is  qualitative.  Intuitively  it  would  seem  that  stable  stratifi¬ 
cation  ought  to  have  somewhat  the  same  effect  as  a  rigid  barrier,  whieb 
would  cause  the  vortices  to  decelerate  and  spread  apart.  The  simplest 
derivation  of  'be  Brunt- Vaisala  frequency,  moreover,  involves  the  idea 
of  a  blob  ot  fluid  bouncing  under  tbs  action  of  buoyancy.  It  is  hard  to 
abandon  the  rigid  barrier  aad  bouncing  blob  as  intuitive  models,  but  that 
is  what  Eqn  (4-37)  ask  ue  to  do.  Throughout  the  development  we  have  as¬ 
sumed  that  $  -  NT  is  small,  m  other  words,  that  stratification  exerts  a 
weak  (but  sustained)  influence  jo  the  local  dynamics  of  the  vortex  pair. 
The  rigid-barrier  analogy  might  work  for  NT  >>  I  ,  and  the  textbook 
bouncing -blob  irgumeri  is  tailored  for  Sousslaesq  waves  having  NT  -  I 
by  definition.  Perhaps  it  should  not  be  too  surprising  that  the  flow  is 
qualitatively  different  whoa  NT  <<i  . 

The  physical  model  implied  by  the  theory  can  He 
appealing,  eeee  cae  accommodates  toe  idea  that  the  vortices  contract 
together. 
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Figure  4-18.  Wake  Trajectory  for  the  Case  T  =  9000  ft.  /sec 
2sq  =110  ft.,  N  =  0.  035  sec”*  .  The  dashed  line  applies 
for  zero  stratification. 


Figure  4-19  illustrates  the  density  field  and  streamlines  near  the  con¬ 
tracting  recirculation  cell.  The  white  area  represents  fluid  captured  at 
the  instant  of  wake  formation,  and  the  shaded  area  represents  the  den¬ 
ser  fluid  at  the  present  wake  location  y{t)  .  The  buoyant  upwash  is  not 
represented  (cf.  Figs.  4-14  through  -16)  but  is  squeezing  the  recircu¬ 
lation  cell,  forcing  light  fluid  to  drain  xipward  near  the  rear  stagnation 
point.  Vorticity  arises  at  the  interface  between  light  and  heavy  fluid, 
bixt  the  trailing  vortices  sweep  the  vorticity  back  toward  the  rear  stag¬ 
nation  point*  The  interfacial  vorticity  becomes  the  boundaries  of  the 
drainage  filament,  as  described  by  Scorer  and  Davenport  (1970).  As 
far  as  the  present  theory  is  concerned,  the  drainage  filament  is  an  in¬ 
terior  detail  of  the  buoyant  upwash  that  eliminates  the  integrable  singu¬ 
larity  T)  (o)  of  the  drift.  The  wake  accelerates  downward  but  loses 

captured  fluid  and  impulse  to  the  buoyant  upwash.  Ultimately  the  im¬ 
pulse  due  to  the  aircraft  radiates  away  as  gravity  waves,  while  the 
trailing  vortex  cores  overlap  and  annihilate  each  other. 

A  preliminary  study  suggests  that  the  interfacial  vorticity  may  con¬ 
tribute  significantly  to  the  induced  velocity  at  the  trailing  vortex  cores. 

If  so,  the  contribution  should  be  added  to  the  induced  field  of  the  buoy¬ 
ant  upwash  treated  in  this  report.  The  quantitative  predictions  of 
Figure  4-18  may  change,  but  the  qualitative  conclusion  that  the  vor¬ 
tices  draw  together  and  accelerate  downward  will  be  reinforced.  The 
effect  of  interfacial  vorticity  will  be  incorporated  in  a  further  report. 
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5.  VORTEX  DECAY 

5.  1  GENERAL 

It  appears  that  there  are  in  general  two  major  types  of  vortex 
decay,  (I)  the  catastropic  instabilities  associated  with  Crow  linking 
and  vortex  breakdown,  and  (2)  the  regular  dissipative  decays.  It  is 
believed  that  in  most  cases  the  instabilities  terminate  the  vortex  life. 
These  are  discussed  in  detail  in  Sections  5.  Z,  5,  3.  For  complete¬ 
ness  we  discuss  here  the  general  dissipative  modes  although  (except 
in  special  circumstances)  they  are  probably  less  important  than  the 
instabilities. 


Calculations  for  the  regular  laminar  dissipation  of  a  single  vis¬ 
cous  vortex  show  that  the  rate  of  decay  is  more  gradual  than  any  de¬ 
cay  observed  in  flight  tests.  This  causes  one  to  expect  that  turbu¬ 
lent  processes  must  come  into  play,  although,  as  pointed  out  by 
Chigier  (1972)  and  Owen  (1970),  there  are  grounds  to  expect  that 
the  core  flow  field  will  tend  to  suppress  turbulence  and  even  to  re- 
laminarize  the  flow.  We  see  that  there  are  three  mechanisms  which 
can  introduce  turbulence  into  the  vortex  core  dynamics;  the  first  is 
tlie  turbulence  of  the  ambient  field  it  is  immersed  in,  the  second  is 
the  turbulence  generated  by  other  elements  of  the  aircraft  --  for  ex¬ 
ample,  the  propulsion  system  --  and  the  third  is  the  mechanically 
generated  turbulence  by  motions  in  the  vortex  cell.  The  aircraft  tur¬ 
bulence  can,  for  these  purposes,  be  regarded  simply  as  increasing 
the  ambient  turbulent  field  in  which  the  vortex  is  immersed,  so  we 
confine  our  attention  to  ambient  and  mechanical  turbulence.  The 
general  effect  of  these  is  presumably  to  increase  dissipation  rates 
resulting  in  increase  in  core  and  cell  size  and  reduction  in  descent 
velocity.  These  effects  have  been  fully  discussed  in  the  section  on 
vortex  transport. 
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We  now  discuss  the  other  factor  which  may  be  expected  to  affect 
the  general  core  dissipation  ~~  this  is  the  presence  of  organized  vor- 
ticity  in  the  external  field.  Such  a  situation  occurs  in  a  uniformly 
sheared  crosswind.  It  would  be  expected  that  the  presence  of  this  am¬ 
bient  vorticlty  (of  one  sign)  would  have  opposite  effects  on  the  two  core 
vorticity  distribution;  causing  the  cores  to  develop  differently,  and  pos¬ 
sibly  causing  the  pair  to  tilt.  This  phenomenon  has  been  fully  discussed 
in  the  section  on  vortex  transport.  It  is  believed  that  this  is  the  effect 
causing  the  appearance  of  the  solitary  vortex  which  is  frequently  excep¬ 
tionally  long  lived. 

1  He  otner  process  by  which  organized  vorticity  is  developed  in 
the  ambient  fluid  is  buoyancy.  Although  the  effects  of  buoyancy  are 
highly  controversial,  it  is  clear  that  the  buoyancy  field  will  generate 
vo’-ticuy  both  in  the  out^r  flow  due  to  non-uniform  vertical  displace¬ 
ments  of  the  stra*ified  density  layers,  and  on  the  cell  boundary  itself 
due  t  i  the  density  discontinuities  between  the  cell  flow  and  the  external 
flow.  A  new  approach  ti  th’s  problem  has  been  described  and  analyzed 
in  Section  4.  4. 


I./:;.;;. 


As  postulated  ir„  the  encounter  hazard  section,  it  is  believed  f:hat, 
with  the  exception  of  the  solitary  vortex,  the  regular  vortex  dissipa¬ 
tion,  by  whatever  process,  is  sufficiently  gradual  that  its  effects  on 
alleviating  the  encounter  hazard  are  not  important.  In  most  cases  the 
vortex  pair  retains  its  danger  potential  until  it  id  destroyed  oy  one  of 
the  two  major  instabilities.  These  -stabilities  are  discussed  in  de¬ 
tail  in  the  following  sections. 


5.2  LINKING  INSTABILITY 

Introduction 

A  theory  for  the  stability  of  a  pair  of  trailing  vortices  has  been 
developed  by  Crow  (1970),  This  theory  examines  the  time  dependent 
displacements  of  vortex  lines  under  an  initial  hypothetical  perturbation 
in  the  wake  of  an  aircraft  Under  the  mutual  influence  of  the  velocity 
fields  associated  with  the  perturbed  vortex  lines,  the  displacements 
become  unstable  and  grow  with  time.  Sinuous  oscillations  develop  on 
the  vortex  lines,  which,  due  to  the  instability,  can  grow  sufficiently 
large  to  cause  the  two  vortex  lines  to  touch  and  subsequently  link  to¬ 
gether.  This  linking  occurs  periodically  in  the  downstream  direction, 
z  ,  at  distances  corresponding  to  the  wavelength  of  the  sinuous  oscill¬ 
ation  of  the  vortex  pair.  At  each  link  point,  each  original  vortex  line 
is  severed,  the  two  originally  parallel  vortices  join  together,  and  a 
closed  vortex  loop  or  ring  is  formed.  Such  a  vortex  ring  is  formed 
between  each  link  point,  resulting  in  a  train  of  vortex  rings  as  the 
wake  of  the  aircraft  continues  to  develop.  The  coherent  flow  in  the 
wake,  characteristic  of  the  energetic  flow  around  the  original  organ¬ 
ized  trailing  vortices,  abruptly  changes  in  character  once  the 
vortex  ring 8  have  formed. 

It  is  shown  by  Crow  that  several  modes  of  instabilities  can  occur, 
although  it  is  argued  that  the  structure  of  atmospheric  turbulence,  acting 
in  the  capacity  of  providing  forced  turbulent  excitation  of  the  instability, 
imposes  an  overwhelming  bias  in  favor  of  long  waves,  which  are  unstable 
only  in  the  symmetric  mode.  Actual  observations  have  confirmed  this 
long  wave  interaction  behavior,  where,  from  the  Crow  theory,  the  wave¬ 
length  is  appi oximately  8.  6  times  the  original,  undisturbed  vortex  spacing 
(Smith  and  Beesmer,  1959;  Chevalier,  1973). 


Vortex  lines  move  in  accordance  with  the  velocity  field  of  the  surround-  . 
ing  fluid.  Atmospheric  turbulence  provides  a  continuous  disturbance  of 
the  vortex  lines  and  convective  turbulent  excitation  of  the  resultant  in¬ 
stability.  Once  initially  distorted,  the  vortex  lines  also  move  in  accord¬ 
ance  to  their  own  mutual  and  self-induced  velocity  fields,  which  was  the 
condition  originally  treated  by  Crow  in  the  analysis  of  the  instability. 

No  attempt  was  made  to  account  in  detail  for  the  effects  of  atmospheric 
turbulent  excitation,  except  to  state  that  such  turbulence  is  the  probable 
factor  which  excites  the  long  wave  instability. 


The  present  analysis  examines  the  effects  of  atmospheric  turbu¬ 
lence  as  an  input  forcing  function  for  the  instability.  The  time  for  the 
instability  to  grow  to  the  point  at  which  vortex  linking  occurs  is  found 
as  a  function  of  the  vortex  strength,  po,  and  spacing,  b  ,  and  the  at¬ 
mospheric  turbulence  dissipation  rate,  s  . 


Structure  of  the  Atmospheric  Turbulence 

The  vortices  from  the  generating  aircraft  are  originally  undisturbed, 
but  under  the  influence  of  atmospheric  turbulence  they  are  continuously 
deformed.  The  instability  is  thus  initiated  by  turbulent  convection,  and 
further,  the  turbulence  acts  continuously  as  a  source  of  input  energy  to 
drive  the  instability. 

The  structure  of  the  atmospheric  turbulence  is  assumed  to  be 
stationary  with  respect  to  the  coordinates  of  the  downward  moving  vor¬ 
tex  pair.  That  is,  the  turbulent  eddies  are  assumed  large  enough  that 
any  time  variation  of  turbulence  is  small  as  the  vortices  move  down¬ 
ward  a  distance  i  =  vT^  where  v  -  l”/2nb  ,  and  T^  equals  the  time 
to  linking.  Typically,  v  =  .  5-3m/sec  ,  and  T,  is  of  the  order  50  sec 
for  light  to  moderate  turbulence.  Thin  requires  a  turbulent  "cell"  {in 
which  the  time  variation  of  turbulence  is  small)  of  dimensions  on  the 


The  resulting  behavior  of  the  vortex  motion  under  the  influence 

of  the  turbulent  field  can  be  analyzed  as  a  time- invariant  linear  filter 

interacting  with  the  turbulent  energy  input  spectrum  integrated  over 

all  wave  numbers,  with  e  constant.  The  frequency  response  function 

of  the  filter  is  obtained  from  the  dynamics  of  the  vortex  instability, 

with  a  peak  at  k  ,  the  wave  number  associated  with  the  maximum 
o 

amplification  rate  of  the  instability. 


From  Crow,  k  =  ft  / b  ,  where  b  is  the  undisturbed  vortex- 
o  miax 

pair  spacing,  and  |3  is  the  dimensionless  wave  number.  gmax  is  of  the 
order  one  so  that  .02  s  k  s  .  10  radians /meter,  which  is  within  the 
inertial  subrange  commonly  encountered  in  geophysical  flows 
(MacCready,  1962).  Although  turbulence  supplies  energy  to  the  insta¬ 
bility  at  all  wave  numbers,  the  filter  action  of  the  dynamics  of  the  vor¬ 
tex  instability  responds  only  to  the  energy  contained  in  the  spectrum 
at  or  near  kQ  .  Wave  numbers  not  in  the  inertial  subrange  will  be 
sufficiently  far  to  either  side  of  kQ  that  the  effects  of  the  filter  function 
will  have  attenuated  their  contribution.  Their  exact  analytical  form  is 
thus  not  critical,  and  the  assumption  of  an  energy  spectrum  described 
by  that  of  the  inertial  subrange  will  be  adequate  over  all  wave  space. 


With  these  assumptions,  the  input  turbulent  forcing  function  for 
the  instability  is  taken  as  the  one-dimensional,  transverse  energy  spec¬ 
trum  in  the  Komolgorov  inertial  subrange.  This  energy  spectrum  is 
given  (page  273,  Tennekes  and  Lumley,  1972)  as: 
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with  <x  =  1.5.  The  actual  turbulent  forcing  function  that  drives  the 
instability,  F(k)  ,  will  be  some  percentage  of  F22(k)  ,  given  by: 


F(k)  a  yF22(k) 


(5-2) 


y  will  be  determined  by  the  geometrical  relationship  which 
couples  F22(k)  with  the  dynamics  of  the  instability. 

Geometry  of  Vortex  Interaction  with  the  Turbulent  Field 


In  order  to  determine  the  ratio,  y  ,  between  the  total  available 
one -dimensional,  transverse  energy,  F22(k)  ,  and  that  which  actually 
couples  with  the  instability,  F(k)  ,  the  relationship  of  the  available  tur¬ 
bulent  energy  to  that  contained  in  the  oscillating  vortices  must  be  estab¬ 
lished  with  respect  to  the  interaction  between  the  geometries  of  these 
two  flows. 

With  reference  to  Figure  5-1,  the  instability  for  the  symmetrical 

mode  develops  at  each  vortex  in  a  plane  inclined  at  Q  ■=  48°  to  the  horl- 

8 

zontal.  At  any  downstream  position,  z  ,  the  perturbation  velocities  de¬ 
velop  in  the  x-y  plane  with  velocity  components  u  ,  and  v  respectively. 


Let: 


y,  z,t)  sU^x,  y,  z,  0}  =  ( z,  0)  +  e  v^z,  0} 

U2(x,  y,  z,  t)  a  U2(x,  y,  s,  0)  =  e^u^z,  0)  *  eyv2(z,  0) 


(5-3) 

(5-4) 


be  the  turbulent  velocities  at  the  positions  of  vortices  1  and  2,  respect¬ 
ively.  The  velocities  at  these  two  positions  have  the  same  direction  in 
the  x-y  plane,  so  that  Uj/v^  =  Ug/v^  *  alt^ou8*1  they  have  a  differ¬ 
ent  magnitude. 


If  a  velocity  in  wave  space  is  defined  as: 


W  ^  ii- - 

/  bn(k,  x,  v,  0)  e1*2  dk  =  Un(x.  y,  *.  0)  (5-5) 


(n  equals  1  and  2),  then,  from  Eqns.  (5-3  and  (5-4)  follows  the  defin¬ 
ition  of  ^n(k,  o)  and  &  (k,  o)  in  a  like  manner. 


Since  th.»  ratios  v  /v  and  ‘i  /u  for  any  constant,  arbitrary  wave 
number:  k  ,  and  at  the  same  distance,  z  ,  remain  constant  for  n  =  1 
and  n  =  2  ,  it  is  possible  to  use  these  transformed  velocities  as  forcing 
functions  to  obtain  the  relative  geometrical  response  of  the  two  vortices 
in  the  x  and  y-directions.  Differentiating  Eqns,  (9)  of  Crow  for  the 
symmetrical  instability,  and  substituting  Eqns.  (8),  with  the  trans¬ 
formed  velocities  added  as  forcing  functions  ,  the  relationships  between 
the  components  of  velocity  resulting  from  the  turbulence,  and  the  vel¬ 
ocity  components  associated  with  the  symmetrical  instability  mode  are 
established: 

8&  r 

ir =  77% (i-  * +  *  <VV  {5“6) 
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r 

i  j  + 

2*b 


v-p2u>)xs  +  (Vj  +  v2) 


(5-7) 


where  the  dimensionless  interaction  functions  Y(0),x(£),  and  a)(g)  are 
defined  by  Crow. 


The  instability  develops  in  a  plane  inclined  at  an  angle  8  to  the 
horirontal.and  the  vertical  and  horizontal  velocities  of  the  developing 
instability,  v#  and  ug  ,  are  coupled  by  vg/ug  -  tan  .  Therefore, 
the  equation  of  motion  of  the  system  is  a  linear  combination  of  Eqns. 
(5-6  and  (5-7)  differentiated  with  respect  to  time.  Because  of  the 
linear  relationship,  it  is  sufficient  to  deal  with  only  i  single  component 


S- 


a 


O 


o£  the  equation  of  motion  for  tho  purposes  of  establishing  geometrical 
relationships  between  velocity  components. 

Differentiating  Eqn.  (5-6)  and  substituting  Eqn.  (5-7)  results  in 
the  x- component  of  the  equation  of  motion: 


8fc 


r  2 

(— Sr)  <l-4»+e2u){l+x-P2u>)*0 

2»b  8 


with: 

and 


+  — %  (1  -  4*  +  P2w)(v.  +  v,)  , 

Zto  1  2 

A  8<c  !k,0) 

*s(k’0)  =  0*  -te— ' 


2  A  M 

a  (k)  a  ( — -Sj)  (1  -  iii  +  0  u)  (1  +  x  “  02(j)  • 
2arb 


(5-8) 


(5-9) 


A  solution  exists  of  the  form: 


£  (M>  *  A  *  B  ea(k*  ♦  C  <ra*k)t 
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(S-iO) 


Substitution  of  Eqn.  (5-iO)  in  Eqn.  (5-8)  with  the  initial  condi¬ 
tions  (5-9),  and  from  Crow,  Eqn.  (12),  the  expression  for  3  : 

3 


tan  0 


*  1  -  £  * 


result*  in  the  determination  of  the  coefficients  A,  B,  and  C. 


*  i  tape. 


B  = 


"sm" 


and 


C  = 


*V°2> 


, A  A  %  {°t+02} 

<W*-tsir- 
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Za(k) 


At  large  times  Eqn.  (5-10)  becomes: 


*fl(k.t)  —  B  e^1*1*  =  ~|gpea(k?t 


from  which 


B  b 


«s<k,  t) 
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If  new  the  radial  displacement  of  one  ef  iiie  vortices  from  is* 
original,  undisturbed  position  i*  given  by  * .  thru  by  defLiUioa- 
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Similarly  to  Es|eu  (5-11),  at  large  times: 
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from  which 


.a  A  » 
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(Uj-U,)*  ti£j- 


4  sin  6 
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t*a  0  =  Un  48  =  i.  11  ,  so  that; 
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4  sin  d , 


(5-14) 


If  long  waves  are  considered,  that  is,  L>  > y  »  b  .  tlw  turbulent 
velocity  field  within  the  "cell”  of  dimension  L  would  result  in  velocity 
component*  at  each  cortex  position  such  that 

A  A  A 

V ui  <<:  VV2 


and,  from  Eqrw  (5*14) 
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Thu*.  a  coupling  feetweee  the  vertical  vertex  motion  and  the  at- 
maapherie  turbulence  will  produce  maximum  energy  Input  to  the  de¬ 
veloping  instability.  The  ratio  of  the  energy  >R  the  oscillating  vortex 
io  that  available  from  the  turbulent  field  i*; 
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*rd.  from  £»p*.  (5*1)  tad  Eqa.  ($-*): 
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Dynamics  of  the  Vortex  Motion 

The  radial  displacement  o {  a  vortex  line  from  ita  original,  un¬ 
disturbed  position  is  given  by  r(z,t)  .  Let: 

5»  . 

r(s,  t)  «  |  r(k,  t)  e  dk  . 

-as 


An  exponential  solution  in  time  for  r(k,  t)  ,  similar  to  that  for 
Eqn.  (5-101  exists: 


rfk,  t)  =  A' 


♦  B’ 


e  +■  c 


e-a(k}t 


which  at  large  times  becomes: 


rik,  t) 


afk)t 

B'  e^  '  a 


ea<kH 

~*&r 


Vk)  . 


(5-17) 


(5-18) 


The  mean  square  displacement  of  a  vortex  element  is  defined  by: 


<r£(E,t})  c  (r2) 


ss  «  ((k^k’Js 

o  1  /  {£fk,  t)  ?|k\tj>e  dkdk!  . 


(5-19) 
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(xfM)  tfk'.t)) 


(5-20) 
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Define  a  function  R(k,  t): 


<r2>  =  J  J  R(k, t)6(k+k' Je1 '^+k' ^zdkdk'  (5-21) 

-oo  -oo 


which  has  the  property  that: 

2  00 

(r  )  =  J  R(k,t)dk  ,  at  k'  =  k  .  (5-21a) 

-u0 


Thus,  from  Eqns.  (r  :<))  and  (5-2 J ) 

R(k,  t)6(k+k')  =  (r(k,t)  r(k',t))  . 

Also  define; 

F(k)6(k+k')  =<Vr(k)$r(k')}  . 
Combining  Eqns.  (5-20),  (5-22),  and  (5-23)  gives: 


(5-22) 


(5-23) 


(5-24) 


-  .  ‘A  f  >T 


The  wavelength,  \  ,  is  related  to  the  wave  number,  k  ,  by; 


X  »  -- 


Therefore, 


a(k)  =  a(-k)  , 


r  ea(k)fcY 


F(k)  . 


(5-25) 


Substituting  Eqn.  (5-25)  in  Eqn.  (5-2 la)  gives: 


,  *  [a(k)tf 

!l  >  ■  L  nrarj F(k)  dk ' 


(5-26) 


This  integral  is  significant  only  near  the  wave  number  of  maxi¬ 
mum  amplification  rate,  kQ  ,  for  the  symmetrical  mode  of  instability. 
The  term  in  brackets  represents  the  filter  function  of  the  forced  oscill¬ 
ation,  and  F(k)  represents  the  turbulent  forcing  function. 


Equation  (5-26)  becomes,  upon  squaring  the  term  irj  brackets: 


<r2>  f“  FfSl  e2“<k>‘  dk 

-oo  a  (k) 


(5-27) 


Expanding  a(k)  in  the  exponential  term  in  a  Taylor  series  and  substi¬ 
tuting  in  Eqn,  (5-27), 


a(k)  =  a(kQ)  +  (k-ko)a'(k0)  + 


(k-k  > 


a"(ko) +  , 


a'(kQ)sO  , 


so  that: 


,  2F(k0)  2ta(ko)  «  -tl.,'(k0)l(k-k0)2 

(I  >  =  -j—  •  /  •  dk 

a  (kQ)  o 


(5-28) 


a"(k  )  is  negative  so  that  t  a"(k  )(k-k  )2  =  -t|a"(k  }|(k~k  )2 
0  o  o  |  o  j  o 

Equation  (5-28)  has  a  solution  of  the  form: 


,  2,  2F<ko>  2ta<ko>  r~t r - 

\T  )  ~  ~7 - e  i J - 

aZ(ko)  ftia"(ko)l 


(5-29 


From  the  definition  from  Crow: 


a(k)  =  — y  a(|3)  ,  |3  =  kb  , 


and 


*(k0)  =  — 

27rtj 


G<«W  = 


2rrb 

r 


T  U. 

2  max 


a"(k  )  =  a"(0  )  =  —  a" 

o  2v  'Kmax'  2tt 


max 


The  "lifetime"  of  a  vortex  pair  is  that  time  T.  at  which  <  r 


(b/2  sin  0)  .  Substituting  this  condition  in  Eqn.  (5-29),  and  recalling  that 


F(k  )  =  F(ig2£)  =  e2/3  {.^£)'5/3 

55  sin  6 


?/,  h  5/3  2Tiia 
4.86a  e2/--(~L_)  e  2x?max 


Thus  the  condition  of  time  to  linking  for  the  sinuous  instability  of 

a  vortex  pair  has  been  stated  in  terms  of  the  vortex  parameters  T0  and 

b  ,  and  the  turbulent  dissipation  rate,  e.  The  relationship  between 

these  parameters  involves  the  dimensionless  amplification  rate 

and  wave  number  f3max  •  From  Crow,  for  the  symmetrical  mode, 

B  =  0.  73  and  a(g  )  =  0,  83  .  A  parabolic  approximation  of  q  vs  p 
rnsix  inEX 

gives  a"((3  „  )  =  -3.  12  .  With  a  =  1.5,  Eqn.  (5-31)  becomes: 
m«*x 


ti  =  0,  00271 


(5-32) 


from  which  a  "universal  lifetime"  plot  may  be  constructed  in  terms  of 

the  dimensionless  parameters  r\  and  t  •  Such  a  plot  is  shown  in  Fig. 

(5-2).  Figure  (5-3)  shows  in  dimensional  form  the  time  to  vortex 

1/3 

linking,  T^,  as  a  function  of  turbulence  dissipation  rate,  e  for 
various  aircraft. 
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Correlation  of  Present  Theory  with  Experiment 


The  predictions  of  time  to  linking  for  various  turbulence  levels 
follow  straight  from  the  theory,  and  involve  no  fitted  constants. 

There  is  scant  experimental  information  with  which  to  compare  the 
present  theory.  However,  a  series  of  experiments  in  which  motion 
pictures  were  taken  of  the  trailing  vortices  from  a  light  plane  were 
performed  {Tombach,  1972).  The  linking  phenomenon  was  observed 
in  some  of  these  experiments,  and  the  time  to  linking  was  measured 
from  the  motion  pictures  for  various  le\  fls  of  atmospheric  turbulence. 
The  results  of  these  experiments  are  presented  in  Figure  (5-4)  along 
with  the  present  theory  in  which  the  vortex  parameters,  ro  and  b  , 
characteristic  of  the  light  plane  used  (Cessna  17 0)  were  included. 

The  agreement  between  experiment  and  the  present  tl  ory  is 
quite  good.  The  scatter  in  the  data  are  well  within  reason,  consider¬ 
ing  that  the  experiments  were  performed  under  field  conditions,  with 
random  effects  undoubtedly  influencing  the  limited  amount  of  data  ob¬ 
tained. 

When  future  experimental  data  becomes  available  with  respect  to 
the  vortex  linking  phenomena  for  different  classes  of  aircraft  subjected 
to  varying  degrees  of  atmospheric  turbulence,  the  constants  in  Eqn. 
(5-32)  can  be  modified  (if  required)  in  a  rational  and  consistent  man¬ 
ner  to  account  for  such  effects  as  non-elliptical  wing  loading,  flap  de¬ 
flections,  wing /non-wing  mounted  engines,  etc. 

The  present  theory  was  developed  for  assumed  "large"  times. 

For  this  reason,  the  behavior  of  the  predictions  of  the  theory  for  small 
times,  say  less  than  10-20  sec,  is  not  as  reliable  as  for  larger  times. 
Modifications  to  the  theory  can  be  performed  that  will  extend  the  range 
of  the  theoretical  predictions  into  the  small  time  region.  However,  ex¬ 
amination  of  Fig.  5-3  shows  that  for  the  class  of  large  airplanes  that 
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(Eqn.  5- 32V 


would  present  a  wake  hazard,  this  small  time  region  occurs  in  the  area 
of  atmospheric  turbulence  that  is  described  as  moderate-to-severe.  It 
is  expected  that  for  such  levels  of  turbulence,  other  vortex  transport 
and  decay  mechanisms  may  dominate  the  wake  behavior.  In  any  event 
extreme  accuracy  in  predicting  wake  behavior  at  the  short  lifetime  end 
is  not  required  since,  from  an  operational  viewpoint,  practical  factors 
other  than  wake  turbulence  require  aircraft  spacings  of  greater  than  10 
to  20  seconds. 

Range  of  Validity  of  Theory 

Three  assumptions  were  made  in  the  development  of  the  present 
theory  which  bear  directly  on  the  validity  of  the  results.  These  assump¬ 
tions  were  that  the  disturbance  energy  spectrum  could  be  represented 
analytically  as  the  two  dimensional,  transverse  turbulent  energy  spectrum 
in  the  inertial  subrange,  that  the  turbulence  is  isotropic,  and  that  vortex 
linking  is  a  'long  time"  phenomenon. 

The  wave  number  of  maximum  amplification  rate,  ,  is  defined  by 

K  =  0  /b  where  9  the  dimensionless  wave  number,  is  of  the  order 
o  max  max 

one,  and  b  is  the  vortex  spacing.  For  the  clrss  of  airplanes  that  would 

present  a  wake  hazard,  .  02  £  Kq  £  ..  10  radians/meter.  This  range  in 

wave  number  is  within  the  atmospheric  inertial  subrange.  Turbulent 

excitation  of  the  instability  at  wave  numbers  cutside  of  this  range  would  be 

greatly  attenuated  iue  to  the  filtering  effect  of  the  vortex  dynamics.  Thus, 

the  inertial  subrange  model  of  the  turbulent  energy  spectrum  over  all 

wave  numbers  is  valid  since  only  those  wave  numbers  at  or  near  K  (which 

o 

are  actually  within  the  inertial  subrange)  contribute  significantly  to  the 
growth  of  the  instability,  and  the  model  is  accurate  here. 

The  value  of  the  turbulent  dissipation  rate,  e  ,  is  only  a  scaling 
parameter  and  does  not  affect  ..he  analysis.  The  value  of  e  enters  into 
the  end  result  and  in  this  respect  it  must  be  chosen  properly.  The  figures 
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shown  assume  that  $  has  been  chosen  for  the  case  of  isotropic  turbulence. 
This  is  realistic  for  situations  whore  the  aircraft  wake  exists  in  an  un¬ 
bounded  medium  (at  altitude).  Near  the  ground,  the  turbulence  is  not 
isotropic,  and  an  effective  e  should  be  used.  The  section  on  Meteorology 
gives  guidelines  for  appropriate  corrections  for  e  near  the  ground.  How¬ 
ever,  near  the  ground,  the  vortices  interact  with  the  ground  plane  itself. 
This  interaction  produces  velocity  fields  at  the  vortices  which  alters  the 
nature  of  the  dynamics  assumed  in  the  present  analysis.  Thus,  near  the 
ground  the  theory  may  not  be  directly  applicable  or  at  best  the  constants 
may  differ  from  those  quoted. 

However,  for  the  isotropic  turbulence  situation,  the  theory  agrees 
very  well  with  the  limited  amount  of  experimental  data  available.  For  the 
case  of  e  -*  0,  the  theory  predicts  an  infinite  time  to  linking.  The  figures 
of  link  time  vs.  e  are  drawn  with  logarithmic  axes  so  that  this  limit  is 
not  immediately  obvious  from  the  graphs. 

The  theory  was  developed  for  assumed  large  times  in  which  the  first 
order  exponential  term  in  the  expression  for  the  radial  displacement  at 
the  vortex  lines  dominates.  Exactly  what  constitutes  "long  times"  depends 
upon  the  type  of  airplane  which  is  generating  the  wake.  This  is  so  be¬ 
cause  the  amplification  rate,  a(k),  contained  in  the  exponential  term  de¬ 
pends  upon  the  aircraft  parameters  and  b  .  As  an  example,  for  a  727 
aircraft,  it  is  estimated  that  the  error  in  the  theory  is  lees  than  10%  for 
a  time  greater  than  about  30  seconds.  Thus,  for  this  particular  aircraft, 
the  theory  does  not  give  realistic  results  for  times  less  than  30  seconds. 

The  theory  can  and  will  be  modified  to  account  for  the  "small  time" 
region.  How-ever.  with  respect  to  the  eneouciet*  ha  sard ,  the  small  time 
region  is  of  little  practical  significance  because  factors  other  than  wake 
turbulence  dictate  aircraft  spacings  much  greater  than  would  exist  at 
"small  time*.  " 
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5. 3  CORE  BURSTING 


General 

Gae  major  form  of  vortex  instability  is  core  bursting  or  vortex 
breakdown.  This  mode  consists  of  a  sudden  abrupt  widening  of  the 
vortex  core  and  in  the  case  of  a  smoke -marked  aircraft  vortex,  the 
disappearance  of  the  tracer  elements.  This  phenomenon  has  been 
observed  in  the  laboratory  (Sarpkaya  1970)  and  in  flight  tests  (Tombach 
1972). 


The  Qualitative  details  of  the  breakdown  are  still  quite  obscure. 

In  the  experiments  of  Sarpkaya,  the  first  effect  seemed  to  be  an 
axi  symmetric  disturbance  under  which  the  core  expanded  and  contra '■ted 
smoothly.  Downstream  of  this  was  a  distinct  spiral  disturbance  which 
was  then  followed  by  a  disorganised  but  roughly  axieymmef  ric  core 
widening,  the  final  breakdown.  In  Tombach *s  experiment#  the  same 
sequence  of  events  seems  to  take  place,  although  the  initial  smooth 
axi symmetric  bulging  is  not  as  distinct. 

Several  explanations  of  vortex  breakdown  have  been  proposed 
and  are  Hated  by  Bilanln  and  Widnall  (1973).  These  are  the  stagna¬ 
tion  of  the  axial  velocity  1966;  Bossei,  1969),  the  conjugate 

jump  theory  {Benjamin  1962),  and  a  stability  approach  (Ludwig,  1962). 
None  of  thee*  is  entirely  satisfactory,  yet  all  contain  common  elements 
and  all  se*m  partially  supported  by  experiment.  It  appears  agreed 
Shat  the  rapsd  enlargement  of  the  vortex  care  is  accompanied  by  axial 
pressure  gradients,  a  ad  that  this  breakdown  can  only  occur  when  the 
flaw  ai?p.  jaches  s  certain  critical  combination  of  axial  and  tangential 
3U;w  profiles.  Apparently  the  magnitude  of  the  triggering  adverse 
pressure  gradient  required  depends  upon  the  proximity  of  the  flow  to 
its  critical  state,  it  is  generally  agreed  that  dissipation  and  core 

development  will  cause  the  core  to  approach  the  critical  state  so  that 
ia  most  cases  the  core  will  eventually  develop  to  a  state  capable  of 
breakdown. 
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It  appears  likely  that  this  critical  state  caa  be  described  crudely 
as  a  function  of  the  axial  velocity  on  the  center  line  and  the  swirl  ratio. 
This  swirl  ratio  is  the  ratio  of  the  maximum  tangential  velocity  to  the 
freeetream  speed.  A  recent  paper  by  Mager  (1972)  gives  the  critical 
conditions  in  fair  agreement  with  Benjamin's  results.  Qualitatively, 
both  of  these  papers  claim  that  as  the  axial  velocity  reduces  so  docs 
the  critical  swirl  ratio. 

If  the  critical  axial /tangential  flow  combination  could  be  deter¬ 
mined,  then  the  breakdown  prediction  problem  would  reduce  to  test¬ 
ing  for  the  critical  state  after  computing  the  core  development. 

It  appears  that  the  core  bursting-phenomenon  is  not  related  to 
interaction  between  the  left  and  right  elements  of  the  horseshoe  vortex 
pair,  but  rather  to  the  development  of  the  core  itself  and  this  development 
is  a  function  of  the  kinematic  viscosity  controlling  the  core  development. 
Observations  by  Tombach  ( 1973)  have  shown  that,  at  the  scale  of  his  flight 
test,  in  which  vortices  were  generated  by  a  single  engine  light  airplane 
(Cessna  170)  core  bursting  occurred  first  in  light  atmospheric  turbulence, 
while  in  high  ambient  turbulence  the  Crow  Instability  invariably  terminated 
the  vortex  life. 

Thus  ; t  is  possible  that,  at  the  ~*me  turbulence  level,  the  time 
scales  for  core  bursting  and  Crow  Instability  are  related  to  some  function 
of  the  airplane  Reynolds  Number,  v-feich  ecTd  be  defined  as  Ub/V.  wh#re 
U.  b  are  (light  speed  and  span  and  v  the  kinematic  viscosity.  On  this 
basis  the  vortices  due  to  small  Slow  aircraft  might  be  expected  to  sSisetpaie 
principally  due  to  core  bursting  while  those  due  to  large  fast  aircraft  due 
to  Crow  Instability.  Some  supportive  e^een  s  for  this  sr  that  moss  small 
scale  lab  type  teste  exhibit  core  bursting  (although  Crow  'ns  tab1  '"ity  can 
certainly  be  excited),  while  very  large  scale  flight  test*  usually  ahow- 
sinuou*  instability. 


Thus  it  is  possible  (although  not  definitely  substantiated)  that  Crow 
Instability  is  the  most  significant  mode  of  decay  for  vortices  characteristic 
of  large  transport  aircraft.  However  tinder  unusual,  or  artifically  per¬ 
turbed  circumstances,  core  bursting  may  be  important. 

The  development  of  the  core  is  principally  a  function  of  time, 
vortex  strength,  ambient  turbulence,  viscosity  and  wing  drag.  Many 
attempts  have  been  made  to  analytically  determine  axial  and  circum¬ 
ferential  velocities  in  the  core.  The  significance  of  such  a  calculation 
is  that  if  core  bursting  is  related  only  to  local  details  of  core  velocity; 
then  breakdown  could  be  predicted  if  the  criteria  on  velocity  profiles 
were  known  and  the  profiles  could  be  analytically  determined.  Conse¬ 
quently  we  discuss  some  of  the  global  aspects  of  the  core  development 
problem,  and  specifically  state  the  invariants  involved. 

Force  and  Moment  In  Global  Terms 

Considering  the  steady  flow  behind  a  lifting  ving,  it  is  clear  that 
in  any  plane  normal  to  the  flight  path  that  the  total  drag  r  nd  lift  must 
be  conserved.  For  this  analysis  we  will  ignore  the  propulsive  system, 
noting  that  usually  this  does  not  contribute  to  the  lift  and  that  its 
axial  contribution  will  be  a  thrust  associated  with  the  slipstream  due 
to  the  jet  engines  or  propellors.  Thus  we  will  find  that  considering 
the  wing  only,  application  of  the  momentum  theorem  in  an  axial 
direction  gives  a  drag  on  the  wing  associated  with  its  profile  and  in¬ 
duced  drag. 

If  we  consider  a  large  circuit,  along  the  centerline  and  around 
one  side  of  the  flow  system  (Fig.  5-5)  we  see  that  on  the  outer  bound¬ 
aries  there  can  be  no  viscous  torque  since  the  flow  is  that  of  a  poten¬ 
tial  dipole  having  velocity  gradients  decaying  like  the  inverse  cube  of 
the  distance.  On  the  center  line,  even  if  vorticity  is  present,  the  shear 
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is  zero  by  symmetry,  so  there  can  be  no  force  due  to  laminar  viscosity. 
Thus  no  torque  can  be  applied  to  the  circuit  C  ,  so  the  impulse  moment 
(loosely  called  the  "moment  of  momentum")  is  conserved,  implying  that 
the  angular  momentum  on  each  side  is  conserved.  This  observation 
seems  first  to  have  been  made  by  Betz  (1933)  by  considering  a  set  of  free 
trailing  vortices.  We  digress  temporarily  to  show  that  the  polar  moment 
of  vorticity  referred  to  by  Betz  is  directly  related  to  the  impulse  moment. 

Relationship  Between  Impulse  Moment  and  Polar 

Vorticity  Moment 

It  is  shown  in  Lamb's  Hydrodynamics  that  all  potential  motions  can 
be  described  by  an  impulsive  wrench  applied  to  the  fluid.  This  wrench 
consists  of  both  force  and  couple  constituents.  The  force  term  is  most 
familiar  and  is  usually  called  the  impulse.  We  note  that  the  impulse  is 
not  strictly  equal  to  the  momentum  of  the  system,  which  in  many  cases 
is  mathematically  indeterminate.  However,  the  impulse  is  frequently 
loosely  referred  to  as  the  "momentum.  "  In  our  analysis  we 
refer  to  the  moment  of  the  impulse  as  the  "moment  of  momentum,  " 
recognizing  again  that  this  term  is  mathematically  indeterminate,  but 
that  the  impulse  moment  is  a  well-defined  quantity.  Here  we  show  how 
this  is  connected  to  the  polar  vorticity  moment  developed  by  Betz 
(1933). 


Take  a  two  dimensional  field,  with  a  curved  line  of  vorticity, 
representing  the  shed  vortex  wake.  The  wake  vorticity,  y  >is  given 
by  dF/ds  where  T  is  the  bound  vorticity  on  the  wing  and  ds  the  element 
of  length  in  the  wake  plane.  The  potential  <}>  in  the  wake  plane  is 
directly  proportional  to  the  circulation  F  ,  with  scale  constants  which 
do  not  concern  us.  We  assume  the  wake  is  symmetrical  about  the 
vertical  centerline  and  consider  only  one  side  as  shown  in  Fig.  5-6. 
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The  first  term  vanishes  at  the  limits  of  integration  since  at  the 
centerline  r  =  0  and  at  the  tip  =  0  .  Now  we  note  that  dr/ds  = 
008(9-0*),  so  obtain 

M  =  -2  M 

o  p 

This  shows  that  the  Betz  formulation,  which  is  frequently  simpler 
analytically,  does  in  fact  express  the  impulsive  moment.  The  impulsive 
force  and  moment  can  be  transferred  to  other  centers.  When  taken  about 
the  centroid  of  shed  vorticity  the  inpulsive  moment  corresponds  to 
the  angular  "momentum"  of  the  flow  about  that  point  and  is  a  measure 
of  the  "swirl"  in  the  vortex  core  and  surrounding  flow. 

Global  Invariants 

We  have  discussed  the  determination  of  lift,  drag,  and  moment 
of  momentum  for  a  viscous  vortex  system.  These  three  quantities  must 
be  invariant.  We  note  further,  that  if  vorticity  has  not  yet  reached  the 
centerline,  as  is  the  case  during  the  early  core  development,  then  this 
total  vorticity  on  each  side  is  conserved  even  though  it  diffuses  radially 
under  viscous  and  turbulent  influences. 

If  the  total  vorticity  on  one  side  is  conserved,  it  can  be  shown 
by  arguments  in  the  Trefftz.  Plane  (Munk,  1924)  that  the  centroid  of 
this  vorticity  (or  the  mean  vortex  span)  must  be  fixed  for  the  lift  to 
remain  constant.  Thus,  we  can  summarize  the  results  to  state  that 
the  following  four  properties  must  be  conserved  along  the  wake  on 
ru.ch  side  of  the  centerline  at  least  within  a  few  hundred  spans  down¬ 
stream  of  the  wing: 

a.  the  total  vorticity, 

b.  the  centroid  of  vorticity, 

c.  the  moment  of  momentum, 

d.  the  axial  force. 
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As  explained  previously,  after  mixing  or  interaction  between  the 
vorticity  generated  on  each  aide  of  the  wing,  then  only  conditions  (c) 
and  (d)  will  still  hold  as  an  expression  of  drag  and  moment  invariance. 
Although  the  lift  will  still  be  conserved,  it  can  no  longer  be  expressed 
by  the  statements  (a)  and  (b). 

It  is  of  interest  to  note  that  as  vorticity  diffuses  through  the  cell, 
it  ie  annihilated  at  the  centerline.  Then,  for  lift  to  be  conserved,  the 
centroid  of  the  remaining  vorticity  must  move  outboard.  We  Bee  that 
this  is  consistent  with  the  requirement  of  conserving  vorticity  polar 
moment,  since  although  the  total  vorticity  is  reduced,  it  occupies  a 
larger  area. 

In  general  terms  the  drag  (or  axial  force)  can  be  determined 
by  applying  the  momentum  equation  to  a  normal  plane  a~d  the  impulse 
moment  calculated  by  taking  the  second  moment  of  vorticity.  To  de¬ 
termine  the  axial  force  both  pressure  and  velocity  must  be  known.  How¬ 
ever,  since  the  flow  ie  in  general  not  homoenergetic  (because  of  viscous 
dissipation)  the  pressure  cannot  be  inferred  from  the  velocity.  Thus 
both  pressure  and  velocity  must,  in  principle,  be  known  for  the  momen¬ 
tum  theorem  to  be  applied. 

This  difficulty  is  avoided  by  a  convenient  approximation,  valid 
when  the  core  is  still  fairly  compact  and  circular  in  shape. 

Representation  of  the  Prut;  auu  Moment  Integrals 

for  Vortex  Cores 

The  standard  technique  of  simplification  is  to  assume  the  region 
to  be  divided  into  two  portions,  an  outer  portion  where  the  flew  is 
ir rotational  so  that  the  total  head  is  conserved  and  where  the  axiai 
perturbation  ie  negligible,  and  an  inner  core  region  where  it  is  assumed 
that  the  flow  field  is  axisyromctric.  consisting  of  a  tangential  and  an 
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axial  component.  In  this  core  region  there  is  appreciable  dissipation 

due  to  viscosity  and  the  Bernoulli  equation  cannot  be  used  to  determine 

pressure  from  velocity.  Fortunately,  the  assumption  of  axisymmetric 

flow  permits  one  to  couple  the  tangential  velocity  at  a  given  radius  to 

2 

the  radial  pressure  gradient  by  use  of  the  equation  dp/dr  =  v  /r  , 
where  v  is  the  tangential  velocity  at  radius  r  . 

It  will  be  seen  that  the  assumption  of  irrotationality  in  the  outer 
flow,  and  of  axisymmetry  in  the  core  flow,  permits  us  to  infer  the  pres 
sure  from  the  velocity  in  both  regions.  It  is  instructive  to  consider  the 
consequence  of  this  assumption,  which  clearly  shows  the  coupling  of 
the  axial  and  tangential  velocities  due  to  viscous  effects  as  analyzed  by 
Batchelor  (1964).  Figure  5*  7  shows  the  inviscid  velocity  distribution, 
v  ,  due  to  a  potential  vortex  with  a  linear  core  distribution.  In  the 
inviscid  outer  field  we  assume  that  the  axial  flow  W  is  at  the  free 
stream  level,  and  the  total  head  is  conserved.  Thus  the  static  pres¬ 
sure  can  be  calculated  directly  from  the  rotational  velocity  compon¬ 
ent.  For  the  core  region,  we  assume  that  the  streamlines  are  cir¬ 
cular;  thus  here  the  pressure  can  be  determined  (subject  to  a  constant) 
by  centrifugal  equilibrium,  depending  only  on  the  core  tangential 
velocity.  The  constant  is  determined  by  matching  with  the  outer  flow. 
This  is  shown  in  Figure  5-7, 

We  now  assume  that  the  total  head  in  the  core  is  nearly  the  same 
as  that  in  the  outer  flow,  as  might  be  the  case  during  the  early  core 
rollup.  We  then  find  from  Bernoulli's  equation,  that  the  axial  velocity 
in  the  core  must  increase  to  conserve  total  head.  This  will  create  a 
jet -like  axial  wake,  that  is  a  velocity  away  from  the  wing  if  die  wing 
were  moving  in  still  air,  /.a  the  tots!  head  in  the  core  is  reduced  due 
to  dissipation,  this  jet-like  flow  will  reduce  and  finally  become  a  vel¬ 
ocity  deficit  or  wake-like  flow. 


Superimposed  on  this  flow  is  the  wake-like  flow  associated  with 
the  profile  drag  of  the  wing, and  represented  by  the  boundary  layer  shed 
from  the  trailing  edge.  It  is  assumed  that  this  spanwise  vorticity  and 
the  streamwise  lifting  vorticity  are  "attached"  to  the  same  particles 
in  the  flow  leaving  the  wing,  and  is  thus  swept  into  the  core  during  the 
rollup  process.  Consequently,  reduction  in  total  head  is  present  even 
before  there  is  any  dissipation  due  to  motion  in  the  core.  Thus  the 
profile  drag  of  the  wing  (a  wake-like  component)  will  reduce  or  even 
eliminate  the  jet-like  flow  at  the  center  of  the  vortex  core. 

In  general  terms,  then,  the  swirling  flow  due  to  lift,  and  repre¬ 
sented  by  induced  drag,  develops  the  initial  jet-like  flow;  while  the 
profile  drag  develops  the  wake-like  flow.  As  the  core  develops,  dis¬ 
sipation  produces  reductions  in  head,  causing  the  axial  flow  to  become 
more  and  more  wake-like. 

When  the  core  becomes  very  large,  the  centrifugal  pressure 
gradients  redr  e  and  the  final  flow  state  is  one  of  a  constant  pressure 
wake,  with  both  the  profile  and  the  induced  drag  appearing  as  a  velocity 
defect  wake. 

Of  course,  the  magnitude  and  rate  of  change  of  axial  flows  depends 
on  the  ratio  of  induced  to  profile  drag  and  the  dissipative  process.  It  is 
not  certain  whether  the  latter  is  laminar  or  turbulent.  However,  we  see 
that  this  can  qualitatively  account  for  the  observations  of  a  jet-like  flow 
near  the  wing  developing  later  into  a  normal  wake-like  axial  flow. 

This  model  so  far  conserves  lift  by  conserving  circulation  and 
vortex  span.  We  note  now  that  if  we  assume  vorticity  and  total  head 
loss  are  uniform  in  the  core,  then  the  additional  two  global  relations 
of  impulse  moment  and  drag  arc  sufficient  to  determine  the  core  radius 
and  the  head  lose.  Thus,  in  a  crude  fashion,  one  can  obtain  the  major 
characteristics  of  the  viscous  vortex  field,  which  are  the  swirl  number 


and  scale  of  the  axial  flow,  Unfortunately,  such  an  analysis  can  not 
give  any  streamwise  variation  in  the  core  parameters.  To  obtain  an 
insight  into  the  core  development  with  time  (or  axial  distance),  more 
parameters  must  be  introduced  to  relax  the  core  assumptions  of  uni¬ 
form  vorticity  and  uniform  head  loss. 

An  interesting  approach  to  this  was  made  by  Mager  (1972).  He 
considered  a  single  vortex  only,  and  ignored  all  but  the  core  flow. 

Then  he  assumed  a  normalized  distribution  for  the  tangential  velocity 
with  no  free  parameters,  and  a  normalized  distribution  for  the  axial 
velocity  with  one  free  parameter,  <X  ,  scaling  this  velocity.  In  this 
model,  of  course,  angular  momentum  of  the  core  is  not  conserved, 
but  is  continuously  reduced  by  viscous  torques  on  the  circumference, 
taken  by  the  author  to  be  laminar.  Thus  a  core  development  with  ax¬ 
ial  distance  can  be  predicted. 

Mager  showed  that,  according  to  his  assumptions,  for  given 
initial  conditions  of  angular  momentum  and  axial  drag,  there  were 
in  general  two  possible  flow  states,  provided  certain  conditions  were 
met.  The  one  state  was  characterized  by  a  larger  core  and  lower  ax¬ 
ial  flow  speed  than  the  other.  Moving  downstream,  subject  to  viscous 
attenuation,  these  states  converged  to  a  single  solution.  Beyond  this 
point  no  solution  was  possible.  He  defined  the  point  at  which  no  solu¬ 
tions  (in  his  assumed  similarity  form)  existed  as  the  critical  point, 
and  derived  a  curve  of  critical  axial  velocity  as  a  function  of  swirl 
number.  This  curve  has  a  very  similar  character,  and  is  aiso  quite 
close  numerically,  to  that  given  by  Benjamin. 

In  order  to  continue  the  flow  beyond  this  critical  point  (at  which 
regular  solutions  vanished),  Mager  introduced  a  further  parameter 
into  the  tangential  velocity  profile.  With  this  parameter  added,  it 
was  found  that  the  flow  beyond  the  critical  point  exhibited  a  very  large 
increase  in  core  size  (about  doubling)  and  a  severe  reduction  in  axial 
velocity,  such  that  the  flow  near  the  centerline  was  close  to  stagnation. 
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Mager  postulated  that  this  model  accounted  for  the  observed  vor¬ 
tex  breakdown  phenomena,  the  two  solutions  In  the  regular  laminar  state 
representing  the  bulging  and  contraction  of  the  axiaymetric  bubble,  the 
vanishing  of  the  axisymmetric  solution  representing  the  spiral  insta¬ 
bility,  and  the  final  axisymmetric  solution  (with  the  large  core  and 
large  axial  flow  deficit)  representing  the  vortex  breakdown. 

It  is  clear  that  the  actual  critical  values  obtained  by  Mager  are 
dictated  by  the  form  of  the  axial  and  tangential  profiles  assumed.  How¬ 
ever,  the  forms  used  are  properly  continuous  and  appear  fairly  reason¬ 
able,  thus  other  selections  would  probably  not  greatly  affect  the  numeri¬ 
cal  results.  However,  the  vanishing  of  the  solutions,  a  key  element  in 
hypothesizing  the  spiral  instability  followed  by  breakdown,  apparently 
occur  because  of  the  limited  number  of  parameters.  Additional  profile 
parameters  could  be  introduced,  which  would  continue  the  real  solutions 
beyond  Mager's  critical  value.  It  is  believed  that  it  would  he  valuable 
to  make  an  extension  of  Mager's  technique,  extending  the  analysis  from 
a  single  axisymmetric  vortex  core  to  the  vortex  pair  associated  with 
an  aircraft  vortex.  Here  one  would  take  into  account  the  force  and 
moment  contributions  of  the  outer  inviscid  flow  and  investigate  numeri¬ 
cally  the  significance  of  further  parameters  on  the  disappearance  of 
solutions.  It  should  be  noted,  too,  that  Mager’s  solution  does  not  con¬ 
sider  unsteady  flows  so  any  wave-like  instabilities  are  excluded. 

In  this  light  we  refer  to  a  new  result  by  Bilanin  and  Widnaii  |19?>) 
in  which  the  unsteady  core  instability  is  treated.  Here  the  authors  state 
that  the  Crow  instability  induces  unsteady  axial  pressure  perturbations 
which  in  turn  cause  the  s^ortex  core  to  undergo  breakdown.  According 
to  this  theory  the  breakdown  will  occur  where  the  two  vortices  are 
furthest  apart. 

This  phenomenon  has  been  observed  in  lab  tests  and  also  appears 
to  occur  la  some  of  Tombaeh’s  flight  tests.  It  is  believed  that  this  is 
another  unsteady  way  of  exciting  vortex  breakdown  and  is  thus  not  in¬ 
consistent  with  the  other  theories  quoted. 
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We  note  that  this  does  not  imply  that  sinuous  instability  is  a 
necessary  condition  for  core  bursting.  In  fact,  the  movies  of  Tombach 
(19?3)  show  many  cases  where  core  bursting  is  clearly  occurring  on 
a  trailing  vortex  which  is  essentially  rectilinear. 

Conclusion 

The  nature  and  precise  mechanism  of  vortex  breakdown  is  still 
controversial.  It  is  generally  served  that  the  breakdown  in  always 
associated  with  adverse  pressure  gradients,  apparently  these  may  be 
either  cause  or  effect.  The  conditions  for  breakdown  to  occur  are  re¬ 
lated  principally  to  the  swirl  ratio  and  the  magnitude  of  the  axial  flow. 
No  general  agreement  on  this  critical  function  has  been  reached  but  both 
Benjamin  and  Mager  give  similar  results  which  can  be  expressed  as 
Figure  5-8.  Here  we  have  expressed  the  axial  velocity  as  the  ratio  of 
the  mean  core  axial  velocity  to  the  freestream  flow.  We  stress  that 
although  scales  are  given  this  ie  only  a  representative  sketch. 

It  should  be  noted  that  both  supercritical  and  subcritical  states 
are  driven  by  dissipation  towards  the  critical  condition  so  that  in 
general  most  cores  will  approach  the  breakdown  state.  We  note  that 
Mager  defines  a  further  dividing  line  in  the  R,  W^/W  diagram  above 
which  no  breakdown  is  possible. 

If  proper  conditions  for  breakdown  were  knoun,  then  it  would 
still  be  necessary  to  be  able  to  calculate  core  development  to  deter¬ 
mine  the  core  state. 

One  of  the  controversial  factors  of  this  calculation  is  the  viscous 
transfer  cons-iant  and  whether  it  is  turbulent  or  laminar.  This  is 
certainly  extremely  important  in  any  prediction  scheme  since  turbu¬ 
lent  dissipation  would  cause  breakdown  an  order  of  magnitude  more 
rapidly  than  laminar  flow. 
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6.  METEOROLOGICAL  ASPECTS 


6. 1  GENERAL 


In  this  section  of  the  report  w a  describe  flow,  stability,  and  tur¬ 
bulence  characteristics  in  the  layer  up  to  abotit  ?0m.  Formulas  inter- 
relsting  factors  in  thif  layer  are  presented  and  discussed.  Then  a 
data  system  is  presented  which  c:  n  be  a  starting  point  for  the  design 
of  %  S',  stern  to  monitor  the  atmosphere  for  the  evaluation  of  techniques 
Ut  monitor  and  fore  ast  wake  behavior.  An  eventual  operational  sys¬ 
tem  will  be  simpler  than  the  research/cvuxuation  system. 

The  70m  height  ar  ses  from  considering  that  the  greatest  danger 
from  vortex- wake  encounters  occurs  at  low  altitudes  where  there  is 
too  little  height  (time)  for  the  encountering  »  rcraft  to  recover  from 
an  unusual  altitude.  For  the  standard  medi  m  and  large  jet  transports, 
it  is  estimated  that  thi?  critical  heigh;  can  reach  uo  to  about  70m.  aad 
of  course  me  danger  gets  greater  as  the  air-raft  height  above  ground 
decreases.  Any  exact  assessment  i.  ‘  the  hazard  is  complex  and  depends 
on  many  quantifier..  This  height  ale.  coincides  '•oximataiy  with  the 
height  along  tha  glide  slope  *'  the  mb  ne  nut^»r  of  »>  jj  $  system. 
Concentration  on  such  xn  area  as  a  test  rea  is  .use  u.«.  because  many 
vortex  sources  are  ava  lab’  >.  and  the  area  be*  th-  hazardous  region 
in  the  landing  approach  when  the  vortices  from  lo'w-fiyulg  and  slow- 
flying  generating  aircraft  are  strong  end  the  haissrd  to  subsequent  siow- 
ftyiag  aircraft  is  high. 


The  AtmosphesU  Environment 


Tne  lowest  1  km  a,  the  3  mo;iohc re,  and  especially  if*e  lowest 
100m,  have  been  intensive’  v  studied  Usa-aswscb  as  this  covers  the  re¬ 
gion  where  hi-.ruans  live,  w/*ere  pollutant*  arc  initially  dispersed,  and 
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'\hi-ri’  i In-  vertical  fluxes  of  momentum,  heat,  and  moisture  arc  strong¬ 
est.  Out  of  (he  studies  there  heve  evolved  various  formulas  for  the 
mean  profiles  and  turbulent  characteristics,  formulas  usually  base*)  on 
a  theoretical  treatment  tailored  by  empirical  results.  To  summarize 
the  situation,  the  overall  understanding  of  the  characteristics  of  these 
lowest  layers  can  be  deemed  "good”.  In  the  simplest  situations  of  ho¬ 
mogeneous  terrain,  and  constant  meteorological  conditions,  which  do 
not  have  strongly  stable  lapse  rates,  there  are  satisfactory  semi-em- 
pirical  formulas  available  to  represent  conditions.  In  the  more  com¬ 
mon  complex  conditions  the  application  of  the  formulas  is  more  diffi¬ 
cult  and  the  apparent  success  in  applying  the  formulas  probably  re¬ 
lates  more  to  the  fact  that  great  accuracy  is  not  required  by  the  prob¬ 
lem  nor  are  the  data  available  to  show  what  accuracy  was  attained. 


A  General  View 

The  planetary  boundary  layer  is  the  layer  near  the  surface  where 
momentum  and  heat  flux  can  be  large,  i.  e.  ,  where  there  is  a  flux  link 
between  the  upper  flow  and  the  earth's  surface.  The  boundary  layer 
height  can  be  under  0.  1  km  or  over  10  km,  with  1  km  being  a  common 
height  in  daytime.  In  most  typical  form,  the  boundary  layer  is  divided 
into  two  parts:  (1)  the  surface  boundary  layer,  across  which  the  verti¬ 
cal  fluxes  are  constant  (this  constancy  provides  simplifications  which 
yield  rather  simple  equations  for  mean  and  turbulent  properties),  and 
(2)  the  total  boundar/  layer,  whose  top  is  found  where  the  fluxes  be¬ 
come  small. 

The  height  (cm)  of  the  surface  boundary  layer  below  which  the 

stress  magnitude  varies  by  less  than  20%  is  put  by  Lumley  and  Panofsky 

2 

(1964)  at  2000  j  (where  t  is  the  surface  stress  in  dynes/cm  ).  Since  j 
is  typically  of  the  order  1-10  dynes/cm,  h  is  typically  between  2*10^- 

4 

2*10  ern,  or  20  and  200m. 
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6.  2  THE  QUANTITIES  OF  INTEREST 


The  Various  Quantities 

Present  thinking  suggests  that  the  primary  meteorological  vari¬ 
ables  of  importance  in  the  transport  and  decay  of  a  vortex-wake  system 
are:  the  mean  horizontal  wind  field  which  causes  horizontal  drift,  and 
the  turbulence  and  stability  fields  which  play  a  role  in  determining  de¬ 
scent  and  decay  of  the  vortex-wake.  These  wind,  turbulence,  and  sta¬ 
bility  conditions  need  to  be  known  throughout  the  region  occupied  by 
the  wake  during  its  whole  life  until  it  has  decayed  to  the  point  where  it 
cannot  const. ..ute  a  hazard.  Another  meteorological  quantity  will  in 
some  cases  also  prove  to  be  of  primary  importance  --  the  mean  verti¬ 
cal  motion  of  the  local  atmosphere.  Cne  final  property,  the  vertical 
shear  of  the  horizontal  wind,  across  the  aircraft  flight  path,  may  be 
of  importance. 

Evidence  is  accumulating  that  stability  of  the  atmosphere  may 
have  only  secondary  importance  in  either  vortex  decay  or  vertical  de¬ 
scent.  Thus  we  do  not  consider  stability  a  primary  variable,  but  meas¬ 
ure  it  because  of  its  effect  on  determining  turbulence,  drift,  or  the  lo¬ 
cation  of  vertical  shears. 

When  we  talk  about  the  mean  wind,  horizontal  or  vertical,  we  are 
thinking  of  "mean"  as  referring  to  events  taking  place  for  a  minute  or 
two,  the  lifetime  of  the  wake.  Viewed  on  a  longer  scale,  the  same  mo¬ 
tion  might  well  be  considered  as  turbulence.  Such  considerations  of 
scale  and  averaging  are  important  in  the  design  of  an  operational  system. 


Turbulence 


Turbulence  covers  all  wavelengths,  and  multiple  correlations, 
in  the  three  directional  components,  jo  one  major  problem  here  is  to 
decide  what  wavelength  range  and  component  is  most  appropriate  for 
the  problem  of  vortex-wake  decay.  One  powerful  simplifying  assump¬ 
tion  which  has  been  used  by  AeroVironment  Inc.  is  that  the  *uvbulence 
of  significance  is  in  the  inertial  subrange  of  eddy  size  .  Kolmop  orov1  s 
similarity  hypothesis  results  in  the  conclusion  that  all  statistical  pro¬ 
perties  of  turbulence  within  the  inertial  subrange  relate  only  to  e  , 
the  equilibrium  rate  of  turbulent  eddy  dissipation.  This  picture  of 
turbulence  is  so  simple  it  permits  many  important  formulas  to  be  de¬ 
rived  easily  by  dimensional  analysis  and  it  also  permits  useful  em¬ 
pirical  relations  to  be  made  even  when  all  the  physical  connecting 
links  may  not  be  understood.  Kolmogorov's  similarity  hypothesis 
starts  from  a  very  simple  picture  of  turbulence  mechanisms.  It 
pictures  that  turbulence  enters  a  system  continually,  primarily  at 
large  wavelengths.  The  large  eddies  break  down  to  smaller  eddies, 
which  eventually  break  down  to  molecular  motions  through  viscosity 
effects.  The  smaller  eddies  are  so  decayed  from  the  large  "mpu". ' 
eddies  that  the  smaller  eddies  have  "forgotten"  their  an  ^stry  and  sc 
cannot  be  aware  of  direction  --  statistically  their  energy  must  be  iso¬ 
tropic.  There  is  an  "inertial  subrange"  within  the  range  of  eddies 
which  exhibit  isotropy;  the  inertial  subrange  covers  all  tho1  e  except 
the  tiny  ones  which  are  strongly  affected  by  viscosity.  In  the  inertial 
subrange  tho  only  quantity  to  which  any  statistical  property  of  the  flow 
can  be  related  is  s  ,  the  equilibrium  rate  at  which  energy  enters  the 
system,  cascades  through  the  inertial  subrange  of  eddy  sizes,  and  is 
removed  (as  heat)  by  viscous  effects.  Xu  the  atmosphere  well  above 
the  ground  the  inertial  subrange  of  eddy  sizes  extends  typically  from 
a  few  centimeters  to  many  hundreds  of  meters.  Even  though  it  is  an 
idealized  concept,  the  relations  derived  using  it  aie  found  to  be  sur¬ 
prisingly  ^cod  representations  oi  the  atmosphere. 
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The  atmospheric  wavelengths  under  discussion  for  vortex-wake 
nreakup  are  probably  no  larger  than  the  primary  wavelength  for  the 
Crow  looping  instability  {300m  for  a  large  jet),  and  may  even  be  con¬ 
siderably  smaller.  Thus  they  fit  into  the  inertial  subrange  and  e  is 
the  appropriate  turbulent  quantity  to  use.  The  Crow-Kolmogorov 
quantitative  instability  theory  considers  that  vortex  motions  in  a  V- 
shaped  (approximately  90°)  trough  along  the  lint  of  flight  are  the  pri¬ 
mary  instability  initiators,  and  assumes  the  strength  of  such  motions 
to  be  related  to  g  ,  The  theory  assumes  that  symmetrical  perturba¬ 
tions  along  this  specific  trough  surface  are  as  likely  as  motion  of  a 
single  line  in  any  single  plane.  It  is  certainly  possible  that  other 
inodes  of  motion,  including  correlations  between  vortex  lines,  are 
fundamental  to  the  initiation  of  Crow  instability.  However,  ail  inodes 
are  within  the  ineriial  subrange,  and  so  e  is  the  turbulent  quantity  to 
observe  and  to  fit  into  theories  about  vortex  lifeapon  as  related  to  at¬ 
mospheric  turbulence. 

If  we  were  concerned  with  wakes  only  at  altitudes  0;f,  say  300m 
or  more,  thus  we  could  probably  comfortably  use  €  as  the  most  ap¬ 
propriate  turbulent  quantity  and  avail  ourselves  of  the  simplicity  in¬ 
herent  in  inertial  subrange  concepts.  However,  we  are  actually  inter¬ 
ested  in  wakes  at  70m,  50m,  and  often  even  lower.  Studies  have  shown 
(MacCready,  1962,  and  the  discussion  by  Lumley  and  Panofsky,  1904, 
pg.  166-167)  that  near  the  ground  the  spectrum  laws  stem  reasonably 
valid  even  at  wavelengths  about  twice  the  height  of  observation  (and 
still  longer  in  unstable  cases).  It  may  seem  surprising  ;hat  wave¬ 
lengths  greater  than  z  can  actually  show  statistical  energy  isotropy 
and  agree  with  the  simplified  spectrum  laws,  but  it  should  be  remem¬ 
bered  that  in  our  turbulence  formulas  we  are  dealing  with  energies 
and  velocities,  not  final  displacements,  and  so  isotropv  of  energy  for 
the  three  directional  components  for  100m  horizontal  wavelengths  may 
be  reasonable  at  a  height  of  only  30m. 
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In  summary,  the  strength  of  the  inertial  subrange  cone  pt  is  so 
great  in  the  wake  area  that  the  monitoring  or  predicting  of  its  sole  vait- 
able,  e  ,  may  be  quite  useful  for  any  evaluation  program  of  monitoring/ 
predicting  systems.  However,  we  must  realize  the  possible  increasing 
inadequacy  of  the  inertial  subrange  concept  for  wake  decay  pre<  etions 
as  we  approach  the  ground,  and  therefore  also  consider  other  formu¬ 
lations  for  turbulence  there,  especially  the  vertical  turbulence  compon¬ 
ent  at  wavelengths  near  the  Crow  looping  stability. 


6.  3  BOUNDARY  LAYER  RELATIONSHIPS 


The  Mnnin-Obukhov  Similarity  Theory 

The  Monin- Obukhov  similarity  theory  develops  useful  relation¬ 
ships  in  the  surface  boundary  layer  based  on  certain  quantities  which 
are  constant.  Lumley  and  Panofsky  (1964)  provides  a  critical  review 
of  the  concepts. 

Throughout  the  surface  boundary  layer  the  momentum  flux  (and 
hence  the  shearing  stress  T  )  is  constant,  and  the  heat  flux  (H)  is  con¬ 
stant.  Using  these,  together  with  ^  (where  g  is  gravity  and  T  is 
temperature)  which  is  another  dimensional  "constant",  it  is  possible 
to  develop  a  velocity  (u*),  a  length  (L),  and  a  temperature  (T*)  which 
are  also  constant  and  are  three  convenient  dimensions  for  scaling. 
When  the  principal  variables,  such  as  temperature,  ”  ind,  and  height, 
are  expressed  non-dimensionally  as  fractions  of  these  quantitites,  a 
series  of  non-dimensional  equations  results  that  are  of  general  valid¬ 
ity  in  the  surface  boundary  layer  for  conditions  not  too  far  from  neu¬ 
tral.  For  simplicity,  we  neglect  humidity  in  this  discussion  since  its 
effects  are  ordinarily  small. 

The  equations  reduce  to  simple  laws  for  the  neutral  stability 
case  (H  -  0),  but  turn  out  to  be  rather  complex  and  awkward  to  apply 
when  H  /  P  .  Much  attention  has  been  paid  to  empirical  approxima¬ 
tions  to  the  basic  equations,  with  the  final  equations  having  a  general 
form  which  relates  well  to  the  fundamental  theory. 

The  friction  velocity,  u*  .  is  related  to  the  shearing  stress  T 
and  density  o  by 

U*  *  (T/p)*  (6-1) 

Since  we  are  in  the  "constant"  stress  region,  T  is  the  same  as 
?  .  the  surface  stress.  The  wind  direction  is  constant  in  the  surface 

boundary  layer. 
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The  scaling  length  L  is  given  by 


u*J  C  pT/(kgH) 
P 


where  is  the  specific  heat  at  constant  pressure,  and  k  is  a  dimen¬ 
sionless  coefficient  (the  voa  Karman  constant, K=  0.4),  The  dimension- 
less  height  —  turns  out  to  be  especially  useful  for  quantifying  the 
effects  of  non-neutral  lapse  rates.  L  ,  a  length  which  can  be  negative 
as  well  as  positive,  is  basically  a  convenient  lapse  rate  parameter. 

1*^1  ^e^nes  height  below  which  mechanical  turbulence  is  dominant. 


•production  rate  of  convective  enerc 


zkgH 


L  production  rate  of  mechanical  energy  u*CpPT  ' 

Because  H  is  usually  not  measured,  another  quantity  L'  is  often 
used  rather  than  L  . 


L'  =  L 


where  Km  is  the  eddy  viscosity  (exchange  coefficient  for  momentum! 
and  that  for  heat  flux.  K^/K^  is  near  uni.y,  being  larger  in  un¬ 
stable  conditions  and  smaller  in  stable  conditions. 


.  „  8v 
U*T^ 


where  0'  its  potential  temperature.  All  the  quantities  in  Eqn.  (6-5) 
are  fairly  readily  measurable 


Another  measure  of  "non- neutralness"  is  Rf  ,  the  flux  Richard¬ 
son  Number,  which  is  sometimes  involved  in  derivations.  It  if  defined 


as  the  ratio  of  the  production  of  turbulent  energy  due  to  buoyancy  forces 
to  that  due  to  Reynolds  stresses  (mechanical  forces). 


A  gradient  Richardson  Number  is  defined  by 


R.  =  R  K  /K 
i  g  m'  n 


(6-6) 


R.  is  the  ratio  of  buoyancy  to  inertia  forces. 


86' 

R  .  £ 

'  T  $2 


(6-7) 


R^  can  be  ascertained  from  mean  wind  and  temperature  measurements 
at  two  levels  of  a  tower. 

By  various  derivations,  the  logarithmic  wind  profile  is  found  for 
neutral  conditions: 


8u  _  u* 

8z  “  k(  z+  z  ) 
o 


(6-8) 


where  z^  is  the  roughness  height.  For  z»  z  ,  substituting 
Eqn.  (6-8)  in  (6-7)  and  comparing  with  Eqn.  (6-5)  shows  us 


5/L*  =  Ri 


(6-9) 


for  cases  very  close  to  neutral.  The  relation  is  more  complex  as  we 
move  away  from  neutral.  The  more  complete  relations  are  shown  in 
Table  6-1.  For  R^  greater  than  (more  stable  than)  about  0.  1  the  ex- 
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change  between  layers  diminishes  so  much  that  similarity  laws  are  not 
applicable.  For  R^  less  than  (more  unstable  than)  about  -1  the  region 
of  free  convection  is  reached  and  L'  no  longer  is  a  useful  scaling  para¬ 
meter. 


TABLE  6-1.  Ri  -  z/L1  RELATIONSHIP 


z/L' 

Q 

R. 

L 

2 

TT 

R./(1-7R.) 

Si 

rt 

0.  1 

0.  3 

”  -  —  -OJ 

0.  01 

0.  01 

R. 

Neutral 

0.  00 

0.  00 

1 

—  —  —  — -  j 

-0.  01 

-0.  01 

1 

I 

-0.  1 

-0.  07 

1 14 

Ri/(l-18Rl)w* 

Si  j 

3 

03 

-0.  3 

-0.  z 

e  1 
D  ! 

-L0 

-0.  4 

The  Wind  Profile 

The  gradient  foi  n  of  the  wind  profile  valid  for  neutral  conditions 
ha»  already  been  given  in  Eqn.  (6-8).  Integrated,  this  becomes 


z+z 

In (-r-£) 


(6-10) 


f;  • 
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which  is  usually  compromised  to 


U  =  ^  cn  (6-11) 

K  2 

o 

since  z>>  2^  where  the  formula  is  usually  applied. 

Fiedler  and  Panofsky  (IbTZ)  summarize  information  o".  the  val¬ 
idity  of  the  logarithmic  wind  profile.  They  point  out  that  Eqn.  (6-11) 
may  be  valid  only  up  to  30m  or  so.  and  offer  a  more  accurate  form 
which  should  serve  up  to  100m: 

U  -  t  144  fz  (6- 11a) 

k  2 

a 

where  f  is  the  Coriolis  parameter  (f  -  2X  radians  per  second  rotation 
of  earth  X  sine  of  latitude).  New  evidence  suggests  the  constant  144 
may  be  large.  In  any  case,  in  the  present  review,  for  simplicity  we 
will  use  Eqn.  (6-11)  rather  than  (6-1  la). 

The  surface  from  which  s  is  measured  is  the  ground,  for  z 

o 

small,  but  should  be  elevated  when  z  is  large.  For  estamole.  in  a 
forest  the  "surface”  may  be  set  at  about  2/3  the  height  of  the  trees  -- 
about  where  obstruction  density  is  maximum. 

If  U  vs  z  is  known  at  two  heights,  u4  and  z  can  both  be  found 

o 

from  Eqn.  (6-11).  For  accuracy,  the  calculation  is  often  made  from  a 
more  complete  wind  profile.  ?o  is  a  property  of  the  local  surface  -- 
a  coustant  for  all  conditions  (although  £t  a  given  point  it  will  vary  with 
wind  direction  since  it  relates  to  upwind  conditions).  With  a  known. 
Eqn.  (6-11)  gives  the  whole  V  v*  s.  profile  throughout  the  layer  cor¬ 
responding  to  one  velocity  U  at  a  particular  measurement  height  *  . 


For  non- neutral  conditions  there  are  various  approaches  to  modi¬ 
fications  to  Eqn-  {6-1 1)  based  on  the  Monifi-ObukLov  concepts,  yielding 
a  log-linear  curve.  The  most  complete  form  of  the  equation  is 

U  =  [  in  — —  <£t)  }  (6-12) 

"o 

2 

where  the  universal  function  ¥(-£j)  is  presented  graphically  by  Lurnley 
and  Panofsky  (1964)  in  Figure  6-1  reproduced  on  the  following  page. 

Equation  (6-12)  covers  adequately  from  the  most  stable  condi¬ 
tions  for  which  the  M.oain- Obukhov  theory  can  hold  (about  R.  -  0.  1)  to 
strongly  unstable  cases,  say  <  -5  . 

Mote  that  for  conditions  very  close  to  neutral  (  R.  <  0.  01)  y  is 

i 

given  explicitly  in  the  legend  of  Fig.  6-1.  For  greater  stability,  up  to 
R.  -  0.  1,  the  explicit  equation  for  ¥(^>)  is  still  4.  3  •  although 

there  is  some  controversy  about  the  coefficient  which  may  be  as  large 
as  7.  0  instead  of  4.  5. 

The  above  relations  are  derived  for  forced  convection,  where 
the  vertical  transfer  of  heat  and  momentum  is  more  from  mechanical 
turbulence  l, Van  by  heat  convection.  Other  results  are  obtained  for 
free  convection,  when  the  relative  effect  of  heat  convection  ‘s  larger, 
say  R.  -  -l  .  when  mechanical  and  heat  energy  are  produced  at  about 
Use  same  rate.  Eqn.  (6-12)  is  still  a  fair  interpolation  formula  for  this 
condi U.-n.  Heat  energy  results  in  much  more  efficient  vertical  transfer 
of  properties  than  is  the  case  for  mechanical  ets?rgy;  the  connective 

eddies  are  larger,  with  more  vertical  continuity.  For  free  con v*c-  j 

Ron  the  Vfonin- Obukhov  approach  is  net  used  fcust  a  tieffcrenl  dimaa- 
tcooi  argument  (mix i-ag  length  theory)  gives? 


d-n 


23?-? 


)*)^.T1jU#Xyr 'VXl  r 


Figure  6-1.  Universal  Function  *  for  the  Integrated  Vind 

Profile.  Main  Graph  for  Unstable  Stratification;  Insert, 
for  Stable  Conditions,  For  |z/L>'|  or  |R.|<0. 01, 
1=4,5  z/L'  =  4.  5  R.  . 
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(6-13) 


-i/3 

)  * 


■4/3 


if  K.  /K  is  constant  with  height,  Eqn.  (6-13)  shows 
“  m-i/3 

and  U  "  z  ,  predictions  which  are  reported  to  fit  wind  observations 
on  convective  days  with  moderate  wind  speeds  quite  well.  In  strong  con¬ 
vection  with  little  wind,  the  U  profile  is  rather  indeterminant. 


Another  approach,  the  power  law  method,  has  been  found  con¬ 
venient  for  describing  wind  profiles,  especially  for  stable  conditions. 
It  has  no  theoretical  basis,  out  has  a  flexibility  and  simplicity  which 
makes  tt  useful  for  tailoring  to  empirical  data. 


U 

"°T 


t“) 


(6-H) 


where  the  subscript  (I)  refers  to  the  conditions  at  a  reference  level 
end  p  is  an  exponent  varying  between  0  and  1.  Panafsky.  Biaekadar. 
ana  McVehil  (1940)  shoved  p  could  be  derived  from  s.  L*.  and  z 

Q 

ann  that  she  resulting  values  of  p  were  consistent  with  numerous  ob¬ 
served  values,  p  is  not  a  constant  with  height  of  the  layer  concerned; 
it  increase*  with  height  in  stable  cases  and  decreases  in  height  in  un¬ 
stable  case*.  It  i*  found  to  deper.-.  little  upon  wind  speed;  it  in¬ 
creases  slowly  with  increasing  roughness,  and  varies  primarily  with 
fapsc  rate  {decreasing  as  the  lapse  rate  becomes  more  unstable). 


Pasquill  (1962)  suggested  certain  stability  classes  which  have 
proven  convenient  for  wind  profile  and  diffusion  studied, 
cefaaed  lr  Table  fe-2  on  the  following  p&g©. 


TABLE  6-2 


RELATION  OF  TURBULENCE  TYPES 
TO  WEATHER  CONDITIONS 


A —  Extremely  unstable  conditions  D —  Neutral  conditions* 

R —  Moderately  unstable  conditions  E — Slightly  stable  conditions 
C — Slightly  unstable  conditions  F — Moderately  stable  conditions 


Surface  wind 
speed,  m/eec 

Daytime  insolation 

Nighttime  conditions 

Thin  overcast 

or  *V% 

cloudinesst  cloudiness 

Strong 

Moderate 

Slight 

<2 

A  • 

A-B 

B 

2 

A-R 

B 

C 

E 

F  r 

4 

B 

n-c 

C 

D 

E 

6 

C 

C-D 

D 

D 

D 

>(> 

C 

D 

D 

D 

D 

•Applicable  to  heavy  overcast,  day  or  night. 

tThe  degree  of  cloudiness  is  defined  as  that  fraction  of  the  sky  above 
the  local  apparent  horizon  which  is  covered  by  clouds. 


These  classes  can  be  used  to  summarize  data  as  follows: 


TABLE  6-3.  DATA  SUMMARY 


Pasquill  Class 

6T/6z 

«L 

<P 

A 

^-0.  6  °C/100' 

25° 

.  15 

B 

-0.  5 

20° 

.  17 

C 

-0.  4 

15° 

.  20 

D 

-0.  2  to  -0.  3 

10° 

.26 

E 

-0.  1  to  0.  3 

5° 

.  39 

F 

0.  4  to  1.1 

2.  5° 

00 

• 

(G) 

M.2) 

n.7°) 

The  6T/6s  criteria  come  from  USAEC  Safety  Guide  23  (1972),  and  are  r 
lated  to  Pasquill  classes  for  moderate  winds.  They  are  measured  be¬ 
tween  10  and  45  m.  The  6,^  (vertical-axis  vane  RMS)  values  are  from 
Safety  Guide  23  and  Gifford  (1968),  and  Slade  (1965).  The  p  values 
come  from  DeMarrais  (1959)  and  apply  to  the  rough  Brookhaven  site 
{zq  ~  0.  5  to  1 ,  0m). 

Golder  (1972)  provides  a  more  rational  way  of  relating  Pasquill 
classes  to  the  quantities  on  which  the  Moriin-Obukhov  laws  are  based. 

He  shows  how  the  Pasquill  classes  are  functions  of  z^  and  L  (read  L' 
or  rather  .1  /IT  for  convenient  plotting).  Figure  6-2  is  taken  from 
Golder' s  report. 


Figure  6-2.  1/L  as  a  Function  of  Pasquill  Classes  and 
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Since  use  of  zq  and  L1  permits  one  to  draw  conveniently  on  the  ration¬ 
al  and  versatile  equations  based  on  Mc-nin-ObuSchov,  i*  seems  prefer¬ 
able  to  use  them  instead  of  the  empirical  power  laws  which  vary  with 
site  and  altitude. 

g  vs  Height 

The  dissipation  rate  g  ,  under  steady  state  conditions,  is  equal 

to  the  mechanical  production  of  energy  u*  6U/hz  plus  the  buoyant 

production  gH/{C  pT)  less  the  flux  divergence  l/p(5F „  /bz) 
p  Is* 

m 

where  Fg  is  the  upward  flux  of  total  turbulent  kinetic  energy  per 
m 

unit 


There  is  considerable  controversy  about  magnitude  of  the  flux 
divergence  terrr.  The  observations  and  theory  cited  by  various  auth¬ 
ors  for  unstable  conditions  vary  from  the  concept  that  the  flux  diverg¬ 
ence  term  fully  cancels  the  buoyant  production  term  to  the  concept 
that  the  flux  divergence  term  is  negligible.  Thus  all  formulations  fit: 


(6-15) 


with  ti  '  value  of  the  coefficient  being  between  0  and  1.  G.  Briggs  (at 
NOAA  Environmental  Research  Laboratories,  Oak  Ridge,  Tennessee) 
has  recently  reviewed  all  the  available  evidence  and  concluded  that 
Cj  "1/3  gives  the  best  fit  with  observations,  and  that  with  this  value 
Eqn.  (6-15)  works  reasonably  up  to  2/3  of  the  total  depth  of  the  mixed 
region. 


Using  the  methods  discussed  by  Hanna,  Hutchison,  and  Gifford 
(1969),  and  setting  =  1/3,  Eqn.  (6-15)  can  be  transformed  to 
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(6-16) 


1-T  R.(l-18R.r 

(1-18R.)* 

It  should  be  noted  that  the  turbulence  intensity  function  common- 
1  /3 

ly  of  interest  is  e  ,  not  e  ,  which  minimizes  the  accuracy  demands 
on  the  e  vs  z  equations.  Taking  the  neutral  approximation,  with 
IL  =  0  ,  Eqn.  (6-16)  becomes 

£3  =u*k_i  z~*  (6-17) 

If  we  have  already  ascertained  zq  and  we  measure  at  height  z  , 
then  from  Eqn.  (6-11) 


u*  =  Ujk/(in  z^/zq) 


(6-18) 


and  so 


3 

U.k3 


z3fn  z,  /z 
1 '  o 


(6-19) 


The  inverse  relation  between  e  and  z  is  not  inconsistent  with  the 
e  observations  from  many  sources  summarized  by  Ball  (1961).  Ball 
examined  data  from  10cm  to  10  km,  covering  various  roughnesses, 
winds,  and  stabilities,  and  although  his  summary  plot  shows  a  wide 
scatter  the  ground-related  turbulence  points  nave  a  z  tendancy. 


t  In  tiie  next  section,  where  we  consider  total  vertical  energy, 
there  is  further  discussion  about  obtaining  in  non-neutral  con¬ 

ditions.  The  conclusion  is  that  Eqn,  (6-17)  is  probably  adequate  in 
conditions  well  away  from  neutral. 
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An  excellent,  recent  review  of  e  vs  z  Is  given  by  Pasquill  (1972). 
One  of  his  main  points  is  that  surface  roughness  variability  in  the  hori¬ 
zontal,  which  makes  even  many  field  research  sites  "non-ideal",  can 
have  a  bigger  effect  on  the  e  vs  z  relationship  than  the  difference  be¬ 
tween  specific  formulas. 


Energy  vs  Height 

Lumley  and  Panofsky  (1%4)  state  that  for  neutral  conditions  in 

the  surface  boundary  layer,  the  total  turbulent  energy  e  should  be 
2 

proportional  to  u*  and  thus  be  given  by 

e  =  CU"V(fn  z/z0)2  ((>-2.0) 

whore  the  dimensionless  constant  C  is  of  order  one.  The  few  avail¬ 
able  data  are  not  in  exact  accord  with  Eqn.  (6-20),  and  the  authors 
note  that  although  the  nearby  z  may  determine  the  vertical  motions, 
the  z  far  upstream  can  have  an  important  bearing  on  the  other  com¬ 
ponents. 

As  to  the  vertical  component  of  turbulent  energy,  0^  ,  contain¬ 
ing  all  frequencies,  the  Monin-Obukhov  theory  indicates  for  a  range 
of  stabilities  that 


a  =  Au*  (6-21) 

w 

where  A  has  a  value  probably  near  1.3  in  the  neutral  case,  but  slightly 
different  values  for  other  IE  or  z/L  .  Combining  Eqns.  (6-8)  and 
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(6-21)  gives,  for  neutral  situations, 

kAU 
w  ~  in  z/z 

o 

In  stable  air  decreases  with  height;  in  unstable  air  it  slowly  in¬ 
creases. 

The  value  of  A  found  by  Monin  (1962),  corrected  by  X  1.  5  as  dis¬ 
cussed  by  Lumley  and  Panofsky  (1964),  is  given  for  approximately  vari¬ 
ous  stabilities  in  termsof  z/L'  .  The  form  shown  by  Panofsky  (1970) 
is  somewhat  different,  but  since  it  represents  more  data  we  will  use 
it.  Putting  it  in  equation  form: 

A  =  1.  3  -  0.  3  z/L,1  (6-23) 

This  is  a  fair  approximation  for  -1.2<  z/L'  <0.  3  , 

If  one  uses  a  vertical  angle  vane  (horizontal  axis)  to  measure 
turbulence,  its  RMS  vaiue  0g  in  radians  will  be  related  to  the  RMS 
of  vertical  velocity  by 

o.  so  /U  (6-24) 

©  w' 

if  we  assume  angles  are  small  so  the  sine  is  the  same  as  the  angle. 
Then  from  Eqna.  (6-22)  and  (6-24)  we  see,  for  the  neutral  case, 

Oq  =  0.  5/fn  z/zq  (6-25) 


0.  5U 
in  z/z 


(6-22) 


6-21 


In  neutral  stability,  (or  any  strong  wind  condition)  is  independent 
of  U  .  Also,  note  that  by  measuring  Og  at  a  certain  height  z^.z  is 
readily  obtained  from  Eqn,  (6-25). 


Equations  which  require  knowledge  of  l/  may  not  be  particular¬ 
ly  convenient  since  L/'s  calculation  from  requires  observing  both 
wind  and  temperature  at  two  heights  and  performing  long-term  aver¬ 
aging  to  obtain  the  needed  gradients.  In  the  practical  case  it  may  be 
more  suitable  to  estimate  H  or  l/  or  R.  .  This  is  what  AeroViron- 
ment  Inc.  did  in  developing  its  General  Concentration  Model  whose  ap¬ 
plication  is  described  by  the  AV  staff  (1972).  They  used  the  method 
suggested  by  Panofsky  and  McCormick  (i960),  who  postulated  on  the 
basis  of  similarity  theory  as  applied  to  the  constant  flux  layer  that 
°w  should  be  a  function  of  height,  the  rate  of  energy  supplied  by 
mechanical  turbulence,  and  the  rate  of  supply  of  convective  energy. 
This  approach  has  the  further  advantage  of  covering  the  free  convec¬ 
tion  range  as  well  as  stabilities  closer  to  neutral.  The  basic  equation 


a  =  A, 

w  2 


z(u*2|y  +  -MH 

0z  pCpT  'J 


(6-26) 


where  A2  and  6  are  constants  (assumed  to  be  -  1.25,  6  -*•  2.4). 

Using  Eqns.  (6-12)  and  (6-26),  one  than  has  <r  as  a  function  of  H  , 

Y(z/L'),  zq,  and  U  .  Next,  various  forms  of  stability  dependent  Y 

as  proposed  by  Panofsky  et  al  (I960)  and  McVehil  (1962)  were  assumed 

in  Eqn.  (6-26).  The  scale  length  L1  is  related  to  heat  flux  and  friction 

velocity  (eqn.  (6-2)).  The  friction  velocity  further  couples  the  mean 

wind  and  L1  through  the  diabatic  wind  equation  (Eqn.  6-12)),  The 

unique  way  in  which  these  parameters  are  related  to  one  another  makes 

it  possible  to  determine  0  as  a  function  of  wind  for  given  values  of  z 

w  o 

and  H  ,  although  not  in  closed  form.  The  results  were  plotted  by  AV 

in  the  form  of  o  vs  U  at  10m  for  families  of  z  and  H  values.  The 
w  o 
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roughness  dependence  was  removed  by  using  a  roughness  reduced  vel¬ 
ocity  Uy  =  Uo(£filO/zQ)/(2n  10/zr)  where  Ur  is  the  wind  speed  at  10m 
corresponding  to  the  roughness  z^  .  The  "universal  roughness  re¬ 
duced  curves"  provided  a  simple  summary  for  computer  storage  giving 

■  /3 

quantitatively  a  =  DUj,  for  the  appropriate  strong  winds,  a  =  CH' 

W  W3  w 

for  H>  0  and  U  <BH  (positive  H  but  very  low  winds),  and  a  -  0 

r  1/3  w 

for  H  <  0  and  U  <  -AH  (negative  H  and  very  low  winds).  For  U 

”  3 

and  0  in  m/sec,  z  in  m,  and  H  in  cals/cm  /min,  the  constants 
w  r 

are  A  =  13. 2,  B  =  5.  7,  C  =  0. 87,  and  D  =  0.  125. 


The  similarity  in  the  derivations  of  Eqns,  (6-21)  and  (6-26), 
(and  the  uncertainty  in  each)  suggests  that  one  will  not  have  much 
larger  errors  if  one  sets 


\  A3U* 
£  3  =  "X  T" 
k3  g3 


(6-27) 


where  the  coefficient  is  assumed  to  be  a  weak  function  of  z/L1  . 
The  variation  of  A3  should  be  less  than  the  cube  root  of  the  variation 
of  A  with  z/L'  (Eqn.  (6-23))  since  the  coefficient  in  Eqn.  (6-15) 
is  so  much  smaller  than  the  6  in  Eqn.  (6-26).  A^  =  1  at  z/L7  =  0, 
the  neutral  case.  The  implication  here  is  that  A^  will  be  within  5%  of 
unity  even  for  z/L'  =  -1  .  Since  the  controversy  about  the  magnitude 
of  constants  involves  indeterminancies  much  greater  than  5%,  the 
temporal  variability  of  e  certainly  exceeds  10%,  and  zq  has  local 
variations,  for  practical  purposes  well  into  the  free  convection  range 
we  can  probably  set  A^  =  1  .  Eqn.  (6-17),  the  neutral  approximation, 
is  deemed  adequate  in  non-neutral  cases;  it  may  be  as  suitable  as  any 
other  simple  equation. 


For  the  RMS  lateral  velocity  n  ,  or  the  RMS  angular  variations  of  a 

direction  vane  with  a  vertical  axis,  cl  »  the  situation  becomes  more 

<P 


complex  because  of  the  increasingly  large  energies  at  increasingly 
large  wavelengths  (or  periods)  which  are  suppressed  in  the  vertical 
turbulence  covered  above. 


For  the  neutral  case,  various  data  suggest  we  can  use  as  a  first 
approximation 


a  =  2  u* 
v 


(6-28) 


in  analogy  to  Eqn.  (6-21),  but  even  in  smooth  terrain  we  should  use 
a  large  zq  to  infer  an  "effective"  u*  such  as  1m  because  large  eddies 
dominate  the  quantity.  There  is  always  the  problem  that  oy  depends 
on  the  sampling  time  used  to  derive  it  --  if  longer  wavelengths  (times) 
are  included,  measured  ov  is  larger.  In  non-neutral  conditions  every¬ 
thing  gets  very  complex.  There  is  no  complete  theory  relating  to 
U,  zq,  R^,  primarily  because  the  instability  aloft  above  the  surface 
layer  can  have  a  large  effect  on  the  large  scale  lateral  motions  down 

low.  Table  III  tries  to  summarize  some  of  the  data  (note  0  ~  o  / U ), 

cp  v 

and  many  other  compilations  are  available  (see,  for  example,  Slade 
1968,  Pasquill,  1962,  and  Lumley  and  Panofsky  1964). 


To  summarize  from  the  last  reference:  "o  increases  with  in- 

v 

creasing  wind  speed,  at  constant  stability,  particularly  in  stable  air. 
It  further  is  generally  much  larger  in  unstable  air  for  the  same  wind 
speed,  with  the  exception  of  low- wind  speed  inversions  in  which  a 
gradual  drift  of  wind  direction  may  produce  large  standard  deviations. 


For  the  longitudinal  component,  0^  ,  in  anology  to  Eqns.  (6-21) 
and  (6-28),  we  find  in  the  neutral  case  o-  -u*,  with  the  constant  of 
proportionality  empirically  set  at  2.  5.  The  relative  variation  of  a 
with  stability  is  a  bit  less  than  given  for  ow  in  Table  III.  Panofsky 
(1970),  provides  the  latest  summary  plot  of  the  constant  of  proportion¬ 
ality  for  o  ~  u#  and  a  -  u*  ,  which  we  summarize  in  a  later  section. 
’  u  v 
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Spectra  and  Scales  of  Turbulence 

The  high  frequency  end  of  the  atmospheric  turbulence  horizontal 
spectrum,  the  inertial  subrange,  in  neutral  and  turbulent  conditions, 
fits  the  "-5/3  law": 


i 

S(k1)  =  be3 


(6-29) 


where  is  wave  number  in  radians  per  unit  length,  e  is  the  dissi¬ 
pation  rate,  S(k^)  is  the  energy  per  unit  wave  number,  and  b  is  0.  50 
for  longitudinal  spectra  and  0.  66  for  lateral  spectra.  If  k^  is  the  wave 
number  in  cycles  per  unit  length,  b  is  0.  15  for  longitudinal  spectra. 
(See  the  review  by  Panofsky,  1970,  and  Pasquill,  1972).  There  is  still 
controversy  about  the  value  of  constants.  Franzen  (1973)  considers 
b^  =  0.  55  as  the  best  estimate,  and  suggests  that  the  relation  between 
b^  and  von  Karman's  k  makes  the  latter  take  the  value  0.  35. 


With  the  shape  of  the  high  frequency  end  of  the  spectrum  known, 

it  is  possible  to  derive  the  intensity  factor,  e  ,  from  any  measurement 

which  depends  only  on  wavelengths  within  the  inertial  subrange.  For 

1/3 


from  the  Og  or  signals  where  only 


example,  we  can  obtain  e 

high  frequencies  (in  the  inertial  subrange)  are  used  in  the  "sigma 

meter"  giving  the  averages  of  8  or  cp  resulting  in  Oq  or  (here  we 

will  call  it  0  ).  The  filter  characteristic  time  T  must  be 

To  i 


T3<  0.44  z/U 


(6-30) 


(see  MacCready  and  Jex,  1964);  for  accuracy  it  is  safest  for  t  3  to  be 
even  considerably  shorter  than  Eqn.  (6-30).  Then  the  authors  find 


e>  =  0.62  U*  r7  o  (6-31) 

3  T3 
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in  degrees. 


2-3 

for  e  in  cm  sec  ,  U  in  m/sec,  T,  in  seconds,  and  a 

5  t3 

e1/?3  can  also  be  derived  from  high  frequency  wind  speed  fluctu¬ 
ations  (see  Franzen,  1965).  The  sensor  must  have  a  very  short  re¬ 
sponse  distance  to  permit  the  measurement. 

For  the  whole  spectrum  of  vertical  velocity  up  to  z  =  50m,  Pan- 
ofsky  (1970)  notes  for  neutral  and  moderately  unstable  cases  that  the 
spectra  can  be  normalized  to  the  following  equation  for  neutral  anu 
moderately  unstable  air: 


kIS(klj  _  3. 36f 

u*2  1+1  Of5/3 


(6-32) 


?  2. 

where  f  =  z  ,  and  u*  can  be  replaced  by  3.  6  .  Note  here  k^ 

is  in  terms  of  cycles  (not  radians)  per  unit  length,  and  f  is  a  frequency 
normalized  by  height. 


kjS(kj)  has  a  maximum  value  of  0.43  u*^  for  k^  z  ~  0.  3  ,  that  is, 
for  a  wavelength  about  3  times  the  height.  In  stable  air  the  peak  shifts 
to  higher  values. 


More  recently,  Pasquill  has  examined  the  S(kj)  spectrum  in 
much  greater  detail.  He  considers  the  scale  X  ,  the  wavelength  of  the 
peak  of  a  kjS(kj)  cur  e,  where  1^  is  the  wavelength  of  maximum  power 
input,  expressible  also  in  corresponding  wave  number  by  1^  ~  (Kjm) 
The  situation  can  be  summarized  as  follows: 


a.  In  neutral  flow  \  /z  is  between  2  and  4,  and  effectively 
constant  with  height  in  the  first  20m  or  so.  The  best  evidence  tends 
toward  the  value  2. 
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b.  The  effect  of  thermal  stratification  is  to  increase  or  de¬ 
crease  the  scale  in  unstable  or  stable  conditions,  and  in  effect  to  in¬ 
crease  or  decrease  the  height  range  with  effectively  linear  increase. 

Pasquill  also  discusses  the  Eulerian  length  scale,  £_  ,  derived 
from  correlation  coefficients.  Since  this  depends  rather  strongly  on 
the  large  wavelengths,  for  which  measurement  (averaging)  is  diffi¬ 
cult,  the  relative  variability  of  is  even  greater  than  that  for  X 

Summa'v  11 ;  above  holds  for  as  well  as 

£ 

In  final  summary,  for  maximum  accuracy  the  scale  constants 
in  Eqn.  (6-32)  should  be  adjusted  as  a  function  of  stability.  However, 
we  feel  that  for  our  practical  application  Eqn.  (6-32)  as  it  stands  is 
probably  adequate. 


The  Total  Boundary  Layer 

Hanna  (1969)  examined  many  different  methods  ior  estimating 
h.j.  ,  the  height  of  the  mixed  region  constituting  the  total  boundary 
layer.  He  found  it  impossible  to  relate  h^.  to  quantities  measured 
simply  near  the  surface,  such  as  u®,  U,  !the  geostrophie  wind), 
together  with  the  Coriolis  parameter.  The  best  formula  for 
larger  than  about  150m  was  found  So  be  that  ol  Laikhtman  (i961). 
with  the  constant  of  proportionality  altered  to  fit  observations: 

,  -1 

b^U^.aO')  »  0.75  Us(^~~)  *  (6- 33) 


where  is  the  average  vertical  gradient  of  potential  temperature 

through  the  boundary  layer.  Quoting  from  Hanna,  when  a  complete 
vertical  temperature  sounding  is  available,  "h^.  is  the  lowea*  level 
which  tie  vertical  gradient  of  temperature  exhibit*  a  discontinuity. 
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For  example,  during  well-mixed  afternoon  periods  when  an  adiabatic 
layer  near  the  ground  is  capped  by  a  relatively  stable  inversion  layer, 
h,p  corresponds  to  the  level  of  the  base  of  the  inversion  layer.  During 
a  clear,  calm  night-time  period  where  there  is  a  ground  based  invers¬ 
ion,  h^,  corresponds  to  the  top  of  the  strong  inversion  layer."  For  the 
altiludes  studied,  A9  1  was  always  positive  (conditions  always  stable 
over  this  layer). 

The  conclusion  is  that  temperature  measurements  aloft  are  needed 
to  establish  h^,  with  any  reliability.  With  such  measurements,  h^  can 
be  taken  directly  from  the  sounding.  Eqn.  (6-33)  may  be  of  some  slight 
help  if  a  partial  AQ ' / A  z  is  given. 


Empirical  Values  of  zq 


The  surface  boundary  layer  equations  are  based  on  a  picture  of 
flat  terrain  with  the  roughness  everywhere  the  same.  The  equations, 
however,  are  used  for  approximating  reality  in  actual  terrain,  which 
can  be  considered  as  an  assemblage  of  areas  each  with  a  different  re¬ 
presentative  zq  ,  all  being  averaged  together  to  give  a  final  net  effect. 
One  effect  of  a  change  of  roughness  can  be  seen  from  the  analysis  by 
Panofsky  (1968).  He  found  that  if  wind  blows  over  a  surface  whose 
roughness  length  abruptly  changes  to  zq  (at  x  =  0),  then  below  the 
height  Zj  where 


0.8 

=0.8  zq  (x/zq) 


(6-34) 


the  boundary  layer  similarity  laws  can  be  applied  based  on  the  new 

z  .  This  whole  subject  of  changing  and  patchy  roughness  is  exten- 
o 

sively  reviewed  by  Pasquiil  (1972).  Where  small  zq  is  involved,  it 
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takes  about  1  kui  for  il  e  influence  to  cover  the  height  of  concern  to 
us  (70m)  for  neutrai  conditions,  ar  J  much  more  for  stable  conditions. 
For  larger  zQ  ,  or  for  unstable  eases,  the  influence  establishes  it¬ 
self  more  quickly.  For  pitchy  roughness,  the  situation  is  obviously 
complex,  and  the  reader  is  referred  to  Pasqu-ll's  paper. 

Some  representative  values  of  roughness  are  given  below: 


Smooth  snow 

Mown  grass 

Long  grass  (60-70m) 

Open  Country  (O'Neill) 
Cotton  Field 
San  Jose,  suburb 
London,  urban 
Philadelphia 

Broekhaven  tower  (in  woods) 


»  .  005  to  0.  I 
c.  0 

6 

7,8 

SO 

fO 

78 

"90  to  240,  depending 
on  wind  direction 

100 


6.4  A  REVIEW  OF  OPERATIONAL  EQUATIONS 


Preliminary.  Find  z  ,  R.,  L',  u* 
— - .  ■  ■  ■  ■•  —  ■«-  o  \ 


Assume  one  measures  U  and  T  at  two  heights,  the  subscript 
O)  referring  to  the  lower  at  Zj  and  the  subscript  (2)  referring  to  the 
higher  at  •  U  and  T  designate  averages,  averaged  over  about  10 
minutes. 

1.  Tn  near  neutral  conditions,  preferaLly  with  strong  winds, 
calculate  z 

o 


From  Uj  at  z^  and  at  z^ 


u*  , 

U  =  y"  lD  z/zc 


(a) 


leads  to 


U1 

trin  via  Ei 


In  z  = 


v^tr-rr 


it) 


Second  methad: 


From  ctj  at  any  height 


In  e  =  In  2 
o 


(c) 


o  -  ov 


From  0  at  any  z  and  U 
w 


In  %  -  in  z 

o 


0.  5U 
a 


w 


Third  method: 


Obtain  e 


1/3 


from  U  and  m  (or  o  )  with  a  very  short 
8  cp  ’ 


filter  characteristic  T,  where  T,,  <  0.44  z/U  .  Then 

o  o 


e3  =  0.62  U3t7  a. 

O 


2. 


Obtain  z  from 
o 


fn  z  =  in  z  - 
o 


2 


e  3z3 


Use  all  three  methods,  and  note  agreement.  Also  note 
systematic  variation  of  with  azimuth  angle  cp  . 

Calculate  R.(z) 


0  t  ^  0  1 

R;(z)  =  (g/T)  — - l—z-z£n  z?/z. 

1  (U^Ujr  1 

which  gives  the  best  finite  difference  approximation  to 
the  formula  for  Richardson  Number  where  "z  =  (z^z^1 

is  the  geometric  mean  of  the  heights  and  z^  . 
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3.  Calculate  L'  at  a  using  the  appropriate  equation  shown 
on  Table  6-1. 

4.  With  this  L'  ,  and  using  the  i'fz/L')  relation  from  Fig¬ 
ure  6-lr  calculate  u*  from  the  height  data  using 

TJ  =  T"  z/z0  "  ^(z/L1)]  (h) 

Then  calculate  u*  from  the  height  2  data.  The  t^o  u* 
values  should  be  fairly  close.  Take  their  average  as  the 
best  estimate  of  u*  .  One  could  weight  the  calculation 
from  the  lower  height  more  strongly  since  the  z/L1  ef¬ 
fect  is  less  at  low  heights,  but  the  calculation  from  the 
higher  height  may  be  more  suitable  as  a  basis  for  extrapo¬ 
lation. 

U  vs  z 

1.  Eqn.  (h)  gives  U  vs  z  for  all  heights  in  the  surface 
boundary  layer. 


e  vs  z 

1.  Calculate  R.  for  all  the  heights  of  interest  from  z/L1 
from  the  Table  6-1  equations 

1  /3 

2.  At  these  heights,  calculate  e  by 


[1-1/3  Ri(l-18Rj)1/4]l/3 
p.i/^x/3  (l_18Ri)1/12 


(i) 
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3. 


As  an  operational  alternative,  ignore  the  R^  effects 
(which  tend  to  be  small)  and  just  assume 


k3z3 


which  is  the  equivalent  of 

=  i  uik’  . 

vF  JxTzTTz 

2  3  1'  o 


(j) 


(k) 


If  U  is  measured  at  z  ,  and  z  is  known,  Eqn.  (k)  gives 

1  1 /3  1  ° 

e  vs  z  .  If  only  c  is  measured  at  z^,  Eqn.  (j)  then 

shows 


e 


i 

3 


1 


(1) 


o  vs  z 

w _ 

1.  Use  =  (1.3-0.  3  z/L')u*  (m) 

0.  5U 

2.  For  neutral  cases,  thus  °w  =  Yi\zjz~  (n) 

'  o 

and  i6  constant  with  height. 

3.  Equivalently,  a  =  0. 5/in  z/z  (o) 

S  '  '  o 

4.  If  H  can  be  estimated  rather  than  R.  or  L1  ,  use  the 

l 

algorithm  provided  in  the  discussion  below  Eqn.  (6-26). 
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a  and  o  vs  z 
v  u 


1.  Consider  ~  u*  and  ~  u*  ,  with  the  constants  of  pro¬ 
portionality  being  functions  of  stability  as  summarized  on 
Figure  1  of  Panofsky  (1970).  The  constant  is  about  2.  5 

for  all  stabilities  for  o  ,  while  for  o  it  is  about  1.  5 

u  v 

for  stable  cases  and  exceeds  3  for  unstable.  The  rela¬ 
tions  for  any  given  case  are  found  to  be  poor  because 
the  energy  yielding  and  0^  is  not  closely  tied  to  local 
small-scale  zq  and  so  similarity  laws  are  really  inap¬ 
plicable.  If  accuracy  of  and  0^  is  required,  they 
should  be  measured  directly  at  one  point  and  extrapo¬ 
lated  upward  and  downward  by  physical  reasoning  and 
empirical  results  from  various  prior  studies.. 


Vertical  Turbulence  Spectra 

1.  For  neutral  and  moderately  unstable  cases, 

kl  S(kl)  _  3.  36 f 

u*2  1  +  1  Of^ 


(P) 


where  k  is  wave  number  in  radians  per  cm,  f  =  k  z  , 

1  2  * 

S(k.)  is  the  power  spectral  density,  and  u*  can  be 

*  2  2 
replaced  by  0.  6aw  ,  k^S(k^)  has  maximum  of  0,  43u* 

at  kjZ  ~0.  3  .  In  stable  air  this  peak  shifts  to  higher 

values  of  kj  . 


6.5  AN  APPROACH  TO  AN  OPERATIONAL  DATA  SYSTEM 


The  system  must  be  fully  automatic,  requiring  no  manual  inputs. 
For  eventual  operational  application  it  should  provide  information  on 
mean  air  flow  profiles  and  turbulence  throughout  a  "target"  volume 
about  1  km  long,  300m  wide,  and  of  thickness  tapering  from  ~  70m  at 
one  end  to  50m  at  the  other.  For  initial  research  purposes  where  one 
is  evaluating  monitoring  and  forecast  schemes,  the  volume  of  inter¬ 
est  can  be  a  bit  smaller  (especially  as  regards  length). 

For  the  system  which  is  to  aid  in  wake  forecasting,  the  atmos¬ 
pheric  characteristics  it  presents  are  estimates  of  future  conditions 
where  "future"  means  a  few  seconds  up  to  several  hundred  seconds. 
Conceptually,  the  forecast  could  be  rather  accurate  if  an  extensive 
observational  system  were  set  up  and  an  associated  extensive  compu¬ 
tational  system.  The  observational  setup  would  have  to  cover  all  the 
air  which,  for  the  period  of  prediction,  would  move  into  the  target 
volume.  For  a  strong  wind  of  20m/sec,  a  forecast  period  of  200 
seconds,  and  considering  the  vertical  and  horizontal  variations  of  the 
wind,  for  one  wind  direction  observations  as  far  as  2000m  from  the 
target  would  be  needed,  and  from  heights  up  to  about  500m  and  cross- 
wind  dimensions  of  1  km  or  so.  To  cover  all  wind  directions  and 
speeds,  obviously  a  dense  three-dimensional  network  extending  out 
several  km  in  all  directions  would  be  required.  Even  with  such  a 
network  and  an  elegant  computational  scheme,  small-scale  turbu¬ 
lence  could  not  be  forecast  in  detail. 

For  simplicity,  atmospheric  motions  are  usually  divided  into 
two  scales  of  motion:  (1)  mean  motion,  explicitly  described,  and 
(2)  turbulent  motion,  described  only  in  the  form  of  statistical  aver¬ 
ages.  Actually,  the  two  scales  form  parts  of  one  continuum.  For 
vortex-wake  transport  and  decay,  it  is  customary  to  consider  that 
the  air’s  mean  motion  provides  the  medium  for  specifying  drift, 
while  the  turbulence  determines  the  decay  mechanism  and  speed. 
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However,  the  whole  motion  of  the  wake  evolves  m  some  tens  of  sec¬ 
onds  or  a  hundred  seconds  or  so,  which  can  be  a  time  down  in  the 
turbulent  eddy  scale  rather  than  the  mean  scale.  To  be  specific,  in 
unstable  conditions  the  mean  vertical  wind  at  50m  over  a  flat  surface 
may  be  zero  and  yet  an  individual  vortex  m^y  spend  its  minute  of  life 
in  air  rising  at  2m/sec.  Similarly,  there  may  be  no  crossflow  com¬ 
ponent  of  vertical  shear  in  the  mean,  but  for  the  evolution  of  an  indi¬ 
vidual  vortex  system  the  shear  may  be  considerable. 

Such  considerations  force  one  into  some  important  practical 
conclusions  with  respect  to  forecasting  vortex-wake  characteristics. 

1.  Forecasting,  to  be  of  practical  value,  must  be  for  a  peri¬ 
od  of  10  minutes  or  more  --  enough  time  so  total  terminal  traffic  flow 
can  actually  be  adjusted  to  the  anticipated  conditions.  Since  there  is 
no  chance  of  forecasting  individual  eddies  over  such  a  time  interval, 
the  meteorological  forecast  is  limited  to  providing  only  mean' values  -- 
especially  the  most  important  limits  of  mean  flow  conditions  and  tur¬ 
bulence  statistics.  The  forecasting  data  base- must  be  at  least  10  min¬ 
utes,  and  experience  at  the  site  may  show  30  minutes  or  even  60  min¬ 
utes  to  be  better. 

2,  For  the  monitoring  of  conditions  in  the  ILS  localizer  beam 
near  the  middle  marker,  one  wants  to  establish  the  mean  flow,  turbu¬ 
lence,  and  thermal  stability  throughout  the  Lagrangian  parcel  of  air  in 
which  the  observed  vortex  moves  and  did  move.  The  prohibition  against 
research  aircraft  and  high  towers  in  this  region  means  direct  measure¬ 
ments  cannot  be  made.  Indirect  probing,  say  by  doppler  lidar  for  vel¬ 
ocity  and  passive  1R  for  temperature,  are  not  deemed  practical  at  the 
present  stage  of  development.  Thus  one  must  do  the  best  one  can  by 
extrapolations  from  surface  (or  low  tcwer)  measurements. 
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All  such  considerations  seem  to  force  the  design  of  an  operational 
system  along  the  following  lines.  The  suggested  system: 

a.  is  directed  especially  toward  the  evaluation  of  wake-vortex 
transport-decay  forecasts, 

b.  involves  some  redundancies  since  the  system,  applied  to  the 
specific  site,  utilizing  the  not-entirely- satisfactory  theoretical 
relations  of  the  surface  boundary  layer,  should  in  itself 
permit  crosschecks, 

c.  should  itself  receive  some  evaluation  and  calibration  by 
specific  pibal  releases  and  measurements  from  a  test  air¬ 
craft  moving  along  the  ILS  path,  and 

d.  is  more  complex  than  the  minimum  which  might  suffice  in 
future  years  for  operational  airport  uses. 

The  system  consists  of  a  network  of  tower-mounted  meteorolo¬ 
gical  sensors  with  appropriate  data  storage  and  real-time  computing 
and  display. 

Tower  #1  will  be  50m  high  (70m  would  be  even  better)  and  hope¬ 
fully  can  be  located  within  2  km  of  the  test  volume  where  upwind  sur¬ 
face  conditions  are  not  dissimilar  from  conditions  upwind  of  the  test 
volume.  An  existing  airport  control  tower  can  perhaps  be  used,  es¬ 
pecially  if  a  small  meteorological  tower  is  superimposed  on  it,  but 
one  must  take  account  of  the  disturbance  of  a  fat  tower  on  U  and  cp 
(azimuth)  measurements. 

At  z  =  10,  20,  35,  and  50m  there  should  be  measurements  of 
U,  cp  and  T  ,  from  which  10  minute  averages  are  derived. 
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Towers  #2,  3,  4,  and  5  will  each  be  just  10m  high,  instrumented 
only  at  the  top,  where  U,  Q,  and  a  (elevation  angle)  will  be  monitored. 
The  towers  should  be  put  at  the  corners  of  a  square  with  500m  sides 
straddling  the  landing  path  with  the  downwind  side  centered  on  the  mid¬ 
dle  marker.  The  measurements  will  be  processed  to  include  10  min¬ 
ute  averages  of  U,  £,  and  a  >  as  well  as  higher  frequency  data  to  be 
discussed  below. 

Desirable,  but  not  absolutely  essential,  measurements  would  be 
9  at  z  =  20m  and  z  =  5m  on  Tower  #1,  and  T  at  z  =  2m  and  10m  on 
Towers  #2-5.  A  monostatic  acoustic  sounder  giving  information  on 
low  altitude  layering  during  stable  conditions  would  also  be  helpful.  It 
would  be  located  someplace  within  the  square  of  Towers  #2-5. 

Tower  #1  is  essential  to  give  meaningful  stability  information,  for 
which  Towers  #2-5  are  presumed  low.  It  also  can  give  this  information 
aloft  (near  the  60m  height  of  the  glide  slope  over  the  middle  markers). 
This  stability  information  tends  to  characterize  conditions  ove-  a  large 
area,  hence  Tower  #1  need  not  be  right  at  the  test  volume.  In  strongly 
stable  conditions,  vertical  shears  aloft  in  speed  and  direction  can  also 
be  monitored  from  Tower  #1. 


^  T  measurements  on  all  the  low  towers  may  be  found  useful  if 
the  relative  sensor  accuracy  is  0.  1C  and  careful  averaging  over  all 
towers  is  done.  This  cannot  give  accurate  R.^  or  L1  measurements, 
but  can  at  least  categorize  the  more  extreme  stabilities  and  instabili¬ 
ties  where  the  neutral-case  equations  may  not  suffice. 


Tower  #2-5  are  used  first  to  derive  mean  r  characteristics  for 

o 

the  area.  Then,  for  operational  forecasts,  they  provide  U  and  Q  in¬ 
formation  from  which  meaningful  average  U  vs  Q  are  derived  to  use 

1  /3 

as  inputs  to  calculate  U  vs  z  ,  e  vs  z  ,  vs  z  ,  and  vertical  tur¬ 
bulence  spectra  vs  r.  .  From  the  high  frequency  ends  of  0  and  0 

8  — 


1/0 

(evaluated  by  sigma  meter)  one  also  can  calculate  c  '  at  10m.  and 
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from  the  low  frequency  sigma  meter  evaluation  of  Qq  one  derives 
at  JOm.  These  provide  useful  redundancy.  Upcurrents  and  downcur- 
rents  are  estimated  from  the  convergences  and  divergences  measured 
using  the  U  and  Q  values  of  Towers  #2-5.  For  the  area  defined  by 
Towers  #2-5,  instantaneous  components  of  wind  perpendicular  to  the 
sides  of  the  square  will  show  the  convergence-divergence  situation. 

A  similar  calculation  for  the  area  defined  by  any  three  of  the  towers 
can  give  clues  as  to  centering  of  the  up  or  down  current,  if  such  is 
desired. 

The  research/evaluation  system  suggested  here  is  presumably 
much  more  complex  than  an  eventual  operational  system  will  be. 

After  the  "calibration"  of  the  site  which  the  use  of  this  system  pro¬ 
vides,  a  final  operational  system  could  conceivably  be  as  simple  as 
(a)  a  single  anemometer  at  5m  or  10m,  (b)  one  measure  of  a  quantity 
related  to  heat  flux  --  a  temperature  between  2m  (or  ground)  and  I  0m  -- 
and  (c)  something  to  indicate  the  height  of  stable  layers,  such  as  an 
automated  short-range  acoustic  sounder.  For  the  critical  stable  cases 
some  direct  observacions  of  stability  and  wind  shear  from  a  tower  to  at 
least  50-60m  is  highly  desirable,  even  if  the  tower  is  a  few  miles  from 
the  operational  site.  The  acoustic  sounder  should  at  least  be  able  to 
show  the  height  above  which  turbulence  is  low,  and  thus  illuminate  cri¬ 
tical  hazard  conditions.  Radiation  measurements  are  so  poorly  coupled 
to  heat  flux  or  stability- instability  that  they  are  not  recommended. 

The  system  detailed  in  this  section  does  not  give  as  much  detail 
as  one  might  want  for  checking  the  accuracy  of  monitoring  system. 

Such  detail  could  only  come  from  indirect  probing,  as  with  a  doppler 
lidar  which  is  considered  not  yet  fully  operational,  or  from  a  network 
of  tall  towers,  which  is  inappropriate  as  long  as  the  tests  are  located 
at  an  airport.  Even  though  this  system  is  not  perfect  lor  case  studies 
of  monitoring  devices,  the  system  will  likely  be  found  satisfactory  for 
that  purpose. 


7.  SUMMARY 


The  general  factors  relating  aircraft  wing  geometry  to  the  vorticity 
shed  from  the  wing  are  presented.  The  rollup  of  the  resulting  vortex 
sheet  into  the  classical  wake  consisting  of  a  pair  of  contra-rotating 
trailing  vortices  is  described,  and  a  set  of  operational  equations  for 
principal  wake  characteristics  is  presented. 

The  wake  encounter  hazard  with  respect  to  a  following  aircraft  is 
discussed.  Axial  a  d  normal  encounter  danger  factors  are  developed. 
These  factors  relate  the  aerodynamic  response  and  control  characteristics 
of  the  following  aircraft  to  the  wake  characteristics  of  the  generating  air¬ 
craft.  It  is  recommended  that  the  axial  danger  factor  D  be  used  as  a 
primary  criterion  for  determining  hazav.  ious  encounter  situations. 

The  general  conditions  relating  to  the  transport  of  a  vortex  pair  are 
reviewed.  The  effects  on  the  vortex  pair  dynamics  du>*  to  vorticity  gen:  ra¬ 
tion  and  re-distribution  are  discussed.  Vorticity  generation  is  caused 
by  density  variation  in  the  fluid.  Vorticity  re-distribution  results  from 
the  ’iffusing  vertex  cores  and  the  entrainment  of  mass  and  the  detrain- 
meat  of  vorticity  at  the  boundaries  of  the  recirculation  cells  surround¬ 
ing  each  vertex.  In  general,  this  diffusion  of  the  vortex  cores  results  in 
departure  from  the  classical  equations  describing  the  . notion  of  the  vortex 
pair.  Hf’sver  it  is  pointed  out  that  this  ;.s  a  relatively  slow  process,  and 
usually,  before  it  becomes  significant,  the  vortex  pair  has  experienced  oi?“ 
of  several  decay  mechanisms  (Unking,  core  bursting). 

The  effects  of  (croaswind)  wind  shea;  are  examined.  In  this  case  the 
geometry  of  the  vortex  recirculation  cells  is  drastically  altered.  The  up¬ 
wind  cell  increases  in  size  and  the  downwind  cell  shrinks.  This  effect  in¬ 
creases  with  shear,  .  nd  th  psmbable  effects  that  wind  shear  might  have  on 
banking  of  the  vortex  pais'  and  core  bursting  are  discussed. 
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Previous  experimental  and  theoretical  works  in  the  area  of  wake 
descent  into  a  stratified  atmosphere  are  discussed.  The  conclusion 
is  that  at  the  present  there  seems  to  be  no  model  that  adequately  de¬ 
scribes  the  phenomenon  of  vorticity  generation.  A  new  theory  for  the 
descent  of  a  wake  into  a  weakly- stratified  atmosphere  is  developed.  For 
the  condition  of  weak  stratification  (large  Brunt  Vaisala  times)  and  fast 
wake  descents  (time  for  wake  to  descend  one  vortex  span)  —  a  situation 
which  normally  occurs  in  practice  —  the  theory  concludes  that  the  vortices 
will  draw  closer  together  with  a  consequent  increase  in  descent  speed. 

Thus  the  effect  of  stable  stratification  alone  (without  entrainment)  is  to 
increase  descent  speed. 

A  discussion  of  the  effects  of  entrainment  alone  (without  buoyancy' 
is  given.  It  appears  that  entrainment  will  definitely  retard  wake  descent, 
however  the  magnitude  and  mechanism  of  entrainment  is  controversial. 

The  general  concepts  regarding  Crow  Instability  are  discussed.  A 
new  theory  is  developed  as  an  extension  of  the  Crow  Stability  Theory 
which  assumes  an  atmospheric  turbulent  input  to  the  growth  of  the  in¬ 
stability.  The  /ortices  are  treated  as  a  pair  of  linear  oscillators  with 
a  random  forcing  function  assumed  to  be  that  corresponding  to  the  on e- 
dimensional,  transverse  energy  spectrum  within  the  Kamolgorov  inertial 
subrange,  The  "time-to-linking"  is  found  from  the  new  theory  as  a  func¬ 
tion  of  the  vortex  parameters,  F^  and  b  ,  and  tne.  turbulence  dissipation 
rate,  e  •  This  correlates  very  well  with  available  flight  test  data. 

Core  bursting  is  discussed  with  respect  to  global  invariant#  of  core 
development.  The  nature  and  precise  mechanism  of  vortex  breakdown 
is  still  controversial.  It  is  generally  agreed  that  breakdown  is  always 
associated  with  adverse  pressure  gradients.  The  conditions  for  break¬ 
down  are  related  principally  to  the  swirl  ratio  and  the  magnitude  of  the 
axial  flow,  botn  parameters  being  a  function  of  the  developing  core  flow. 
The  exact  nature  of  the  dissipative  processes  in  influencing  the  develop- 


ment  of  the  core  is  still  not  known,  but  plays  an  important  role  in  the 
core  flow,  and  hence  in  the  time  required  for  breakdown. 

The  final  section  shows  how  to  determine  the  atmospheric  dynamic 
parameters  from  general  meteorological  information,  using  such  con¬ 
cepts  as  Kolmogorov's  hypothesis  and  Monin-Obukhov  Similarity. 
Operation..!  equations  are  presented,  and  a  projected  design  for  an 
operational  meteorological  data  system  is  given. 


8.  CONCLUSIONS 


The  general  thrust  of  this  report  is  concerned  with  the  still  un¬ 
solved  problems  of  vortex  descent  and  decay  in  a  real  atmosphere.  Con¬ 
siderable  theoretical  progress  has  been  made  in  certain  aspects  of  this 
problem,  in  that  new  analytical  results  have  been  developed  for 

1)  Time-to-linkage  in  a  turbulent  environment. 

2)  Descent  in  a  stratified  flow. 

3)  Cell  shapes  in  sheared  flow. 

These  theoretical  results  can  only  be  validated  by  experiment,  and 
a  major  recommendation  is  that  the  appropriate  data  be  obtained.  The 
analyses  developed  here  and  the  specific  problems  identified,  will  provide 
the  proper  basis  for  sound  experiments  designed  to  resolve  some  of  these 
difficulties. 

For  the  important  area  of  core  bursting,  the  conclusion  is  that 
still  no  definite  operational  criteria  exist  and  that  more  experimental  and 
theoretical  research  is  required.  It  is  possible,  however,  that  natural 
core  bursting  seldom  occurs  in  the  vortex  wakes  ox  large  aircraft  and 
thus  an  operational  criteria  is  not  a  practical  necessity.  Further  flight 
test,  and  examination  of  existing  test  data  is  necessary  to  establish  this 
hypotheei ». 

The  problem  of  the  solitary  vortex,  and  the  possibility  that  this  is 
caused  by  a  sheared  crosswind  is  raised  in  this  report.  This  appears 
to  be  an  important  operational  conside  ration,  more  experimental  data 
here  is  badly  needed.  It  is  likely  that  further  analysis,  continuing  the 
sheared  cell -shape  developed  in  this  report,  would  assist  in  unde  rsiandsng 
ibis  phenomenon. 
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The  turbulent  entrainment  in  the  descending  vortex  pair  is  discussed 
at  length.  Here  again,  experiments  are  needed,  and  the  report  defines 
appropriate  parameters  which  would  be  essential  in  designing  the  proper 
tests  to  unambiguously  determine  the  turbulent  mechanisms. 

For  the  wing  aerodynamics  defining  the  vortex  configuration,  the 
aircraft  dynamics  defining  the  hazard,  and  the  meteorology  defining  the 
atmospheric  dynamics,  this  report  gives  specific  operational  equations 
which  can  be  used  in  predictive  models.  It  is  believed  that  in  these 
fields  the  state  of  the  art  is  sufficiently  well  developed  that  engineering 
calculations  can  be  made  with  reasonable  confidence. 

An  assessment  of  the  situation  with  respect  to  the  prediction  of 
vortex  wake  descent  and  decay  is  as  follows: 

1)  Wing  aerodynamics  defining  vortex  configuration,  aircraft  . 
dynamics  defining  hazard,  and  meteorology  defining  atmospheric  dynamics 
are  reasonably  well  understood. 

2)  Vortex  descent  and  decay  in  a  stratified  flow  is  still  not 
properly  understood. 

3)  Crow  Instab  lity  is  well  understood,  and  equations  presented 
here  give  a  rational  approach  to  predicting  time -to-linkage  as  a  function  of 
turbulence. 

4)  Core -bursting  is  still  very  poorly  understood  and  qualitative 
equations  are  lacking,  however  core  bursting  may  not  be  as  important  as 
Crow  Instability  for  vortices  from  large  aircraft. 

5)  Unsymmetric  effects  causing  vortex  tilting  and  unsymmetric 
breakdown,  resulting  in  a  long-lived  solitary  vortex  have  been  observed; 
and  are  not  understood.  These  may  be  very  important  for  hazard  pre¬ 
diction. 
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6)  Further  analytical  research  and  flight  test  experimentation 
is  very  definitely  required  to  resolve  the  problems  above,  both  at  altitude 
and  in  ground  effect.  This  report  defines  a  basis  for  the  rational  and 
effective  design  of  continued  research. 
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